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SUMMARY 


Project  Themis,  A  Center  for  Remote  Sensing  at  the  University  of 

Kansas  was  initiated  in  September,  1967  as  a  broad  interdisciplinary  effort 

for  the  improvement  of  user  utility  of  remotely  sensed  data.  A  balanced 

program  involving  theoretical  work,  sensor  development,  prc^e&sing  and 

display,  and  data  analysis  for  specific  user  application  was  e^^:ablished. 

During  this  report  period  the  program  has  been  implemented  by  work  on 

broad- spectrum  radar- target  interactions,  on  multi- image  processing  and 

* 

display  systems  (expansion  of  the  IDECS),  and  on  interpretation  of  radar 
images  for  geographic  and  geologic  purposes.  The  desi  n  of  a  4-8GHz 
continuously  swept  radar  to  be  mounted  on  a  "cherry  picker"  truck  is 
almost  complete  and  many  components  are  now  on  hand.  Acoustic  simula¬ 
tion  of  a  5- octave  system  has  'Started  with  equipment  assembled  during 
the  reporting  period.  Design  of  interfaces  between  the  existing  analog 
IDECS  and  digital  equipment  is  well  along,  with  some  units  already  under 
construction  and  most  equipment  ordered.  Adaptive  multi- image  category 
recognition  programs  have  been  tested  on  multi- polarization  radar  data 
from  an  agricultural  area,  and  the  results  look  quite  promising.  Prelimin¬ 
ary  work  on  use  of  radar  to  determine  the  geomorphology  of  the  Wasatch 
Range  in  Utah  is  complete.  A  geologic  and  geomorphologic  study  of  Darien 
Province,  Panama,  based  on  radar  imagery,  is  in  the  field  checking  stage. 

Project  Themis  at  the  University  of  Kansas  is  funded  under  ARPA- DOD 
Work  Order  1079  and  is  monitored  as  Contract  DAAK02- 68-C-0089  by  the 
U.S.  Army  Engineering  Topographic  Laboratories,  the  Geographic  Systems 
Division,  Bernard  B.  Scheps,  Project  Engineer.  The  project  is  supported 
at  the  University  of  Kansas  by  the  following  personnel,  some  of  whom  are 
not  directly  funded  under  Themis,  but  whose  work  is  contributory  to  the 
effort. 

*  Acronym  for  Image  Discrimination,  Enhancement  and  Combination  system, 
an  electronic  scanning  device  which  permits  simultaneous  viewing  and 
combination  of  multiple  images,  with  provision  for  signal  processing 
for  optimal  presentation  of  selected  terrain  phenomena. 


A.  Faculty  Members 


Electrical  Engineering  5 

Geography  1 

Geology  2 

B.  Ph.D.  Candidates 

Electrical  Engineering  14 

Geography  5 

Geology  1 

C.  Masters  Candidates 

Electrical  Engineering  7 

Geography  0 

Geology  3 


In  addition  to  the  personnel  listed,  a  number  of  undergraduate  students  are 
employed  as  technicians  under  the  three  principal  disciplines.  A  working 
arrangement  is  being  developed  with  faculty  members  of  the  departments  of 
Psychology  and  Sociology  in  the  field  of  pattern  recognition  and  false-color 
perception  and  will  undoubtedly  be  formalized  in  the  near  future  as  areas 
of  common  interest  become  more  clearly  defined. 

Since  its  inception  in  1964,  the  Remote  Sensing  Laboratory  at  the 
University  of  Kansas  has  been  actively  engaged  in  a  number  of  studies 
which  relavjj  directly  to  the  work  being  done  under  Project  Themis .  Themis 
is  of  particular  importance  to  the  University  of  Kansas  since,  because  of 
its  broad  objectives,  it  provides  a  focus  for  the  work  being  done  under  the 
other  studies.  As  a  result,  considerable  cross- fertilization  between  pro¬ 
jects  is  occuring  almost  automatically  and  the  total  effort  of  the  Remote 
Sensing  Laboratory  is  emerging  as  a  conesive  entity. 

The  studies  currently  underway  at  the  Remote  Sensing  Laboratory 
in  addition  to  Project  Themis  are: 
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PROJECT  THEMIS:  A  CENTER  FOR  REMOTE  SENSING 
SEMI  ANNUAL  TECHNICAL  REPORT  (U) 


I.  INTRODUCTION 

This  report  summarizes  the  technical  work  which  has  been  accom¬ 
plished  during  the  first  six -month  period  of  contract  no.  DAAK02-68-C-0089 . 
Technical  reports  and  memoranda  which  describe  the  work  to  date  are  refer¬ 
enced  in  the  text  and  are  appended.  The  appended  documents,  for  the  most 
part,  are  working  documents  which  will  be  up-dated  as  the  investigations 
continue.  When  sufficient  progress  has  been  made  that  significant  conclu¬ 
sions  may  be  drawn,  these  working  memoranda  will  be  revised  into  a  form 
suitable  for  scientific  publication. 

The  technical  work  is  discussed  under  the  task  headings  used  and 
defined  in  the  contract  purchase  description  and  further  defined  in  the  Study 
Plan  which  was  issued  as  a  report  by  the  University  of  Kansas  in  January, 

1968.  This  Study  Plan  deviates  somewhat  from  the  purchase  description  in 
subtask  definition.  The  reorientation  of  some  of  the  tasks  was  necessary  in 
order  to  most  efficiently  utilize  the  capabilities  and  interests  of  the  faculty 
members  and  graduate  students  who  are  performing  the  investigations.  Another 
factor  which  was  considered  in  task  reorientation  was  the  low  level  of  the 
manpower  budget  relative  to  the  level  anticipated  when  the  original  proposal 
for  contract  was  written.  Although  both  the  Study  Plan  and  a  subsequent 
proposal  for  renewal  indicate  some  re.  •  ientation,  the  overall  purpose  and  goals 
of  Project  Themis  at  the  University  of  Kansas  remain  unchanged. 


II.  TECHNICAL  REPORT 

TASK  4.1*  TARGET-SENSOR  INTERACTION  STUDIES 

The  emphasis  in  this  task  is  on  theoretical,  modeling  and  experimen¬ 
tal  studies  of  broad -spectrum  radar  systems.  These  studies  encompass  both 
scattering  measurement  and  imaging  systems. 


*  Task  numbers  are  those  of  the  USAETL  Purchase  Description 
dated  16  September  1967. 
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SUBTASK  4.1.1  THEORETICAL  STUDIES 

At  this  time  theoretical  scattering  studies  are  being  performed  at 
the  University  of  Kansas  under  National  Science  Foundation  grants.*  Work 
done  in  these  studies  which  is  applicable  to  broad  spectrum  scattering 
will  be  reported  both  to  the  granting  agency  and  to  USAETL.  The  work 
reported  to  USAETL  will  be  in  the  form  of  substantiation  of  empirical  model¬ 
ing  studies  which  are  being  conducted  under  task  4.1.2. 

SUBTASK  4.1.2  ULTRASONIC  MODEL  STUDIES 

The  equipment  for  ultrasonic  modeling  of  broad-band  radar  systems 
is  assembled  and  calibrated,  and  data  are  being  collected  from  a  variety  of 
targets.  The  system  description  is  appended  to  this  report  as  Appendix  A 
(CRES  Technical  Memorandum  133-2).  This  memorandum  is  one  chapter  of 
an  extensive  report  on  the  system  design,  data  acquisition  and  data  processing 
and  analysis.  The  report,  in  the  form  of  a  Ph.D.  dissertation  by  Mr.  John 
Rouse,  will  be  issued  during  the  summer,  1968. 

The  principal  difficulty  in  data  acquisition  at  this  time  arises  from 
the  broad  frequency  spectrum  (0.1  MHz  to  3,0  MHz)  that  is  being  investi¬ 
gated.  A  single  transducer  cannot  be  used  across  this  spectrum  due  to 
bandwidth  limitations,  so  it  is  necessary  to  use  a  number  of  transducers 
whose  spectra  overlap.  Thus*  Q  continuum  of  specti'um  is  not  possible  in 
the  experimental  setup  and  it  is  necessary  to  sweep  over  a  portion  of  the 
spectrum,  change  the  transducer  and  sweep  the  next  portion.  The  indivi¬ 
dual  transducer  characteristics  must  be  well  known  in  order  that  data  may 
be  corrected  to  produce  a  continuum  of  data  across  the  desired  band. 

SUBTASKS  4.1.3,  4.1.4,  4.1.5,  4.1.6  ELECTROMAGNETIC  INSTRUMENTA¬ 
TION  AND  MEASUREMENT 

These  subtasks  are  combined  into  a  single  task  description  because 
of  the  close  interrelation  of  subtasks  which  prevents  convenient  separation. 


*  NSF  Grants  GK-875  and  GK-1153 


3 


A  simple  octave-bandwidth  radar  system  is  being  constructed  at  the 
University  of  Kansas  which  will  operate  either  as  fixed  or  mobile  equipment, 
and  in  either  a  scattering -measurement  or  imaging  mode.  All  components 
have  been  specified  and  ordered,  and  construction  will  begin  in  the  near 
future  with  fixed  site  operation  planned  for  early  summer,  1968.  This  acti¬ 
vity  is  presently  somewhat  behind  schedule  as  projected  in  the  study  plan. 

A  GFE  boom  truck  has  been  received,  which,  while  in  poor  mechani¬ 
cal  condition,  can  be  repaired  to  be  functional  for  the  early  phases  of 
ground  mobile  operation.  Expert  mechanic's  opinion  advises  that  the  *ruck 
not  be  driven  too  far  from  a  repair  shop.  Effort  is  being  directed  in  parallel 
with  radar  system  construction  toward  fixing  antenna  mounts  atop  the  mov¬ 
able  boom  for  the  mobile  ground  measurements .  The  truck  should  be  opera¬ 
tional  and  the  antenna  mounts  in  place  as  scheduled,  although  some  time 
slippage  may  occur  if  repairs  must  be  made  to  the  boom  hydraulic  operating 
system . 

GRES  Technical  Memorandum  133-1-1,  again  the  first  portion  of  a 
major  report,  is  appended  as  Appendix  B.  This  memorandum  outlines  system 
parameters  in  terms  of  desired  performance  and  of  components  which  are 
readily  available  from  commercial  sources . 


TASK  4.2  SENSOR  AND  PREPROCESSOR  DEVELOPMENT 

This  task  includes  parameter  and  design  considerations  and  develop¬ 
ment  of  prototype  sensors  and  sensor  systems  including  associated  prepro¬ 
cessing  of  sensor  data.  Due  to  fund  limitations,  subtasks  4.2.2  (Carto¬ 
graphic  radar  systems),  4.2.3  (imaging  radars  for  small  satellites)  and 
4.2.5  (semi-focussed  imaging  systems)  are  not  being  pursued  under  this 
contract,  although  preliminary  studies  of  radars  for  small  satellites  are 
being  conducted  under  another  contract.*  Subtask  4.2,1  (evaluation  of 
broad-spectrum  imaging  systems)  will  be  commenced  when  data  from  the 
system  under  development  are  available. 


*  NASA  Contract  NAS9-7175 
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Subtask  4.2.4,  a  study  of  radar  scatterometer  systems  is  presently 
being  implemented  jointly  under  the  auspices  of  the  Naval  Ordinance  Labor¬ 
atories  (contract  #N  62306-67-C-0044)  and  Project  THEMIS  .  Parameter 
studies  are  being  conducted  for  sea -ice  measurements  and  include  consider¬ 
ations  of  roughness,  frequency,  resolution  area,  polarization  and  antenna 
patterns.  The  studies  are  sufficiently  general  in  nature  that  they  can  be 
extrapolated  to  terrain  phenomena.  A  working  memorandum  on  these  studies 
will  be  issued  in  June,  1968.  Included  will  be  preliminary  results  of  a 
comparison  study  between  pulsed,  cw  and  FM  systems. 


TASK  4 . 3  DATA  PROCESSING  AND  DISPLAY 

The  purpose  of  this  task  is  to  improve  existing  techniques  for  data 
display,  such  as  those  employed  in  the  IDECS  system,  and  to  originate  new 
techniques  for  processing  and  display  of  remote  sensor  data  for  improved  inter- 
pretibility.  At  this  time,  the  working  subtasks  under  4.3  fall  into  two  cate¬ 
gories,  1)  improvement  of  the  IDECS  system  (Subtasks  4.3.1  and  4.3.3)  and 
2)  implementation  of  adaptive  and  Bayes  decision  techniques  (Subtasks  4.3.4 
and  4.3,5).  As  work  progresses  in  these  categories,  subtasks  4.3.2  (evalua¬ 
tion  of  color  combinations  of  multiple  images)  and  4.3.6  (new  display  techniques) 
will  be  implemented.  Ultimately,  the  adaptive  and  Bayes  decision  techniques 
will  be  instrumented  and  interfaced  with  the  IDECS  system  in  such  a  way  that 
the  IDECS  operator  will  have  maximum  flexibility  in  enhancing  data  by  direct 
combinations  of  images  along  with  computer  processed  data. 

Representative  of  the  work  accomplished  to  date  under  task  4.3  are 
three  memoranda  which  are  appended  as  Appendices  C,  D  and  E.  The  first 
of  these  deals  with  the  mechanization  of  the  IDECS  system  improvement  and 
is  essentially  the  outline  of  a  working  technical  report  which  will  be  issued 
when  the  work  underway  nev^rs  completion.  The  second  and  third  memoranda 
describe  applications  of  adaptive  and  Bayes  decision  techniques  to  actual 
radar  data  obtained  from  an  agricultural  tost  site. 
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TASK  4.4  DATA  ANALYSIS  AND  APPLICATION 

Although  this  task  specified  analyst  interaction  with  tasks  4.1  and 
4.3,  to  date  the  only  interaction  which  has  occured  is  with  task  4.3,  data 
processing  and  display,  since  no  data  are  yet  available  from  the  broad- 
spectrum  radar  systems  of  task  4.1.  A  review  of  the  subtasks  as  outlined 
in  the  study  plan  follows: 

SUBTASK  4.4.1  EVALUATE  GEOGRAPHIC  AND  GEOLOGIC  INFORMATION 

AVAILABLE  FROM  SINGLE  RADAR  IMAGES 

A  preliminary  study  of  the  geomorphology  of  the  Wasatch  Range, 

Utah,  using  single  radar  images  is  appended  as  Appendix  F.  Originally 
the  study  plan  included  single  image  evaluation  not  only  of  general  sites 
but  also  of  a  specific  vegetation  site  and  the  geology  between  Topeka  and 
Olathe,  Kansas  as  part  of  Subtask  4.4.1.  By  agreement  with  Mr.  B.  Scheps, 
Project  Engineer,  this  portion  of  the  study  plan  has  been  redirected  to  field 
work  ip.  and  the  interpretation  of  radar  imagery  of  Darien  Province,  Panama. 

This  work  has  come  at  a  most  propitious  time  since  the  University 
of  Kansas  was  able  to  assign  to  it  two  graduate  students  who  have  had  pre¬ 
vious  experience  in  interpreting  imagery  from  Darien  Province  and  thus  were 
able  to  take  into  the  field  some  a  priori  knowledge  of  the  area.  Because  the 
work  was  performed  in  conjunction  with  a  "one-time"  data  acquisition  effort 
by  the  Engineering  Topographic  Laboratories,  it  was  imperative  that  support 
be  given  by  the  University  of  Kansas  at  this  particular  time.  The  particular 
qualifications  of  the  involved  University  of  Kansas  personnel  should  assist 
significantly  in  the  understanding  of  the  geoscience  aspects  of  this  important 
area.  A  preliminary  report  of  this  investigation  will  be  issued  during  the  summer, 
1968. 

Subtasks  4.4.2  (relation  of  image  data  to  wavelength,  polarization, 
depression  angle,  resolution  and  bandwidth)  and  4.4.3  (information  increase 
with  the  use  of  multiple  radar  images)  are  not  to  be  implemented  until  broad- 
spectrum  image  data  are  available  and  the  IDECS  system  is  fully  operational, 
respectively.  It  is  expected  that  the  image  data  will  be  available  during  the 
second  year  and  that  the  IDECS  system  will  be  operational  to  the  extent 
described  in  task  4.3,  Appendix  C,  during  the  summer,  1968. 
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It  should  be  noted  that  the  two  memoranda  listed  under  task  4.3 
(Appendices  D  and  E)  required  initial  geoscience  evaluation  of  both  single  and 
multiple  image  data  and  thus  should  properly  be  included  as  reports  under 
task  4.4. 


III.  CONCLUSIONS  AND  RECOMMENDATIONS 

Progress  has  been  made  in  all  phases  of  Project  Themis /  with  Study 
Plan  target  dates,  in  general,  being  met.  Some  delay  has  been  encountered 
in  task  4.1,  Target-Sensor  Interaction  Studies,  in  the  construction  of  the 
broad- spectrum  radar  system.  At  the  conclusion  of  the  next  report  period, 
it  is  anticipated  that  this  portion  of  task  4.1  will  bo  approximately  60  days 
behind  schedule.  The  delay  is  the  result  of  difficulties  encountered  in 
obtaining  components.  Because  of  the  heavy  demands  being  made  on  manu¬ 
facturers  and  suppliers  to  support  the  present  military  effort,  stock  items 
are  in  short  supply  and  numerous  delays  and  back  orders  have  been  encoun¬ 
tered.  Nearly  all  components  have  been  obtained  at  this  time,  however, 

and  construction  has  been  commenced. 

Since  work  is  progressLig  approximately  at  the  rate  anticipated  by 
the  Study  Plan,  it  must  be  concluded  that  the  plan  is  a  feasible  one  and 
that  it  is  within  the  constraints  of  both  funding  and  manpower.  It  is  recom¬ 
mended  that  work  be  continued  at  the  present  level  with  no  foreseeable 
redirection  during  the  next  report  period. 
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BROAD  SPECTRUM  BACKSCATTER  ACOUSTIC  SYSTEM 

J.  W.  Rouse,  Jr. 

INTRODUCTION 

Many  researchers  have  for  some  lime  recognized  that  the  mono¬ 
chromatic  nature  of  radar  hampered  its  performance  of  several  specific 
tasks  and  restricted  the  expansion  of  radar  techniques  into  new  fields  of 
remote  sensing.  Consequently  the  development  of  panchromatic  radar 
systems  is  underway.  The  panchromatic  or  broad -spectrum  radar  system  design 
poses  many  problems .  A  major  difficulty  arises  from  the  fact  that  the 
frequency  dependence  of  backscatter  is  for  the  most  part  undefined.  Such 
questions  as  the  appropriate  center  frequency  and  spectrum  width  for  any 
specific  application  can  only  be  answered  intuitively  based  on  a  compila¬ 
tion  of  data  from  multiple -discrete  frequency  measurements,  which,  for 
the  most  part,  is  an  educated  guess. 

Obtaining  answers  as  to  necessary  system  parameters  for  specific 
applications  is  a  major  concern  to  broad-spectrum  radar  designers.  Even 
aside  from  the  issue  of  system  design,  the  lack  of  definition  of  the  frequency 
dependence  of  backscatter  is  an  issue  of  concern.  The  great  bulk  of  theoreti¬ 
cal  descriptions  of  scattering  from  randomly  rough  surfaces  have  depended 
upon  experimental  results  to  justify  certain  simplifying  assumptions.  In 
addition,  of  course,  the  validity  of  the  individual  works  is  determined  by 
its  conformity  to  experimental  results.  In  the  case  of  angular  dependence, 
considerable  data  is  available  and  consequently  many  theories  lead  to 
excellent  correlation  with  experiment.  Unfortunately  no  such  store  of  con¬ 
trolled  experimental  results  exists  for  frequency  dependence.  The  limited 
number  of  measurements  available  are  either  too  widely  spaced  or  the  terrain 
(target)  is  too  inadequately  defined  to  be  of  much  value. 

The  primary  objectives  of  this  work  are  to  obtain  an  experimental 
description  and  an  analytical  expression  for  the  frequency  dependence  of 
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backscatter  from  randomly  rough  surfaces  .  The  secondary  objectives  include 
providing  experimental  verification  of  theories  predicting  the  frequency  de¬ 
pendence  of  backscatter,  and  examining  the  broad -wpectrum  technique  applied 
to  traditional  radar  problems  such  as  scintillation,  glint  and  the  resonance 
phenomenon. 

The  entire  experiment  program  is  conducted  using  acoustic  waves  in 
water.  It  is  important  that  the  reader  appreciate  the  analogy  between  acoustic 
waves  and  electromagnetic  waves,  and  understands  that  the  results  obtained 
by  this  research  are  applicable  to  broad-spectrum  radar  systems .  Specific 
cases  where  the  electromagnetic  wave  and  acoustic  wave  analogy  is  imperfect 
are  pointed  out  in  the  text,  however  in  general  one  can  assume  that  the  results 
obtained  by  measurements  of  acoustic  wave  returns  from,  for  example,  a  sand 
coated  target  in  a  water  tank  may  just  as  well  have  been  obtained  by  flying 
a  radar  over  a  particular  terrain  region  having  the  same  statistical  properties. 
Polarization  effects  cannot  be  determined  since  the  incident  acoustic  waves 
are  scalar  in  nature. 

Acoustic  simulation  of  electromagnetic  problems  is  well  developed  and 
documented.  Even  so  its  value  as  an  analog  computer  is  either  not  understood 
or  not  appreciated  by  many  radar  engineers .  An  acoustic  simulator  makes  it 
possible  to  perform  many  of  the  experiments  normally  reserved  for  flight  testing 
in  the  laboratory,  and  to  adjust  parameters  of  the  radar  system  in  a  highly 
controlled  and  repeatable  manner.  The  simulator  is  advantageous  because 
both  distance  and  time  are  scaled  by  a  factor  containing  the  ratio  of  acoustic 
propagation  in  water  to  the  electromagnetic  propagation  in  air;  a  factor  of 

2  X  10^.  Therefore  the  radar  problem  is  compressed  to  a  size  easily  accom¬ 
modated  in  a  water  tank,  and  the  electronic  system  is  video  rather  than  rf . 

The  flexibility  is.  further  enhanced  by  the  option  to  scale  distance  in  a  non¬ 
linear  manner.  The  essence  is  that  model  experiments  can  be  performed 
rapidly,  repeatably,  and  at  much  smaller  cost  than  full  scale  experiments,  and 
can  provide  results  completely  analogous  to  a  large  number  of  electromagnetic 
problems . 

This  research  program  is  subdivided  into  four  action  areas:  1)  devel¬ 
opment  of  instrumentation  to  implement  the  broad-spectrum  technique. 
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2)  mpasurement  of  backscatter  from  randomly  rough  targets  of  known  statis- 
tiC'i,  3)  determination  of  the  frequency  dependence  predicted  by  several 
scattering  theories,  and  4)  correlation  of  experimental  and  theoretical  results. 
In  addition,  emphasis  is  placed  on  the  effects  associated  with  implementing 
the  broad-spectrum  technique  in  the  microwave  region.  For  example,  images 
of  a  variety  of  targets  illuminated  with  a  broad -spectrum  signal  will  be  pre¬ 
sented  and  discussed  in  a  later  report. 

The  design  and  construction  of  a  system  providing  broad-spectrum 
illumination  occupied  much  of  the  effort  of  this  program.  Of  particular  con¬ 
cern  is  the  behavior  of  the  ultrasonic  transducers  used  as  the  transmitter 
and  receiver  antennas.  These  devices  are  piezoelectric  crystals  constructed 
for  resonance  operation  at  particular  discrete  frequencies.  Operation  of  these 
devices  off-resonance,  as  is  necessary  in  the  broad-spectrum  mode,  poses 
special  problems  because  of  the  sharp  skirts  of  the  transducer  responce  and 
spurious  resonances  due  to  the  transducer  casement  and  mount.  Evading 
these  problems  and  defining  the  transducer  effects  on  the  data  are  significant 
achievements  of  the  system  effort.  A  thorough  investigation  of  the  transmission 
and  reception  methods,  data  processing,  and  calibration  is  necessary  to 
facilitate  interpretation  of  the  experimental  results  and  to  establish  their 
validity. 

The  measurements  program  deals  primarily  with  randomly  rough  sur¬ 
faces  having  known  statistics.  Measurements  of  backscatter  power  are  re¬ 
corded  for  a  variety  of  target  surfaces  at  various  incidence  angles  across 
a  frequency  range  of  30:1  (X  =  15  mm  to  0.5  mm) .  The  targets  are  designed 
to  have  correlation  functions  which  Fourier  transform  to  functions  almost 
entirely  contained  within  the  illuminating  frequency  spectrum.  The  text 
will  classify  targets  by  their  roughness  density  spectrum  which  is  the  Fourier 
transform  of  the  correlation  function. 

Several  approaches  to  describing  scattering  of  electromagnetic  and 
acoustic  waves  have  been  developed  in  recent  years.*  Each  such  approach  has 


*  Fung,  A.  K.,  CRES  Report  48-3,  Scattering  Theories  and 
Radar  Return,  1965 


a  frequency  dependence  whether  intentionally  or  incidently  introduced.  To 
accomplish  the  comparison  of  experimental  results  with  theoretical  predic¬ 
tions,  a  detailed  review  of  the  pertinent  theories  is  necessary.  In  each 
case  careful  attention  must  be  given  to  the  necessary  surface  condition  imposed 
by  the  development. 

The  strength  of  this  overall  study  rests  in  the  use  of  broad-spectrum 
Illumination  to  determine  the  frequency  dependence  of  backscatter.  The 
technique  provides  information  for  a  broad,  continuous  range  of  frequencies, 
thus  eliminating  extrapolation  errors  and  unknowns.  It  is  an  extremely 
flexible  experimental  approach.  In  essence  it  is  a  new  concept  in  radar 
systems  and  the  overall  study  includes  an  evaluation  of  this  concept. 


SYSTEM  INSTRUMENTATION 

The  system  used  for  this  measurements  program  was  designed  and 
built  by  the  Center  for  Research,  Inc.,  Engineering  Science  Division.  Fund¬ 
ing  was  shared  by  the  Department  of  Defense  under  contract  DAAK02-68-C-0089 , 
the  National  Aeronautical  and  Space  Administration  under  contract  NAS  9- 
7175,  and  the  Center  for  Research,  Inc.  The  need  for  flexibility  heavily 
influenced  the  selection  of  subassemblies  throughout  the  instrument,  and 
in  many  cases  the  subassembly  used  has  many  times  the  actual  capability 
required . 

A.  System  Function 

To  produce  data  necessary  for  precisely  describing  the  frequency  depen¬ 
dence  of  backscatter,  the  experiment  system  ideally  must  transmit  a  uniform 
amplitude,  variable-width,  broad-spectrum  acoustic  wave  in  a  frequency 
independent  beam.  Such  a  system  cannot  be  built;  however, an  approximation 
to  this  system  is  possible. 

The  system  as  built  operates  in  what  is  termed  a  "slow-sweep"  mode. 
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In  this  mode  the  transmitted  frequency  is  swept  from  minimum  to  maximum 

in  a  fixed  time.  The  transmitter  is  gated  on  by  a  narrow  pulse  at  a  repeti- 
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tion  rate  from  10  to  10  faster  than  the  sweep  rate,  hence  the  transmitter 
output  consists  of  a  pulse  train  in  which  each  pulse  contains  a  different 
carrier  frequency.  This  technique  is  illustrated  in  figure  A.l .  The  approach 
is  actually  a  multiple -discrete  frequency  method  where  several  hundred 
discrete  frequencies  are  used  and  the  spacing  is  such  that  the  output 
spectrum  of  one  overlaps  the  spectra  of  neighboring  pulses. 

The  envelope  of  the  average  amplitudes  of  the  returns  associated 
with  pulses  t^  through  t^  is  the  average  return  power  spectrum  for  the  particular 
target  being  tested.  The  handling  of  the  data  necessary  to  obtain  this  spec¬ 
trum  is  complex  and  is  described  in  detail  later. 


Figure  A.l.  Slow-sweep  mode  sampling  technique 
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The  principle  of  operation  of  the  broad -spectrum,  slow-sweep  mode 
is  simple,  however  the  implementation  is  difficult.  The  primary  problem 
is  the  piezoelectric  transducers  used  as  the  transmitting  and  receiving  antennas. 
Their  beamwidth  and  gain  are  highly  frequency  sensitive,  and  careful  atten¬ 
tion  is  needed  to  insure  the  validity  of  the  experimental  results . 

B.  System  Description 

The  block  diagram  of  the  system  is  shown  in  figure  1 ,  It  consists 
of  three  main  sections:  1)  transmitting,  2)  receiving,  and  3)  data  processing. 

1.  Transmitting  Section 

This  section  consists  of  a  function  generator,  swept 
frequency  oscillator,  pulse  generator,  gate,  gated  power  ampli¬ 
fier,  transformer,  and  a  transmitting  transducer. 

The  function  generator  provides  two  outputs .  One  output 
Is  a  triangle  wa  -..orm  on  a  dc  voltage  level.  This  output  drives 
the  swept  oscillator.  The  other  output  is  a  square  wave  and  is 
used  to  gate  the  pulse  generator  "on”  during  the  interval  of 
positive  slope  of  the  triangle  waveform,  and  "off”  during  fhe 
negative  slope  interval. 

The  swept-frequency  oscillator  provides  a  leveled  out¬ 
put  across  the  frequency  range  programmed  by  the  function  genera¬ 
tor.  The  output  is  fed  through  a  gate  and  into  a  gated  power  ampli¬ 
fier.  Each  gate  is  switched  "on”  (closed)  by  the  output  pulse 
of  the  pulse  generator.  The  purpose  of  the  double  gate  is  to  in¬ 
crease  the  "on"  -  "oft"  ratio  and  thus  minimize  the  noise  sub¬ 
sequently  experienced  at  the  receiver  input. 

The  power  amplifier  has  a  gain  of  approximately  30  over 
the  frequency  range  0.1  MHz  to  3.0  MHz.  The  push-pull  output 
stage  in  the  tester  drives  25  feet  of  23  ohm  cable  (two  50  cables 
in  parallel)  which  connects  to  a  3:1  step-up  transformer  and  a 
200  Q  terminating  resistor.  In  addition  a  small  inductor  is  in 
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series  with  the  transducer  to  form  a  low  pass  filter  arrangement. 

2.  Receiving  Section 

This  section  consists  of  a  receiving  transducer,  cable 
driving  amplifier,  low-noise  preamplifier,  and  final  amplifier. 

The  cable  driving  amplifier  is  a  unity  gain  impedance 
transformer  having  a  92  R  output  impedance  matched  to  the  cable . 

The  low -noise  preamplifier  provides  80  db  gain  over  the 
frequency  range  of  0.1  MHz  to  3.0  MHz  with  a  noise  voltage  of 
less  than  30ii  volts .  The  unit  has  a  92  output  impedance  and  is 
divided  into  two  stages  with  variable  attenuation  between  stages. 

The  final  amplification  is  obtained  in  a  Tektronix  oscil¬ 
loscope  preamplifier  (type  CA) . 

3.  Data  Processing  Section 

The  fundamental  unit  in  the  data  processing  section  is  a 
256  channel  pulse  height  analyzer.  This  unit  counts  the  number 
of  pulses  of  equal  amplitude  and  stores  the  counts  in  a  storage 
element  corresponding  to  the  pulse  amplitude.  Each  of  the  256 
storage  channels  represents  a  particular  pulse  amplitude  range, 
and  each  count  in  the  channel  indicates  the  number  of  pulses 
occurring  within  that  amplitude  range . 

In  this  system  the  256  channels  represent  frequency  ranges 
rather  than  amplitude  ranges.  The  data  processing  section  con¬ 
verts  the  return  amplitude  at  a  particular  frequency  into  a  number 
of  counts  at  a  set  amplitude  corresponding  to  the  return  signal 
frequency.  This  is  accomplished  by  converting  the  analog  return  to 
digital  and  causing  the  digital  count  pulse  to  vary  in  amplitude 
in  direct  proportion  to  the  swept  frequency. 

Following  figure  2  .,  the  data  processing  operation  is 

as  follows.  The  return  signal  is  obtained  from  the  VERT  SIG  OUT 
of  the  Tektronix  545B  oscilloscope.  The  signal  is  sampled  in  a 
high  speed  solid-state  switch  controlled  by  the  Tektronix  545 B 
+A  GATE.  The  sample  is  stretched  and  compared  to  a  ramp  function 
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in  a  differential  amplifier.  The  differential  amplifier  gives  a 
-3v  output  for  an  interval  equal  to  the  time  during  which  the 
stretched  pulae  is  greater  than  the  ramp.  This  interval  serves 
as  a  gate  on  a  5  KHz  clock  (5  KHz  is  the  maximum  rate  acceptable 
by  the  pulse  height  analyzer) .  The  clock  output  is  compared 
with  the  triangle  waveform  driving  the  swept  frequency  oscil¬ 
lator.  The  comparison  is  made  by  an  operational  amplifier  that 
changes  gain  proportional  to  the  triangle  input.  The  result  is 
that  the  pulse  height  analyzer  receives  a  pulse  train  of  ampli¬ 
tude  proportional  to  frequency  and  duration  proportional  to  return 
amplitude . 

The  results  of  several  sweep  intervals  are  stored  and  the 
average  spectrum  is  displayed  on  a  CRT.  The  CRT  is  photographed 
and  subsequently  analyzed. 

C.  System  Calibration 

The  calibration  program  is  concerned  with  four  system  sections. 

These  are  1)  power  amplifier,  including  output  transformer  section,  2)  low- 
noise  amplifier,  3)  data  processor,  including  pulse  height  analyzer,  and  4) 
transducers.  Of  these,  the  transducers  occupy  the  prime  attention. 

1.  Low-noise  Preamplifier 

The  low-noise  preamplifier  consists  of  two  40  db  gain 
Integrated  circuit  amplifiers  (RCA  3022) .  Each  unit  has  an  upper 
3  db  point  at  7.5  MHz  with  a  noise  figure  of  appioximately  5  db. 
The  two  stages  are  separated  by  a  0  -  101  db  variable  attenuator. 
The  gain  of  the  unit  is  55  db  +  1  db  across  the  frequency  range 
0  J.  -  3.0  MHz. 

2 .  Power  Amplifier 

The  power  amplifier  consists  of  two  sections.  The  first  is 
a  basic  push-pull  amplifier  with  an  optional  gating  feature.  The 
gain  of  this  stage  is  10  +.  .1  over  the  frequency  range  0.1  MHz 
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to  3.0  MHz.  The  second  stage  consists  of  a  step-up  transformer, 
termination,  and  low -pass  filter.  The  transformer  ratio  is 
2.65  +0.1  over  the  frequency  range.  The  assembly  is  separated 
from  the  amplifier  by  25  feet  of  25  ohm  coxial  cable.  The  sepa¬ 
ration  is  advantageous  to  the  operation  since  it  isolates  the  cable 
capacity  from  the  transducer  capacity.  The  transducer  impedance 
is  highly  capacitive  and  the  variation  of  impedance  from  unit- 
to-unit  is  great.  This  requires  special  attention  in  design  since 
the  maximum  power  deliverable  is  a  function  of  the  particular 
transducer  cloud,  .once  calibration  is  required  individually  accord¬ 
ing  to  the  transducer  in  use . 

3.  Data  Processor 

The  data  processor  calibration  consists  of  determining  the 
input  signal  voltage  required  to  produce  each  addition  count. 

The  voltage  versus  count  curve  is  linear  and  has  a  slope  of  0.066  • 
volts/count  with  an  intercept  of  0.15  volts  at  8  counts.  The 
reference  count  is  six,  i.e.  with  no  signal  in,  a  stx  count  is 
present. 

4.  Transducers 

Measurements  in  the  acoustic  system  are  made  by  means 
of  piezoelectric  transducers .  This  program  used  the  simple  flat 
disk,  barium  titanate  transducer  for  all  data.  These  transducers 
were  designed  for  use  at  a  single,  electrical  and  mechanical  re¬ 
sonance  frequency.  Attempting  to  employ  them  in  a  broad  spec¬ 
trum  system  created  special  problems. 

It  would  be  advantageous  to  use  only  one  pair  of  transducers 
to  record  measurements  over  the  entire  frequency  range  of  interest. 
Unfortunately,  no  transducer  was  located  that  would  meet  the 
requirements.  The  basic  difficulty  is  that  the  transducer  gain 
pattern  is  frequency  dependent  to  the  extent  that  reasonable  size 
beam  patterns  at  3.0  MHz  become  too  large  for  the  limited  size 
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water  tank  at  0.1  MHz.  Additional  problems  arise  with  regard 

to  resonance  and  harmonics. 

In  most  documented  acoustic  wave  experiments,  trans¬ 
ducer  gain  is  measured  by  transducer  pairs  and  seldom  is  the  gain 
of  an  individual  transducer  obtained.  In  this  study  an  attempt  was 
made  to  determine  the  absolute  gain  of  a  single  transducer.  To 
accomplish  this  two  type  measurements  were  necessary, 
a)  Single  Transducer  Measurement 

The  test  setup  for  this  measurement  is  shown 
in  figure  3.  A  single  transducer  is  used  as  both 
the  transmitting  and  receiving  antennas.  The  trans¬ 
ducer  is  oriented  so  that  the  acoustic  wave  propagation 
is  normal  to  a  0.25  inch  thick  plate  glass  target. 

The  operation  for  single  transducer  measure¬ 
ments  Is  as  follows.  The  sweep  frequency  oscillator 
is  set  to  a  fixed  frequency.  The  output  of  the  oscilla¬ 


tor  is  gated  twice  to  provide  a  high  on-off  ratio,  and 

transmitted.  On  transmission,  the  switch  is  in  position 

A.  After  the  transmitted  pulse  interval,  the  switch  is 

thrown  to  position  B,  and  the  return  signal  enters  the 

preamplifier.  The  measurement  is  the  ration  of 

V  The  transducer  gain  is  calculated  from  the 
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equation. 


G  = 


47r  (2h) 
\K 


C.l. 


where  2h  is  the  distance  to  the  target,  \  is  the  wave¬ 
length,  and  K  is  the  reflection  coefficient,  which  is  un¬ 
known.  Equation  C.l .  is  developed  from  the  equation 
for  one-way  propagation  by  assuming  the  receiving 
transducer  is  at  the  image  point  of  the  transmitting 
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transducer  and  that  K  =  1 .  The  correct  value  of  K  is 
determined  later. 

Plate  glass  is  used  as  the  target  to  insure  a 
smooth  surface  to  all  wavelengths  in  the  range  of 
interest.  The  calculated  reflection  coefficient  for 
glass  in  water  is  0.83. 

As  will  be  evident  in  the  data,  using  the  plate 
glass  as  a  target  created  some  unexpected  difficulties. 
The  glass  target  is  approximately  0.25  inches  thick. 

This  width  is  equal  to  one  wavelength  at  about  1  MHz 
(velocity  of  propagation  of  sound  in  glass  is  6  km/sec) . 
Assuming  the  glass  acts  as  a  low-loss  transmission  line, 
the  terminating  impedance  on  the  backside  of  the  plate 
is  reflected  to  the  front  side  and  thus  presents  a  matched 
impedencc  termination  to  the  incident  energy.  The 
matched  condition  should  appear  at  about  0.5  MHz, 

1.0  MHz,  1.5  MHz,  2.0  MHz,  2.5  MHz,  and  3.0  MHz 
based  on  the  thickness  of  the  glass.  The  results  confirm 
the  condition  at  these  frequencies.  Fortunately  the  nulls 
are  very  sharp  for  frequencies  above  1  MHz.  The  two 
nulls  at  1  MHz  and  0.5  MHz  are  generally  broad,  and 
hence  the  data  are  suspect  in  this  range, 
b)  Transducer  Pair  Measurements 

The  test  set  up  for  this  measurement  is  shown 
in  figure  4.  Two  transducers  are  aligned  directly 
facing  each  other,  1.0  meter  apart.  The  measurement 
recorded  is  the  ratio  of  voltage  at  the  receiving  trans¬ 
ducer  to  the  voltage  at  the  transmitting  transducer, 

'^OUT^IN  • 

In  this  arrangement  signals  are  recordable  for 
the  entire  frequency  range  of  0.1  MHz  to  3.0  MHz  for 
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nearly  all  transducers  .  The  gain,  either  receiving 
or  transmitting,  is  calculated  using  the  equation  for 
one-way  propagation. 


C.2. 


where:  =  gain  of  transmitting  (or  receiving)  trans¬ 

ducer 

Gj.  =  gain  of  receiving  (or  transmitting)  trans¬ 
ducer  measured  by  single  transducer  method 

R  “  spacing  between  transducers 

X  =  wavelength 

Although  it  was*  mentioned  that  the  entire  fre¬ 
quency  range  can  be  covered  in  the  pair  measurements, 
can  only  be  calculated  over  small  segments  .  This 
is  due  to  limited  number  of  values  of  G^  (equation  C.2 .) 
available  for  calculation  from  the  single  transducer 
measurements . 

The  data  plots  are  very  similar  in  shape  to  those 
obtained  in  the  single  transducer  measurement,  especially 
above  1  MHz.  However,  there  is  a  difference  in  magni¬ 
tude  of  14  -  15  db.  This  difference  is  due  to  the  fact 
that  the  reflection  coefficient,  K,  was  assumed  to  be 
unity  in  equation  C.l .  The  actual  value  of  K  can  be 
determined  using  both  the  single  and  pairs  data. 


where  h  =  height  above  target  (single) 
■R  =  spacing  (pairs) 
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^IN  '^OUT  =  input  and  output  of  receiving 

transducer  (single) 

Vjj^*  and  =  input  and  output  of  transmitting 

transducer  (single) 

ViN  and  =  input  and  output  of  pairs  measure- 

ment 

The  value  of  K  for  the  3.0  MHz  12®  and  6®  pairs  is 
0.21 .  This  compares  with  a  calculated  value  for  K 
of  0.83.  The  difference  is  most  likely  due  to  a  com¬ 
bination  of  effects  such  as  1)  difference  in  the  power 
delivered  to  the  transducer  mount  as  compared  to  that 
actually  delivered  to  the  water,  2)  difference  in  actual 
properties  of  the  glass  target  and  ideal  properties  of 
glass;  density  equal  2400  kb/cu  3  and  velocity  of  pro¬ 
pagation  equal  6  km/sec,  and  3)  transmission  line 
impedance  transformation  effects . 

Using  a  K  value  of  0.21  the  receiving  gain  is 
computed  for  each  transducer  used  in  the  experiment 
both  as  individual  units  and  in  pairs . 

The  disk  or  piston  transducers  used  in  this  work 
have  an  active  surface  area  measuring  several  wave¬ 
lengths  across  and  produce  a  directive  antenna  pattern. 
The  far  field  pattern  for  a  flat  disk  transducer  operating 
in  an  infinite  baffle  is  given  by 


F  (0)  = 


2  Jj  (k  a  sin  ©) 
k  a  sin  0 


C.  4. 


where  =■■  Bessel  function  of  first  kind  and  first  order 

k  =  wave  number 

a  =  radius  of  flat  disk 

0  =  angle  from  surface  normal 
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D.  Summary 

The  broad  spectrum  acoustic  system  provides  accurate,  high  sensi¬ 
tivity  measurements  over  the  frequency  range  0.1  MHz  to  3.0  MHz.  The 
calibration  tests  confirm  that  no  single  transducer  will  operate  satisfactorily 
over  the  full  range  and  hence  data  must  be  recorded  in  segments  as  dictated 
by  the  transducers  used. 

Acknowledgement;  This  work  was  supported  by  Project  Themis  (USAETL 
Contract  DAAK02-68-C-0089 ,  ARPA  order  No.  1079)  and  NASA  Contract 
NAS9-7175. 


Figure  3.  SINGLE  TTRANSDUCER  TEST 
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ANTENNA  AND  OUTPUT  POWER  CONSIDERATIONS 
FOR  THE  POLYPANCHROMATIC  RADAR  SYSTEM 

* 

by 

W.  P.  Waite 

1.  INTRODUCTION 

The  polypanchromatic  radar  is  a  device  lo/  performing  broad- 
spectrum  electromagnetic  scattering  and  imaging  measurements.  The 
system  is  designed  to  operate  from  a  truck-mounted  boom  and  is 
completely  self-contained  and  easily  transportable  to  s'^rious  target 
areas  of  interest.*  The  system  design  is  necessarily  specialized; 
however/  the  individual  components  of  the  system  are  standard  commer¬ 
cially  available  laboratory  and  test  equipment. 

This  memorandum  examines  the  requirements  and  performance  of 
the  system  both  from  a  general  standpoint  and  in  view  of  the  components 
available  from  industry  and  within  the  budget  limitations  of  the  THEMIS 
contract.  A  specific  antenna  and  transmitter  are  recommended  for  pur¬ 
chase  and  use  as  a  result  of  this  study. 

2.  GENERAL 

The  basic  design  goals  of  the  polypanchromatic  radar  systam  are 

1 .  Operation  over  a  bandwidth  of  at  least  one  octave . 

2.  Measurement  of  backscattering  cross  section  as  a  function  of 

a)  Frequency 

b)  Bandwidth 

c)  Incidence  Angle  (from  0°  to  as  near  grazing  as  possible) 

3.  Production  of  panchromatic  imagery  as  a  function  of 

a)  Frequency 

b)  Bandwidth 

*  The  measurement  program  for  which  this  equipment  is  to  be  used  will  be 
similar  to  the  discrete  frequency  measurements  described  in  the  following 
r©f  ©r©nc0-* 

Cosgriff,  R.  L.,  W.  H.  Peake,  and  R.  C.  Taylor.  "Terrain  Scattering 
Properties  for  Sensor  System  Design" ,  Terrain  Handbook  II,  Engineer¬ 
ing  Experiment  Station,  Ohio  State  University,  Antenna  Lab,  Columbus, 

Ohio,  1959. 
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4.  Maximum  usaga  of  commercially  available  components. 

Figure  1  shows  a  block  diagram  of  the  proposed  system.  This 
system  is  the  same  as  that  shown  in  the  original  THEMIS  proposal 
with  the  exception  that  the  operating  frequency  range  has  been  changed 

from  8,0-12.0  GHz  to  4. 0-8.0  GHz. 

In  the  following  sections  the  desired  resolution  and  range  are 
examined  for  their  effect  upon  the  antenna  configuration  and  the  required 
output  power.  These  effects  are  examined  in  view  of  the  commercial 
equipment  available  and  specific  model  numbers  are  recommended  for 
purchase  as  a  result  of  this  study. 

3.  ANTENNA  CONSIDERATIONS 

The  truck  upon  which  the  system  is  to  be  mounted  will  be 
furnished  as  GFE  and  B.  B.  Scheps  (ETL  Project  Engineer  for  this 
contract)  has  indicated  the  boom  height  above  ground  of  the  available 
vehicle  to  be  42  feet.  This  height  refers  to  the  vertical  distance  from 
the  ground  to  the  bottom  of  the  bucket  when  the  boom  is  fully  extended 
and  in  the  vertical  position.  For  performance  of  the  desired  backsca tier¬ 
ing  measurements  the  return  at  vertical  incidence  must  be  measured. 

In  this  mode  of  operation  the  minimum  range  of  the  system  is  determined 

by  the  antenna  height  above  grour  ’. 

Figure  2  illustrates  the  deployment  of  the  boom  when  performing 
backscattering  measurements.  The  angle  between  the  boom  and  vertical 
should  be  at  least  30  degrees.  This  is  required  not  only  to  position 
the  antenna  over  the  target  area,  but  to  prevent  the  tnick  itself  from 
significantly  contributing  to  the  return  through  the  antenna  sidelobes . 
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Figure  1  —  Polypanchromatic  Radar  System 
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Figure  2  —  Boom  Deployment  for  Backscattering  Measurements 


If  we  assume  the  antenna  to  be  mounted  at  about  the  same  level 
as  the  bottom  of  the  bucket,  the  antenna  height  above  ground  or  the 
minimum  required  range  is  given  approximately  by 

R  ,  =  R  cos  ©“  42  (0.866)  =  36.4  feet 
min 

=  11.1  meters 


This  minimum  range  corresponds  to  a  transmitted  pulse  length  of 


T  = 


2x11.1 
3  x  10® 


=  74  nsec 


where  t  =  transmitted  pulse  length 

R  ,  =  minimum  range 

min 

c  =  velocity  of  light 

To  allow  sufficient  time  for  receiver  recovery  after  the  main  bang 
the  maximum  pulse  length  should  not  exceed  50  nsec  and  operation  with  even 
shorter  pulse  lengths  may  well  be  desirable. 

The  resolution  requirements  imposed  upon  the  antenna  will  be 
distinctly  different,dependent  upon  the  mode  of  operation  of  the  system 
(backscattering  measurement  or  imaging) .  The  azimuth  resolution  of  the 
system  must  be  accomplished  by  beamwidth  limitation  of  the  antenna  in 
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either  mode  of  operation.  The  range  resolution  may,  however,  be  pro¬ 
vided  by  either  beamwidth  limitation  of  the  antenna  or  by  the  pulse  length 
limitation  of  the  ground  illumination.  Figure  3  shows  the  various  types  of 
azimuth  and  range  resolution  geometry. 


a)  Azimuth  Resolution  Geometry 


b)  Range  Resolution  Geometry  for  Antenna  Beam 
Angle  Limitation 


I 

i 


I 


t 
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c)  Range  Resolution  Geometry  for  Pulse  Length  Limitation 
Figure  3a,  3b,  3c  —  Resolution  Geometry 
The  azimuth  resolution  is  given  approximately  by 


Pa  RP  ==  R  ^  ^  cos  e 

where  p  =  azimuth  resolution 
a 

R  =  range  of  resolution  cell 
p  =  antenna  azimuth  beam  angle 
\  =  wavelength 
^  =  length  of  antenna  aperture 
h  =  height  of  antenna 
0  =  incidence  angle 


(1) 


The  resolution  cell  on  the  ground  when  obtained  by  beam  limita¬ 
tion  in  both  azimuth  and  range  will  actually  approximate  an  ellipse  (ex¬ 
cept  at  vertical  where  it  is  circular)  when  the  antenna  aperture  is  circular^ 
such  as  with  a  parabolic  reflector.  The  above  expression  for  azimuth 


I 


-t 


resolution  is  then  the  minor  axis  of  the  elliptical  illuminated  area.  The 
major  axis  dimension  is  determined  from  Figure  3  (b)  and  is  given  by 

=  h  (tan  <^2  “ 

where  ^  ^1“®”^ 

and  p^  =  range  resolution  (length  of  major  axis) 

h  =  height  of  antenna 

0  =  pointing  angle 

a  =  antenna  elevation  beam  angle 

The  range  resolution  when  determined  by  pulse  length  illumination 
cf  the  ground  is  shown  in  Figure  3  (c)  and  is  given  approximately  by 


CT 

2  sin  e 


e  >  0^^  (3) 


where  p^  =  range  resolution 

T  =  transmitted  pulse  length 
c  =  velocity  of  light 
0  =  incidence  angle 

The  above  expression  for  range  resolution  is  accurate  only  for  0  grv  jter 
than  approximately  10°  and  for  ^  «  R.  At  vertical  incidence  the  resolu¬ 
tion  cell  provided  by  pulse  length  limitation  is  circular^ with  a  diameter 
given  approximately  by 


D  =  2  VhcT* 
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where  D 
h 
c 

T 


=  diameter  of  illuminated  area  (at  vertical) 
=  height  of  antenna 
=  velocity  of  light 
=  transmitted  pulse  length 


The  antenna  resolution  requirements  for  each  mode  of  operation 
(back scattering  or  imaging)  will  differ  and  will  consequently  first  be 
examined  separately. 


3  •  1  Backscatterinq  Mode  of  Operation 

In  the  backscattering  mode,  measurements  must  be  made  at 
and  near  vertical  incidence.  Examination  of  Equation  (4)  shows  that 
pulse  length  limited  operation  is  not  practical  in  this  region  of  operation, 
for  if  we  assume 


h  =  3.64  feet 


then 


T  =  10  nsec  (the  minimum  practical 
pulse  length) 


D  =  2  >/hCT  »  38  feet 

and  the  diameter  of  the  illuminated  area  is  far  too  large . 

The  operation  of  the  system  must  be  with  beamwidth  limitation 
at  angles  near  vertical  incidence .  The  diameter  of  the  circular  illuminated 
area  at  vertical  incidence  is  then  given  by  Equation  (1) . 

It  may  also  be  seen  from  Equation  (1)  that  the  beamwidth  is  a 
function  of  the  operating  frequency.  Correction  for  the  change  in  illuminat¬ 
ed  area  must  be  made  during  the  sweep  in  frequency. 

Table  1  lists  maximum  and  minimum  resolution  cell  diameters  at 
vertical  incidence  for  a  variety  of  available  parabolic  reflector  sizes. 


•i 
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TABLE  1  —  Variation  of  Illuminated  Area  at  Vertical  Incidence 
h=  36.4  0=  0° 


Beam 

Angle 

Diameter  of  Illuminated  Area 

Reflector  Size 
(feet) 

Min 

(8.0  GHz) 
(degrees) 

Max 

(4.0  GHz) 
(degrees) 

Min 

(8.0  GHz) 
(feet) 

Max 

(4.0  GHz) 
(feet) 

1 

7.0 

14.1 

4.5 

9.0 

2 

3.5 

7.0 

2.2 

4.5 

3 

2.3 

4.7 

1.5 

3.0 

4 

1.8 

3.5 

1.1 

2.2 

6 

1.2 

2.3 

0.7 

1.5 

8 

0.9 

1.8 

0.5 

1.1 

Table  1  illustrates  the  advantage  of  a  large  reflector,  which  is 
desirable  to  constrain  the  angular  variation  within  the  beam.  However, 
at  this  relatively  short  range  a  large  reflector  reduces  the  total  illumin¬ 
ated  area  to  such  a  small  size  that  only  targets  or  areas  of  relatively 
small  scale  structure  and  size  may  be  adequately  illuminated.  The 
best  compromise  between  these  conflicting  requirements  seems  to  be  a 
reflector  size  of  2  or  3  feet  for  use  near  vertical  incidence. 

Table  2  shows  the  variation  of  azimuth  and  range  resolution  as 
a  function  of  incidence  angle.  Range  resolution  through  pulse  length 
limitation  is  included  for  angles  away  from  the  vertical.  The  figures 
of  Table  2  show  beam  limited  operation  to  work  quite  satisfactorily  out 
to  an  Incidence  angle  of  approximately  60°.  After  this  point,  however, 
range  resolution  must  be  provided  by  pulse  length  limitation. 

Another  interesting  aspect  illustrated  by  Table  2  is  the  variation 
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of  slant  range  with  incidence  angle.  Angles  up  to  85  degrees  may  be 
measured  at  a  range  of  less  than  500  feet^even  with  the  boom  at  maximum 
elevation.  If  angles  even  nearer  grazing  are  desired,  it  would  be  a  simple 
matter  to  decrease  the  boom  height  and  measure  out  to  approximately  89° 
with  little  increase  in  range. 

For  instance,  with  an  assumed  boom  height  of  8  feet^the  slant 
range  required  for  an  incidence  angle  of  89  degrees  is  458  feet.  Thus, 
the  maximum  range  required  of  the  system  in  the  scattering  mode  is  well 
under  1000  feet  and  the  maximum  range  capability  will  probably  be  set  by 
the  Imaging  requirement. 

3.2  Imaging  Mode  of  Operation 

In  the  imaging  mode  it  is  desirable  to  maximize  the  total  range 
covered  in  elevation.  In  this  mode  the  azimuth  resolution  must  again 
be  provided  by  the  beam  limiting  of  the  antenna;  however,  the  range 
resolution  must  be  provided  by  pulse-length  limitation. 

The  difficulty  here  is  one  of  procurement.  The  only  reflector 
shape  compatible  with  the  necessary  broadband  feeds  and  available  as 
a  standard  item  is  a  paraboloid.  The  azimuth  and  elevation  beam  angles 
are  thus  approximately  equal. 

This  offers  another  contradiction  in  requirements.  The  reflector 
size  should  be  large  to  provide  resolution  at  the  maximum  range.  On 
the  other  hand,  the  reflector  size  should  be  small,  with  the  resultant 
large  beam  angle  in  elevation,  to  provide  maximum  coverage  of  the  ground 
from  near  to  far  range . 

The  requirements  of  the  system,  then,  would  best  be  met  by  a 
rectangular  aperture  giving  fan  beam  illumination  of  the  area.  This  may 
possibly  be  obtained  by  locating  a  manufacturer  who  has  produced  a  rec¬ 
tangular  reflector  of  the  required  size  for  some  previous  project  and  would 
supply  another  of  the  same  design  for  a. reasonable  (within  the  budget) 
price.  Alternatively,  it  may  be  possible  to  procure  a  standard  circular 
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dish  with  the  desired  maximum  dimension  and  then  deliberately  spoil 
the  elevation  beam.  This  might  be  ctone  by  adding  absorbing  material 
as  shown  in  Figure  4, 


Figure  4  -  Possible  Means  of  Elevation  Beam  Spoiling 

Virtually  every  major  manufacturer  of  broadband  antennas  has 
already  been  contacted  and  none  of  them  has  available  as  a  standard 
or  semi-standard  item  the  rectangular  shape  desired.  The  modification 
of  a  standard  item  such  as  shown  in  Figure  4  promises  to  be  both  time 
consuming  and  expensive.  It  is  recommended  that  the  imaging  system 
be  designed  for  operation  with  the  standard  parabolic  dish  sizes  while 
continuing  the  search  for  a  more  desirable  aperture  that  might  be  added 
at  a  later  date . 

Imaging  operation  with  a  parabolic  dish  will  restrict  the  elevation 
beam  angle  such  that  the  angle  of  incidence  must  be  very  near  grazing 
if  a  reasonable  range  is  to  be  illuminated.  This  poses  the  problem  of 
operating  in  an  undesirable  range  of  incidence  angles.  This  might  be 
at  least  partially  alleviated  by  selecting  operating  locations  on  steep 
hills  or  bluffs  and  thus  increasing  the  effective  antenna  height  above 
the  ground  being  imaged. 

Table  3  shows  the  variation  of  slant  range  with  incidence  angle 
for  a  variety  of  effective  antenna  heights . 
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Table  3  -  Slant  Range  as  a  Function  of  Incidence  Angle 


Slant  Range 


Incidence  Angle 
(degrees) 

h  =  40  ft. 
(feet) 

h=  100  ft. 
(feet) 

h  =  200  ft. 
(feet) 

h=  300  ft. 
(feet) 

90 

89 

2,290 

5,730 

11,500 

17,200 

88 

1,150 

2,860 

5,720 

8,580 

87 

765 

1,910 

3,820 

5,730 

86 

572 

1,430 

2,860 

4,300 

85 

459 

1,150 

2,290 

3,440 

84 

382 

956 

1,910 

2,870 

83 

328 

821 

1,640 

2,460 

82 

287 

718 

1,440 

2,150 

81 

256 

639 

1,280 

1,920 

80 

230 

576 

1,150 

1,730 

78 

192 

481 

960 

1,440 

76 

165 

413 

825 

1,240 

74 

145 

363 

725 

1,090 

72 

129 

324 

647 

970 

70 

117 

292 

585 

876 

68 

107 

267 

533 

800 

66 

98 

246 

491 

736 

64 

91 

228 

456 

685 

62 

85 

213 

426 

640 

60 

80 

200 

400 

600 

55 

70 

174 

349 

523 

50 

62 

156 

311 

466 

45 

57 

141 

282 

424 

f3( 

2(1 
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Some  estimate  of  the  maximum  range  of  the  system  must  now  be 
made.  We  will  assume  that  the  maximum  usable  resolution  is  100  feet 
(in  either  azimuth  or  range) .  Table  4  lists  the  standard  reflector  sizes 
and  the  maximum  range  corresponding  to  an  azimuth  resolution  of  100  feet. 
The  frequency  used  is  the  lowest  operating  frequency  (4.0  GHz). 

Table  4  -  Maximum  Range  for  100  Foot  Resolution 


Reflector  Size 

Beam  Angle 
(4.0  GHz) 

Maximum  Range 

8 

1.8 

3,180 

6 

2.3 

2,490 

4 

3.5 

1,640 

3 

4.7 

1,220 

2 

7.0 

820 

1 

14.1 

410 

From  a  comparison  of  Tables  3  and  4  several  interesting  conclusions 
may  be  drawn.  First,  even  with  an  effective  antenna  height  of  300  feet/ 
a  range  of  over  approximately  3,500  feet  will  place  operation  in  the  region 
of  85  to  90  degrees  incidence  angle.  This  is  undesirable  for  imaging  of 
real  targets, as  airborne  imagery  is  sparse  in  this  region.  The  images  ob¬ 
tained  in  this  program  should,  as  nearly  as  possible,  simulate  those  which 
might  be  obtained  by  an  aircraft  system.  This  would  require  operation 
below  an  incidence  angle  of  approximately  80  degrees.  For  operation  in 
this  range  of  incidence  cngles,  the  range  is  even  less  for  any  reasonable 
assumption  of  effective  antenna  height. 

Second,  if  incidence  angles  belgw  80  degrees  are  desirable,  the 
use  of  a  smaller  antenna  may  be  advisable  to  maximize  the  angular  elevation 
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coverage.  It  is  possible  to  use  even  a  2  or  3  foot  dish, as  the  range 
corresponding  to  these  incidence  angles  is  sufficiently  short  that  the 
azimuth  resolution  is  still  reasonable. 

From  the  above  considerations  it  is  apparent  that  long  range  cap¬ 
ability  is  not  required  if  the  purpose  of  the  imaging  mode  is  to  investigate 
effects  in  normal  airborne  imaging  systems .  Some  argument  for  increased 
range  capability  might  be  offered,  however,  if  the  system  is  to  be  used 
to  investigate  the  effects  of  broad  bandwidth  on  ground-based  battlefield 
surveillance  radar.  The  frequency  range  of  the  Phase  I  system  (4  to  8  GHz) 
is  probably  higher  than  would  be  desired  for  this  application.  Some  in¬ 
vestigation  of  the  foliage  penetration  and  resonance  effects  in  metallic 
targets  would  be  possible  in  this  area,even  at  short  ranges.  The  de¬ 
creased  range  would  merely  prevent  complete  simulation  of  an  operational 
system,  but  the  parameter  of  interest  in  this  mode  could  still  be  investi¬ 
gated  . 

There  is  another  problem  presented  by  the  short  range  obtainable 
at  lower  incidence  angles, even  for  the  simulation  of  airborne  radar  imagery. 
This  is  that  the  spatial  distribution  of  composite  targets  (such  as  fields, 
grassland,  woods)  is  greater  than  the  system  range  at  normal  incidence 
angles .  Thus ,  the  range  of  an  image  would  only  span  one  or  two  such 
composite  targets.  This  objection  may  be  partially  overcome  by  increas¬ 
ing  the  effective  antenna  height  through  careful  site  selection.  Another 
possibility,  for  agricultural  investigations,  would  be  to  image  a  small 
agricultural  experiment  station  where  the  individual  fields  are  quite  small 
(50  to  100  feet  maximum  extent)  and  the  variety  of  crops  is  great. 

To  summarize,  it  appears  that  even  though  imagery  ol  several 
miles  effective  swath  width  (depth  in  range  for  this  system)  would  be 
desirable,  such  imagery  cannot  be  produced  at  appropriate  angles  with 
a  ground-based  system.  Relatively  narrow-swath  imagery,  however,  can 
be  produced  and  will  be  extremely  valuable  for  investigating  the  broad 
bandwidth  aspects  of  radar  imagery.  The  overall  system  range,  then,  to 
produce  this  type  of  imagery  need  only  be  in  the  order  of  3,000  feet. 
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From  other  considerations#  namely  the  mechanical  aspects  of  ob¬ 
taining  and  mounting  the  antenna  and  azimuth-elevation  positioner#  a  maxi 
mum  diameter  of  6  feet  is  placed  upon  the  reflector  size .  An  antenna  of 
this  size  should  be  procured  for  use  at  maximum  range.  As  mentioned  pre¬ 
viously#  another  antenna  of  smaller  size  and  larger  elevation  beam  angle 
is  desirable  for  short-range  imaging  at  lower  incidence  angles.  The  2  or 
3  foot  diameter  reflector  recommended  for  scattering  mode  measurements 
would  also  be  usable  for  this  application. 

The  foregoing  analysis  has  assumed  throughout  that  the  E  and  H 
plane  radiation  patterns  are  the  same.  This#  of  course#  is  not  necessarily 
the  case;  however#  it  is  proposed  that  the  feed  system  be  capable  of  both 
horizontal  and  vertical  polarization#  so  that  any  advantage  gained 
in  one  case  would  be  lost  when  operating  with  the  orthogonal  polarization. 

Before  drawing  final  conclusions  based  on  the  above#  the  required 
power  output  will  be  studied.  For  the  purpose  of  initial  calculations  a 
2-foot-diameter  paraboloid  will  be  assumed  for  scattering  and  near-range 
imaging  operation  and  a  6  foot  diameter  is  assumed  for  maximum-range 
imaging  operation.  A  tentative  maximum  range  of  500  feet  for  the  scatter¬ 
ing  mode  and  3 #000  feet  for  imaging  operation  is  assumed. 

4.  OUTPUT  POWER  CONSIDERATIONS 

Using  the  tentative  conclusions  of  the  preceding  sections#  the 
required  output  power  will  be  calculated  for  both  the  scattering  and  imag¬ 
ing  mode  of  operation. 

The  required  radiated  power  is  given  by  the  radar  equation  as 

T%  — _ r  _ 


(5) 


where 


R 

G 

X 


a 


=  transmitted  power 
=  minimum  required  received  power 
=  range 

=  antenna  gain 
=  wave  length 
=•  radar  cross  section 


The  minimum  required  received  power  may  be  expressed  as 


Pj.  =  K  T  B  F  S 


where  K 
T 
B 
F 
S 


=  Boltzman's  constant 
=  temperature  ^Kelvin 
=  received  bandwidth 
=  receiver  noise  figure 

=  minimum  required  output  signal-to-noise  ratio 


The  antenna  gain  may  be  expressed  as 


G  = 


47r  €  A, 


where  e  =  antenna  efficiency 
Ap  =  aperture  area 


The  radar  cross  section  may  be  expressed  as 
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where  a  =  radar  scattering  cross  section  per  unit  area 
o 

A  =  illuminated  area 

Using  the  above  relations,  the  required  radiated  power  niay  now 
be  expressed  as 


^  47r  K  T  B  F  S  (9) 

2.2. 

Some  additional  simplifications  may  still  be  made  by  noting 
that  we  are  primarily  interested  in  the  maximum  required  power.  First 
we  note  that 


A  =  ^  (1°) 

^p  4 


where  d  =  diameter  of  paraboloid 


The  radar  scattering  cross  section  will  be  a  minimum  for  composite 
type  targets/  and  may  be  expressed  as 


<7  = 


(11) 


where  =  scattering  cross  section  per  unit  area 
A  =  illuminated  area 

Now,  if  we  approximate  the  antenna  pattern  by  a  square  beam, 
the  illuminated  area  may  be  expressed  as 
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A  =  p^p,  (12) 

where  =  3  7.imuth  resolution 
=  range  resolution 

The  azimuth  resolution  is  provided  by  beamwidth  limitation  in 
all  operating  modes ,  thus 

Pa=R^  (13) 

For  the  purpose  of  calculating  the  maximum  required  power  we 
are  interested  in  near-maximum  range,  which  corresponds  to  high 
incidence  angles.  At  the  higher  incidence  angles,  or  near  maximum 
range,  the  range  resolution  is  provided  by  pulse-length  limitation  of  the 
illuminated  area  in  any  operating  mode.  Thus 


-  c  T 
^r  ~  2  sin  0 


(14) 


where  c  =  velocity  of  light 
T  =  pulse  length 
0  =  incidence  angle 


The  required  radiated  power  may  now  be  expressed  as 


128 

TT 


K  T  B  F  S  \  sin  0 
2  ,3 

e  d  c7q  c  T 


(IS) 


Some  additional  simplification  of  the  above  expression  is  possible 
with  the  following  assumptions. 
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The  input  noise  temperature  of  the  receiver  is  assumed  to  be 
at  the  standard  temperature  of  290°K  (approximately  62  F) .  The  factor 
KT  is  then 


9 1 

KT  »  4  X  10  watts/hertz  of  bandwidth 

The  receiver  bandwidth  is  taken  as  that  compatible  with  the 
design  goa.  for  minimum  pulse  length.  The  design  goal  for  maximum 
operating  pulse  length  is  taken  as  20  nsec.  The  compatible  bandwidth 

is  given  by 


B  ==  Lil  =  =  60  MHz  (16) 

2  X  10® 

Since  we  are  interested  primarily  in  the  maximum  required  power, 
the  wavelength  is  taken  as  the  maximum  across  the  operating  range 

\  (f  =  4  GHz)  =  “  =  0.075  meters  (17) 

Again,  since  we  are  interested  primarily  in  the  maximum  range 
condition  we  assume 

sin  0  ss  1 


The  antenna  efficiency  is  taken  from  manufacturers'  data .  The 
efficiency  figure  used  is  the  lowest  value  given  for  any  feed-reflector 
combination  considered  for  this  application.  Assume 

c  =  0.3 


Of  course  c,  the  velocity  of  light,  is  taken  as 
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g 

c  =  3  X  10  meters/sec 

To  summarize,  the  assumptions  made  regarding  the  parameters 
of  Equation  (15)  are 

—2 1 

KT  =  4  X  10  watts/hertz  of  bandwidth 
•^mln  =  20  nsec 

B  =  =  60  MHz 

T 

\  =  0«075  meters 

sin  0=1 


c  =  0.3 

g 

c  =  3  X  10  meters/sec 
Equation  (15)  may  then  be  expressed 


P^=  2.72  X  10 


-22  F  S  R' 


d  OTqT 


If  the  minimum  pulse  length  is  used,  we  have 


.36  X  10 


-12  F  S  R' 


a  O', 


(18) 


(19) 
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Th6  exprsssion  abovG  Is  for  th©  power  required  at  the  antenna 
port,  in  terms  of  the  power  received  at  the  same  port,  and  does  not  take 
into  account  the  internal  system  losses.  To  arrive  at  the  figure  for  the 
power  required  from  the  power  amplifier  we  must  make  some  assumptions 

regarding  these  internal  losses. 

Referring  to  the  block  diagram  of  Figure  1 ,  we  shall  assume  a 
duplexing  loss  of  2  db  (each  way)  and  a  cabling  and  mismatch  loss  of 
1.5  db  for  both  transmit  and  receive.  The  total  system  internal  losses 

are  then 


L=2. 0  +  2.0  +  1.5  +  1.5=7.0db 

The  power  required  from  the  power  amplifier  is  then  increased  by 
this  ratio,  and  is  given  by 


=  6.8  X  10 


-12  F  S 
a 


(20) 


4,1  Scattering  Mode  of  Operation 

From  the  discussion  of  section  3,1  we  make  the  following  assumptions 
for  calculating  the  power  required  in  the  scattering  mode . 

F  =  10  db 


S  =  20  db 


R  =  500  feet  =152  meters 
max 


d  =  2  feet  =0.61  meters 


(7  =  -40  db 

o 


m 
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With  these  assumptions,  the  required  amplifier  power  given  by 
Equation  (20)  is 

3 

P  =6.8x10”^^  — =  1050  watts 

a 

The  survey  of  traveling  wave  tube  amplifiers  for  use  in  this 
application  has  shown  the  available  amplifiers  to  fall  into  three  general 
power  ranges.  There  are  a  wide  variety  of  tubes  available  in  the  range 
of  1  to  20  watts  CW  power  output.  The  price  range  for  these  tubes  runs 
from  approximately  $2,500  to  $5,500  including  power  supply.  The  next 
range  is  100  watt  CW  power  output, and  the  average  cost  of  these  tubes 
is  $7,000  less  power  supply.  The  highest  power  available  the  desired 
one-octave  frequency  range  is  1  KW.  These  tubes  are  pulse  modulated 
with  approximately  a  1%  duty  cycle  and  alp.  sec  minimum  pulse .  The 
cost  of  these  tubes  is  in  general  over  $10,000  without  pov;er  supply. 

Comparing  these  prices  with  the  THEMIS  budget,  we  see  that 
only  the  low-power  CW  tubes  are  within  the  price  range  allotted  for  this 
component.  This  makes  it  extremely  desirable  to  operate  in  this  power 
range  if  at  all  possible.  There  is  yet  another  consideration  making  this 
still  more  des  Table.  If  a  1  KW  tube  is  used, it  cannot  be  driven  directly 
from  the  mixer  output  of  the  system, as  the  B.  W.  O.  sweep  generator 
has  a  leveled  output  of  10-15  MW  and  the  conversion  loss  of  the  mixer 
is  approximately  10  to  13  db.  This  means  an  intermediate  stage  of  wide¬ 
band  amplification  must  be  provided.  The  lower  power  tubes  (1-20  watts) 
have  enough  saturated  gain  that  they  may  be  driven  directly  without  an 
intermediate  stage.  It  would,  then,  be  extremely  desirable  to  initially 
operate  the  system  with  a  lower  power  tube.  If  more  power  is  required 
at  a  later  date  for  either  ground-based  or  possibly  airborne  operation,  the 
higher  power  output  could  be  added  with  the  previous  low-power  output 
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amplifier  serving  as  a  driver. 

Let  us  again  examine  the  required  power  output  to  see  if  it  may 
be  reduced  from  the  1050  watt  level  shown  previously. 

The  output  signal- to-noise  ratio  of  20  db  may  be  obtained  by 
postdetection  integration  of  several  pulses.  The  improvement  over  the 
predetection  signal-to-noise  ratio  decreases  the  required  power  by  the 
same  factor.  This  averaging  or  integration  is  not  only  desirable  for 
power  output  reduction,  but  is  actually  required  ii  tne  differential  scat¬ 
tering  cross  section  of  composite  targets  is  to  be  measured  without  signi¬ 
ficant  error  due  to  Rayleigh  fading  of  the  return  signal. 

The  radar  return  signal  from  an  area  target  is  made  up  of  component 
signals  re-radiated  from  numerous  small  scatterers  and  simultaneously 
observed  by  the  radar.  The  average  power  received  is  the  sum  of  the 
powers  received  from  the  individual  scatterers.  The  differential  scatter¬ 
ing  cross  section  is  then  an  average  and  can  only  be  calculated  when 
enough  measurements  to  compute  a  good  avera'jo  have  been  made  at  points 
where  the  signals  are  essentially  independent  of  each  other. 

The  fading  nature  of  the  return  signal  described  above  is  evident 
whenever  there  is  relative  motion  between  the  target  and  radar.  In  this 
ground-based  operating  mode,  there  is  no  constant  velocity  difference  as 

t 

with  airborne  systems.  When  measuring  a  single  discrete  target  such  as 
a  car  or  tnick  there  will  be  essentially  no  effect.  However,  when  measur¬ 
ing  a  composite  target  such  as  a  wheatfielc^  the  movement  caused  by  the 
wind  will  produce  a  widely  varying  return  signal.  If  the  effects  of  this 
variation  are  to  be  overcome, the  return  signal  must  be  averaged. 

If  we  assume  the  return  to  be  made  up  of  many  components  with 
relatively  the  same  amplitude  and  random  phase,  the  return  will  be 
characterized  by  a  Rayleigh  distribution.  This  distribution  is  characterized 
by  equal  mean  and  standard  deviation:  that  is,  the  distribution  after  square 
law  detection  is  given  by 
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V_ 

,  "''o 

p(V)=f-®  (21) 

o 

o  =  B  =  Vg 


where  ji  is  the  mean,  V  is  the  postdetection  voltage,  and  a  is  the  standard 
deviation . 

For  relatively  small  numbers  of  independent  samples,  the  standard 
deviation  for  the  mean  of  Rayleigh-distributed  variables  differs  somewhat 
from  that  for  normally-distributed  variables .  As  the  number  of  samples 
Increases,  the  central  limit  theorem  applies, and  Rayleigh  and  normal 
results  are  very  close.  The  standard  variation  of  the  mean  is  given  approxi¬ 
mately  by 


(22) 


Forn  larger  than  20  or  30,  this  approximation  is  quite  good.  For  smaller 
values  of  n,  it  is  off,  both  because  of  the  necessity  to  apply  small-sample 
theory  and  because  of  the  difference  between  the  Rayleigh  distribution 
and  the  Gaussian. 

The  5%  and  95%  range  that  the  mean  can  take  on  is  given  by 


Range  =  2x1. 645 


3 .29  u 


(23) 


when  the  assumptions  of  Equation  (22)  are  made.  Expressed  in  decibels, 
this  range  is  plotted  in  Figure  5.  Clearly,  a  range  of  2  or  3  db  is  easy 
to  achieve,  but  a  range  much  better  than  1  db  requires  a  very  large  number 
of  independent  samples.  The  experiment  design  must  take  this  into  account. 
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and  compromise  is  necessary  between  the  integration  time  allowed  and  the 
desired  range  for  the  mean  value. 


NUMBER  OF  INDEPENDENT  SAMPLES  (n) 


Figure  5  -  Accuracy  of  Averages 

The  required  independence  of  the  samples  to  be  averaged  may 
arise  from  variation  in  either  time  or  frequency.  The  individual  and 
composite  effects  of  both  these  types  of  averaging  will  be  studied  in 
a  later  Memorandum.  We  will  concern  ourselves  here  only  with  determin¬ 
ing  the  number  of  independent  samples  desired.  Examination  of  Figure  5 
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shows  it  to  take  200  independent  samples  to  reduce  the  range  of  the 
mean  to  1  db.  We  will  assume  this  to  be  the  desired  range. 

With  this  assumption  that  200  independent  samples  will  be 
averaged/ we  may  take  advantage  of  the  integration  improvement  of  signal- 
to-noise  ratio.  If  we  assume  postde^^ection  integration^ the  single¬ 
pulse  signal-to-noise  ratio  is  given  by 


s 


1 


(24) 


where  Sj^  =  single  pulse  signal-to-noise  ratio 


s  =  output  signal-to-noise  ratio 
n  =  number  of  pulses  integrated 

X  =  0.5  to  0 . 9, depending  on  integration  efficiency  assumed 


If  we  assume  x  =  0.75  and  n  =  200 


n*=  (200) 53 

=  17.2  db 

Referring  again  to  Table  2,  we  note  that  the  pulse  length  may  be 
increased  by  a  factor  of  2  at  the  higher  incidence  angles  (above  75  de¬ 
grees)  while  maintaining  a  range  resolution  equal  to  or  less  than  the 
azimuth  resolution.  If  we  use  a  40  nsec  pulse  at  the  higher  incidence 
angles, another  3  db  improvement  over  the  result  obtained  from  Equation  (20) 
may  be  realized.  Applying  this  20.2  db  improvement  to  the  1050  watt 
figure  obtained  previously#  we  find  the  power  required  from  the  final  amplifier 
to  be 


=  1050  X 


1 

105 


=  10  watts 


B  28 

To  summarize  again,  with  a  10  watt  final  amplifier  we  may 
measure  the  differential  scattering  cross  section  to  an  incidence  angle 
of  at  least  85  degrees  and  down  to  a  level  of  -40  db.  This  may  be 
accomplished  with  a  2  foot  diameter  paraboloid  and  a  relatively  inefficient 
broadband  feed. 

4.2  Imaging  Mode  of  Operation 

From  the  discussion  of  section  3.2  we  make  the  following  assump¬ 
tions  for  calculating  the  power  required  in  the  imaging  mode . 

F  =  10  db 

S  =  20  db 

Rmax  ~  3,000  feet  =  915  meters 
d  =  6  feet  =1.83  meters 

=  -30  db 

With  these  assumptions,  the  required  amplifier  power  given  by 
Equation  (20)  is 


P  =  6.8  X  10"^^  =  850  watts 

Again^at  the  more  distant  ranges  we  may  use  a  pulse  length  of  at 
least  twice  the  minimum,  which  reduces  the  required  power  to  425  watts. 

Assuming  the  same  postdetection  integration  efficiency  as  before, 
we  must  integrate  the  following  number  of  pulses  to  operate  with  10  watts: 
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0.75  _  425 
"  "10 


42.5 


n=  150 


This  amount  of  integration  may  easily  be  provided  by  slowing  the 
azimuth  sweep  rate.  Some  trade-off  between  signal-to-noise  ratio  and 
the  assumed  value  for  differential  scattering  cross  section  is  also  possible 
and  will  be  studied  in  a  later  memorandum. 

Summarizing  once  again,  it  appears  feasible  to  meet  the  imaging 
requirements  with  a  power  amplifier  output  in  the  range  of  10  watts. 

5.  COMPONENT  SPECIFICATIONS 

A  thorough  survey  of  available  broadband  antenna  and  TWT 
amplifiers  has  been  conducted.  The  available  components  are  now 
evaluated  in  light  of  the  foregoing  analysis. 

5.1  Antenna 

The  initial  survey  of  antenna  manufacturers  has  been  completed. 

On  the  basis  of  price,  delivery  and  completeness  of  line, consideration 
has  been  reduced  to  two  companies;  American  Electronics  Laboratories 
and  Scientific. Atlanta. 

Before  comparing  specific  model  numbers  the  antenna  requirements 
will  be  summarized. 

1.  Frequency  range;  4, 0-8.0  MHz  minimum 

2 .  Polarization;  linear 

The  feed  system  should  be  capable  of  either  vertical  or  horizon¬ 
tal  polarization. 

3.  Reflector  size;  one  6  foot  diameter  paraboloid 

one  2  or  3  foot  diameter  paraboloid 
A  rectangular  aperture  would  be  preferable  for  imaging  operation; 
however,  it  appears  that  none  is  available  at  this  time. 


4.  Movement  range:  -2°  to  +90^  elevation 

+  90^  azimuth 

The  azimuth  movement  must  be  capable  of  sector  scanning. 

Elevation  movement  need  merely  be  remotely  adjustable. 

Both  of  the  firms  under  consideration  have  logarithmic  feeds  with 
linear  polarization  and  a  frequency  range  of  1-11  GHz.,  The  polarization 
of  this  type  feed  is  variable  by  adjustment  of  the  feed  mounting .  Scientific- 
Atlanta  also  has  a  line  of  dual-polarized  feeds.  These  feeds  are  ridged 
waveguide  horns  and  are  available  with  frequency  ranges  compatible  with 
available  sweep  equipment. 

The  choice  between  the  types  of  feeds  is  basically  one  of  trading 
frequency  coverage  for  the  ability  to  measure  the  depolarized  component. 

The  increased  range  (over  just  the  4. 0-8.0  GHz)  is  desired^as  the  system 
will  be  expanded  in  Phase  II  to  cover  at  least  one  additional  band.  With 
the  logarithmic  type  feed  the  same  antenna  could  be  used.  With  the  horn 
type  feed  the  antenna  would  have  to  be  changed  along  with  the  other  RF 
equipment  when  changing  operating  bands.  However,  the  dual  polarization 
of  the  feed  makes  it  possible  to  easily  change  linear  polarization  and  to 
measure  the  depolarized  return.  With  the  addition  of  matched  feed  lines  and 
a  90°  broadband  hybrid  it  is  even  possible  to  transmit  and  receive  circular 
polarization 

Shown  below  are  three  (Comparative  antenna  models. 


Table  5  -  Comparison  of  Antenna  Models 


Manufacturer 

Freq.  Range 

Reflecter  Dia. 

Type  Feed 

Price 

AEL 

1-11  GHz 

6  ft. 

Logarithmic 

$  995 

S-A 

4.0-12.4  GHz 

6  ft. 

Logarithmic 

$1,670 

S-A 

3.95-8.1  GHz 

6  ft.. 

Dual  Polarized 

$1,825 

Horn 
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The  AEL  unit  does  not  include  either  mounting  hardware  or 
feed  line^i  as  do  both  Scientific-Atlanta  units .  The  AEL  unit  may  be  pro¬ 
vided  with  a  remote-control  feed  mounting  adjustment  for  changing  polari¬ 
zation.  This  is,  of  course,  at  additional  cost. 

Primarily  on  the  basis  of  ability  to  measure  the  depolarized  returrv 
and  secondarily  because  of  the  ease  of  switching  polarization,  the  Scientific 
Atlanta  dual-polarized  horn  type  feed  system  is  selected. 

Table  6  shows  a  listing  of  the  antenna  equipment  recomriended  for 
purchase. 


Table  6  -  Recommended  Antenna  Components 


Manufacturer 

Item 

Model  No. 

Quantity 

Price 

S-A 

6  ft.  dia.  Parabolic 
Reflector 

22-6 

1 

$1 

,050 

S-A 

2  ft.  dia.  Parabolic 
Reflector 

22-2 

1 

$ 

350 

S-A 

Dual  Polarized  Horn 
Feed 

28-3.9/6 

1 

$ 

775 

S-A 

Dual  Polarized  Horn 
Feed 

28-3.9/ 

1 

$ 

775 

S-A 

Elevation-over-Azimuth 

Positioner 

5403-1-E 

1 

$3,190 

S-A 

Positioner  Control 

Unit 

4112 

1 

$ 

495 

S-A 

Local  Positioner 

Control 

4112L 

1 

$ 

250 

S-A 

Slip  Ring  Assembly 

-16 

1 

$ 

500 

S-A 

Control  Cables  (60ft) 

5051-12 

2 

$ 

210 

S-A 

Control  Cables  (60ft) 

5051-8 

1 

$ 

92 

S-A 

Indicator  Unit 

4422-12 

1 

$ 

590 

S-A 

Rotary  Joint 

10-1 

1 

300 

$8,577 

An  additional  $700  will  be  spent  for  factory  alignment  of  feeds 

and  patterning  of  each  unit  at  3  frequencies.  $  700 

$9,277 
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5.2  Transmitter 

As  shown  previously,  the  required  power  is  in  the  range  of  10 
watts.  The  power  ranges  and  approximate  price  ranges  of  TWT  amplifiers 
(tube  and  power  supply)  are  shown  in  Table  7 . 

Table  7  -  TWT  Power  Ranges 


Power 

Operation 

Price 

1-20  watt 

CW 

$  2,500  -  $  5,500 

100  watt 

CW 

$10,000  -  $15,000 

100  watt 

pulse 

$10,000  -  $15,000 

1  kilowatt 

pulse 

$15,000  and  up 

From  the  above  it  is  seen  that  the  10*watt  range  not  only  suits 
the  system, but  is  also  the  only  range  that  likewise  suits  the  budget. 
During  Phase  III  the  system  will  become  airborne  and  will  then,  in  all 
probability,  require  a  1  kilowatt  power  amplifier.  The  proposed  10  or 
20  watt  amplifier  for  the  ground-based  system  will  serve  quite  nicely  as 
the  driver  for  a  later  airborne  system.  At  the  upper  end  of  the  price  range 
shown,  it  is  possible  to  obtain  a  20 -watt  amplifier  with  sufficient  gain 
to  be  driven  directly  by  the  swept  oscillator  or  mixer  output.  The  lower- 
priced  amplifiers  are  in  reality  not  adequate  by  themselves,  for  most 
have  a  saturated  gain  of  30  db  or  less, which  is  not  sufficient  for  this 

application  without  a  low-power  driver  stage. 

The  TWT  amplifier  recommended  for  purchase  is  the  Alto  Scientific 
Company  Model  20C4, 0-8. 0-Q60.  This  unit  has  a  maximum  output  of 
20  watts  CW  with  a  saturated  gain  of  50  db.  The  cost  of  the  unit  (TWT 
and  power  supply)  is  $5,055, 

Acknowledgement;  This  work  was  supported  by  Project  Themis  (USAETL 
Contract  DAAK  02-68-C-0089 ,  ARPA  Order  No.  1079). 
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Introduction 


The  IDECS  System  consists  of  real-time  and  near  real-time  image 
processing  equipment.  The  individual  devices  in  the  IDECS  System  are 
designed  primarily  for  versatility  of  operation  and  interconnection  to  allow 
evaluation  of  different  image  processing  techniques.  Th^  basic  capabilities 
of  the  IDECS  System  have  been  added  under  various  research  contracts  with 
the  principal  contributions  as  follows: 

1.  Real-time  video  processing  devices,  color  matrixing,  and 
Initial  research  were  supported  primarily  under  NASA  contract 
#  17-004-003. 

2.  A  six  image  synchronous  flying  spot  scanner  system  and 
general  purpose  real-time  decision  logic  with  color  coding 
capabilities  (the  spectrance  selector)  are  being  supported 
under  USAETL  contract  #  DAAK02-67-C-0435 . 

3.  Improved  operator  control,  digital  interface  and  the  capabilities 
for  automatic  control  card  processing  are  being  supported  under 
project  Themis. 

The  final  configuration  of  the  IDECS  System  as  it  is  now  envisaged  is  shown 
in  Figure  1.  The  subsystems  which  are  being  supported  by  Project  Themis 
are  indicated  by  asterisks  and  are  described  individually  in  the  text. 

I_.  .  ANALOG -TO -DIGITAL  CONVERSION  SUBSYSTEM 

A.  Function.  The  function  of  this  subsystem  is  to  convert  analog 
signals  from  the  scanning  densitometer  or  any  of  the  flying-spot  scanners 
into  digital  information.  Each  analog  value  is  converted  to  a  binary  number 
variable  between  six  and  twelve  binary  digits  and  then  stored  in  one  or  two 
6-bit  words  on  magnetic  tape.  This  subsystem  is  also  to  have  sufficient 
flexibility  so  that  several  channels  can  be  multiplexed  together.  This  would 
allow  the  repetitive  sampling  of  each  flying-spot  scanner  in  a  sequence. 


* 

Acronym  for  Image  Discrimination,  Enhancement 
and  Combination  Svster: 
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B.  Status.  The  analog-to-digital  converter  and  the  magnetic  tape 
unit  are  on  older  and  should  be  received  next  month.  Unfortunately,  the 
orders  were  delayed  for  some  time  because  of  the  involved  procedure  required 
to  obtain  approval  for  digital  equipment.  The  control  logic  and  circuits  have 
been  designed  for  maximum  flexibility  and  are  currently  being  implemented. 

The  analog-to-dlgital  conversion  subsystem  is  currently  scheduled  for  com¬ 
pletion  by  July  1 ,  1968,  and  at  that  time  a  detailed  report  will  be  written. 

A  block  diagram  of  this  subsystem  is  shown  in  Figure  2 . 

II.  DIGITAL-TO-ANALOG  SUBSYSTEM 

A.  Function.  The  function  of  this  subsystem  is  to  convert  digital 
information  from  a  digital  magnetic  tape,  a  small  computer,  or  other  digital 
information  sources  into  an  analog  voltage.  Initially  this  subsystem  will 

be  able  to  convert  two  binary  numbers  of  ten  bits  each  into  two  analog  voltages 
suitable  for  input  into  a  CRT  display. 

B.  Status .  The  two  required  dlgital-to-analog  modules  have  been 
ordered  and  received.  The  necessary  control  circuits  are  being  designed  and 
implemented.  This  subsystem  is  scheduled  for  completion  by  August  1,  1968, 
and  at  that  time  a  detailed  report  will  be  submitted.  The  block  diagram  for 
this  subsystem  is  shown  in  Figure  3. 

III.  DISPLAY  INTERFACE  SUBSYSTEM 

A.  Function.  The  function  of  this  subsystem  is  to  enter  data  at  a  slow 
rate  into  the  video  memory  and  then  display  this  information  at  video  rates 
onto  the  color  display  unit  or  any  black-and“'V/hite  monitor .  Also  this  sub¬ 
system  will  eventually  have  the  capability  of  sampling  the  video  signal  to 
produce  an  analog  signal  in  the  audio  range  which  can  be  transmitted  directly 

over  telephone  lines. 

The  magnetic  disc  video  '■ecordlng  subsystem  will  be  capable  of  the 
following  functions. 

1.  Simultaneous  playback  of  six  s.eparate  signals 

2.  Simultaneous  "real-time"  recording  or  erasure  of  three  data 
tracks ,  this  function  to  be  mechanically  switchable  to  a 
second  group  of  three  data  tracks. 
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3.  Sequential  "slow-scan"  recording  of  any  or  all  of  the  above 
data  tracks,  again  by  means  of  mechanical  switching. 

4.  Provision  for  the  use  of  a  suitable  electrical  input  signal  to 
accomplish  the  functions  of  selective  recording  and  erasure  of 
data. 

5.  A  TV  keying  generator  with  output  signal  compatible  with  the 
requirements  of  (4)  and  generating  a  rectangular  pattern  (as  seen 
on  a  TV  monitor  screen)  of  variable  dimensions  and  position. 

The  same  generator  will  supply  separate,  variable  delay,  vertical 
and  horizontal  drive  outputs. 

6.  A  video  converter  which  accepts  a  standard  "real-time"  television 
input  signal  and  uses  high-speed  sampling  to  produce  a  "slow- 
scan"  output  in  the  audio  frequency  range. 

A  block  diagram  of  this  subsystem  is  shown  in  Figure  4. 

Status.  An  extensive  study  was  made  to  determine  the  best 
method  of  storing  at  least  two  images  with  reasonable  cost.  The  result  was 
that  there  is  no  commercial  unit  readily  available  that  would  satisfy  all 
requirements;  howrver,  one  equipment  manufacturer  is  willing  to  modify  an 
existing  product  to  meet  most  of  our  requirements.  The  video  memory  order 
has  been  placed  but  it  will  not  be  received  unvil  next  fall.  Some  of  the 
control  circuits  have  been  designed,  although  none  as  yet  have  been  built. 
Completion  date  for  this  subsystem  is  tentatively  set  for  January  1,  1969. 

IV.  SMALL  COMPUTER  INTERFACE  SUBSYSTEM  (Optional.  To  be  included 
if  funding  for  the  computer  is  obtained  either  ough  THEMIS  or  other  sources). 

A.  Function.  The  purpose  of  this  subsystem  is  to  interface  a  compu¬ 
ter  so  that  data  may  be  processer  and  various  control  functions  performed  in 
"real  time. " 

B.  Status .  The  analog-to-digital  conversion  subsystem  and  the 
digital-to-analog  subsystem  are  being  designed  so  that  only  a  small  addi¬ 
tional  cost  will  be  incurred  if  these  units  are  connected  directly  to  a  small 
computer.  Also,  techniques  have  been  investigated  for  coding  data  so  that 
one  frame  of  video  information  at  reduced  resolution  could  be  stored  in  an  8K 
computer  memory  by  utilizing  the  redundancy  in  the  data. 
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V.  SCANNING  DENSITOMETER 

A.  Function.  To  process  an  image  on  a  digital  computer,  it  is 
required  that  the  image  information  be  changed  to  an  analog  signal.  This 
can  be  done  by  using  a  scanning  densitometer  which  consists  of  a  rotating 
and  translating  glass  drum  with  a  photomultiplier  tube  inside.  The  image 
is  placed  around  the  glass  drum.  A  light  source  placed  above  the  drum 
shines  through  the  image  into  the  photomultiplier  tube.  The  output  of  the 
photomultiplier  is  an  analog  signal  proportional  to  the  grey  scale  of  the 
image.  The  translation  of  the  image  and  the  spot  size  are  set  so  that  20 
line  pairs/mm  are  produced. 

B.  Status.  The  basic  unit  has  been  designed  and  built;  however, 
the  control  circuits  have  not  been  designed. 

VI.  MASTER  VIDEO  CONTROL 

A.  Function.  To  evaluate  the  functions  and  subsystems  individually 
and  combinationally ,  it  is  necessary  that  the  IDECS  system  have  the  capabi¬ 
lity  to  connect  and  disconnect  any  device.  This  is  the  function  of  the  video 
control.  By  switching  the  video  control,  various  combinations  of  input  and 
processing  devices  are  presented  at  the  output  devices.  The  IDECS  block 
diagram  (Figure  1)  shows  the  input,  processing,  and  output  devices  that  are 
connected  to  the  video  control. 

B.  Status .  The  video  control  has  been  designed  and  constructed  but 
has  not  been  connected  into  the  system. 

VII.  POSITION  AND  SIZE  DEVICES 

A.  Function.  In  order  to  combine  and  process  more  than  one  image, 
it  is  necessary  that  the  images  be  geometrically  congruenced,,  This  can  be 
accomplished  by  controlling  the  raster  on  the  flying  spot  scanning  units . 

The  position  and  size  devices  control  this  raster. 

The  positioning  device  will  translate  and  rotate  the  raster,  causing 
the  image  cn  the  display  to  move  horizontally,  vertically,  and  angularly. 

The  size  device  will  expand  or  contract  the  raster;  therefore,  the  size  of 
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the  image  on  the  display  is  controlled  horizontally  and  vertically.  The  posi¬ 
tion  and  size  devices  are  shown  in  Figure  5. 

B.  Status .  Both  the  position  and  size  devices  have  been  implemented 
for  all  six  flying  spot  scanning  units. 

VJII.  SYNCHRONIZATION  DEVICE 

A.  Function.  For  the  information  on  the  input  and  output  devices  to 
correspond,  it  is  essential  that  these  devices  scan  in  unison.  The  synchro¬ 
nization  device  accomplishes  this  by  sending  pulses  to  the  horizontal  input 
of  the  sweep  circuits  at  a  rate  of  15,750  Hz  and  to  the  vertical  input  of  the 
sweep  circuits  at  a  rate  of  60  Hz.  The  synchronization  device  is  shown  in 
Figure  5. 

B.  Status .  A  synchronization  generator  is  in  operation  to  match  the 
flying  spot  scanning  units  to  the  displays . 

IX.  MONITOR  DISPLAY 

A*  Function.  The  monitor  display  will  exhibit  any  input  or  processing 
device  when  connected  properly  through  the  master  video  control.  Since  it 
will  be  mounted  in  the  operator’s  console,  it  will  be  convenient  for  use  in 
congruencing  several  images .  It  will  also  be  used  in  calibrating  the  flying 
spot  scanning  unit  and  other  circuits . 

B.  Status .  A  10"  cathode  ray  tube  capable  of  20  MHz  has  been 
ordered.  A  Tektronix  oscilloscope  is  operational  as  a  monitor  display  at 
present. 

X.  COLOR  DISPLAY 

A.  Function.  The  color  display  is  used  to  convert  video  signals  pro¬ 
portional  to  grey  scale  on  photographs  or  radar  images  into  a  color  which 
varies  in  intensity  and  hue. 

B.  Status  .  The  color  display  is  operational. 
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ABSTRACT 

Automatic  classification  of  remote  sensed  data 
It  a  necessity  since  satellite  remote  sensors  currently 
on  the  design  board  are  expected  to  gather  more  data  each 
hour  than  could  be  analyzed  by  traditional  methods  In  a 
year.  A  machine  which  could  determine  the  structure  of 
the  sensed  envlronmert  by  grouping  together  similar  data 
signals  and  relabeling  these  signals  with  the  same  label 
would  reduce  both  storage  requirements  and  transmission 
times  for  the  remote  sensor  In  addition  to  doing  pattern 
recognition  In  space  or  on  the  ground. 

The  paper  Is  Introduced  by  a  discussion  concern¬ 
ing  the  adaptive  process  of  how  predictive  models  of  an 
environment  can  be  generated.  This  discussion  leads 
Into  a  mathematical  description  of  such  a  model  for  pattern 
recognition.  The  model  was  programmed  for  the  GE  625 
and  four  possible  combinations  of  transmitted  and  received 
like  and  complementary  polarized  K-band  radar  images  of 
agriculture  In  western  Kansas  were  used  as  data. 

Results  confirm  that  various  categories  of  vegetation 
such  as  alfalfa,  grain  sorghum,  wheat  stubble  and  com  can¬ 
not  be  distinguished  from  one  another  on  the  basis  of  their 
structure  within  the  radar  data.  This  suggests  that  for 
pattern  recognition  purpose  t,  certain  categories  may  be  expanded 
with  little  or  no  Information  loss  and  an  Increase  in  efficiency 
of  the  pattern  recognition  process.  Further  It  also  suggests 
*  that  If  the  Investigator  wants  to  distinguish  between  the 

confused  categories  such  as  alfalfa,  grain  sorghum,  wheat 
stubble  and  com,  an  additional  sensor  must  be  used  which 
'  can  make  the  distinction. 


'  Non-lmaglng  sensor  systems  such  as  gelger  counters  magnetometers,  gravity  gradlometers , 
spectrometers,  photometers,  and  particle  counters  at  present  are  collectively  gathering  a  quarter 
billion  bits  of  Information  per  day.  Some  planned  hlgh-resolution  imaging  systems  are  expected 
together  by  themselves  over  three  billion  bits  of  Information  per  day.  This  makes  automatic 
pattern  recognition  of  remote-sensed  data  a  necessity.  Over  one  hundred  thousand  analog  tapes 
of  remote  sensed  data  are  stored  in  ESSA's  library  and  about  two  hundred  additional  such  tapes 
are  created  dally.  The  only  way  such  a  great  volume  of  information  can  be  utilized  is  by  auto¬ 
matically  recognizing  the  pattern  of  relationships  which  exist  within  the  data. 


A  pattern-recognition  system  Is  one  kind  of  system  which  determines  the  relatlonshl^p  of 
category  cr  type.  We  will  endeavor  to  describe  an  adaptive  pattern-recognition  system,  which 
needs  no  teacher,  is  easily  built,  and  Is  able  to  automatically  classify  data.  The  system  we 
describe  has  been  simulated  and  tested  with  agricultural  radar  images,  and  the  results  of  the  test 
will  follow  the  description. 


A  pattern-recognition  system  takes  patterns  fgr  inputs  and  produces  outputs  which  are  either 
direct  symbols  for  categories  or  binary  M-tupIes  which  form  a  code  for  categories.  The  Input 
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patterns  form  N-tuple  of  measurements  (or  sometlri.es  their  binary  code)  taken  of  some  environment. 
Each  component  of  the  N-tuple  could  represent  the  output  of  one  sensor  system.  Figure  1  Illustrates 
the  pattern  recognition  process. 

We  begin  with  a  few  elementary  definitions.  A  system  Is  a  structured  device  which  produces 
•n' output  related  to  Its  past  and  present  Inputs  according  to  some  deterministic  or  probabilistic 
rule.  The  rule  Is  called  the  structure  of  the  system.  If  the  output  Is  related  only  to  the  present  In¬ 
put,  the  system  Is  memoryless.  If  the  output  Is  related  to  both  the  past  and  present  Input,  the  sys¬ 
tem  has  memory.  If  the  form  of  the  rule  at  time  t  depends  on  some  of  the  Inputs  before  time  t,  then 
the  system  Is  changeable;  otherwise,  the  system  Is  fixed. 

An  example  of  a  memoryless  fixed  system  could  be  a  radio,  the  electromagnetic  field  at  Its 
antenna  being  considered  the  Input  and  the  acoustical  field  generated  by  Its  loudspeaker  being  con¬ 
sidered  the  output.  It  Is  memoryless  because  the  output  does  not  depend  on  previous  Inputs: 
only  on  the  present  Input.  Past  radio  programs  do  not  Influence  the  radio's  present  operation.  An 
example  of  a  fixed  system  with  memory  could  be  an  adding  machine,  the  buttons  pushed  on  the  key¬ 
board  register  being  considered  the  Input  and  the  number  displayed  on  the  carriage  register  or  paper 
tape  be, mg  considered  the  output.  The  adding  machine  has  memory  of  the  sum  of  all  the  past  Inputs. 
Examples  of  changeable  systems  can  be  found  In  any  of  the  self-learning,  self-organizing  or  adap¬ 
tive  systems . 

A  pattern-recognition  system  can  either  have  a  memory  or  be  memoryless,  and  can  be  fixed 
or  changeable.  If  the  system  changes  Its  own  structure  lhv‘  change  will.  If  It  Is  to  be  meaningful, 
depend  on  the  Input-pattern  sequence,  so  that  the  present  state  of  structure  will  provide  some  Indi¬ 
cation  of  past  Input  history;  hence,  a  self -changeable  sysftem  Is  a  system  with  memory.  Fixed 
pattern-recognition  systems  have  limited  use,  since  the  precise  transformation  from  Input  pattern 
to  category  type  must  be  known  beforehand.  With  the  present  state  of  the  art  for  Imaging  sensors, 
prior  knowledge  Is  a  strict  requirement.  Radar  Images,  for  InsUnce,  can  be  calibrated  at  best  to 
only  3  db  from  the  Initial  transmitted  signal  to  the  final  developed  image.  Different  pictures  taken 
with  the  same  sensor  system  of  exactly  the  same  environment  at  two  different  times  can  vary  by  a 
two-to-one  gray  scale  variation.  A  fixed  pattern -recognition  system  would  have  little  use  here, 
because  the  data  has  an  unknown  bias.  All  of  Its  information  Is  relative  to  the  bias  levels.  Thus, 
what  Is  needed  Is  some  kind  of  pattern  recognition  system  which  can  change  or  adapt  according  to 
the  biases  and  stochastic  non-llnearltles  of  the  sensor  system  and  according  to  the  structure  of 
the  environment  sensed,  so  that  the  structure  of  the  pattern-recognition  system  Is  optimum  for  the 
kind  of  Job  It  Is  supposed  to  do. 

The  set  of  all  logically  possible  measurements  we  will  call  measurement  space.  To  each 
point  In  measurement  space  a  pattern-recognition  system  assigns  one  category,  so  that  the  set  of 
points  which  are  assigned  to  each  category  forms  a  cell  in  a  partition  of  measurement  space,  as 
Illustrated  in  Figure  2. 

The  pattern-recognition  problem  Is  usally  presented  as  follows:  The  Investigator  forms  I 
training  sequences  S^,  I  =  1 , . .  .1  of  measurements  such  that  the  measurements  In  the  ith  sequence 

all  represent  measurements  of  objects  which  the  Investigator  knows  are  In  the  I  category.  The  I 

sequence  actually  Is  a  sample  from  the  l^^  population,  which  Is  a  population  of  objects  In  the  I 
category.  It  Is  possible  for  the  same  measurements  to  occur  In  several  sequences  so  that  It  Is  a 
problem  to  decide  what  category  to  assign  to  such  measurements.  Given  the  sequence  S^,  1  =  1,...] 

the  Investigator  would  like  to  find  a  decision  rule  (equivalence  relation)  which  partitions  measure¬ 
ment  spacf'  Into  I  cells  so  that  the  l^^  cell  contains  the  set  of  measurements  which  are  most  repre¬ 
sentative  of  the  l^^  category.  A  pattern-recognition  system  Is  a  system  which  has  such  a  decision 
rule  for  Its  structure.  Since  the  structure  of  the  system  Is  determined  from  the  training  sample  we 
may  say  that  the  system  has  memory  relative  to  the  training  sample. 

With  the  problem  posed  In  this  way,  there  is  no  guarantee  that  the  categories  which  the 
Investigator  makes  up  are  really  representative  of  the  structure  of  the  measurement  space.  Hope¬ 
fully,  the  categories  and  sensors  chosen  present  each  category  as  an  Isolated  cluster  In  the  measure¬ 
ment  space.  However,  this  does  not  always  happen:  sometimes  there  are  two  clusters  In  the 
measurement  space  which  the  Investigator  groups  together  as  the  same  category.  In  this  case, 
structural  Information  of  the  measurement  space  Is  lost,  since  once  categorization  Is  done  there  can 
be  no  further  distinctions  made  between  these  two  clusters.  Sometimes  what  the  Investigator  calls 
two  categories  actually  form  only  one  cluster  In  the  measurement  space.  In  the  first  case  the 
Investigator  Is  losing  Information,  and  In  the  second  case  he  Is  wasting  effort. 

On  the  basis  of  these  Ideas  we  will  restate  the  pattern-recognition  problem.  Given  an  em¬ 
pirical  probability  distribution  on  measurement  space,  we  would  like  to  find  a  rule  (equivalence  rela- 
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tion)  which  partitions  measurement  space  so  that  each  cell  of  the  partition  represents  a  ciuster  or 
similarity  set  in  measurement  space.  Each  point  in  such  a  set  is  similar  to  'very  other  point  in  tie 
set  because  they  all  belong  to  ^e  same  grouping  or  cluster.*  Pattern  recognition,  approached  from 
this  perspective,  becomes  a  method  of  handling  data  in  terms  of  families  of  aggregates,  each  aggre¬ 
gate  being  a  cluster  of  points . 

The  information  which  J  measurement  conveys  to  the  investigator  is  not  contained  in  the  pre¬ 
cise  values  of  the  components  of  the  N-tuple,  for  those  values,  if  perturbed  Just  a  little  would  still 
usually  give  a  measurement  conveying  as  much  information  to  the  investigator  as  in  the  un$)erturbed 
case.  Rather,  the  information  which  a  measurement  conveys  to  the  investigator  is  contained  in  the 
aggregate  family  to  which  it  belongs. 

We  illustrate  this  idea  about  information:  suppose  there  are  five  different  kinds  of  objects 
or  occurrences  in  the  environment  sensed.  Further  suppose  that  to  each  object  or  occurrence  there 
corresponds  an  ideal  measurement,  the  measurement  which  would  be  made  if  the  environment  were 
not  stochastic  and  if  the  sensor  system  were  accurately  calibrated,  completely  linear,  and  noise¬ 
less,  However,  due  to  all  the  conditions  mentioned,  measurements  of  the  same  objector  occurrence 
will  deviate  a  little  from  the  ideal,  distributing  themselves  around  the  ideal  as  illustrated  in  Figure  3. 
Since  these  deviations  were  caused  by  the  random  effects  of  the  environment,  thermal  noise  in  the 
sensors,  etc.,  the  investigator  is  not  interested  in  their  exact  values.  In  a  classification  problem 
the  only  item  of  interest  is:  which  object  or  occurrence  caused  the  measurement  we  made — i,e., 
which  aggregate  family  does  the  measurement  come  from? 

A  pattern-recognition  system  which  would  handle  the  data  in  terms  of  aggregate  families 
would  have  large  application  to  the  remote  sensing  field.  Without  any  knowledge  or  training  from 
the  investigator,  the  system  could  transmit  back  the  occurrence  of  an  aggregate  family  instead  of 
the  actual  measurement  made.  This  would  yield  a  large  saving  of  channel  capacity  and  preserve 
the  information  in  the  structure  of  the  environment  sensed. 

The  problem,  then,  is  how  to  go  about  designing  such  a  system.  The  key  to  this  lies  in  our 
understanding  of  similarity.  Similarity  of  measurements  means  that  the  difference  between  the 
measurements  is  statistically  unimportant.  The  reason  the  difference  is  unimportant  is  that  two 
measurements  which  are  similar  both  belong  to  the  same  cluster  (to  the  same  aggregate  family),  and 
what  characterij:es  the  location  of  the  cluster  also  characterizes  the  location  of  the  measurements 
within  it.  This  is  the  basis  of  similarity. 

Clusters  could  be  located  by  finding  subsets  of  measurement  space  whose  general  location 
also  characterizes  the  measurements  within  it.  Location  implies  a  given  metric,  however,  the  gen¬ 
eral  problem  may  have  a  measurement  space  which  can  have  no  metric  on  it  at  all.  The  generali¬ 
zation  of  location  may  be  easily  done  with  measures  of  concomitant  variation,  an  example  of  which 
follows.  Suppose  we  have  in  mind  a  subset  which  we  think  is  a  cluster.  Every  time  a  measurement 
occurs  which  is  in  the  subset,  we  will  designate  it  "a."  Every  time  a  measurement  occurs  which 
is  not  in  the  subset,  we  will  designate  it  "b."  Here  there  is  a  100%  concomitant  variation  of  the 
measurements  in  the  cluster  designated  by  the  letter  "a."  If  we  wanted  to  remember  these  associa¬ 
tions,  we  could  store  all  those  points  designated  by  the  letter  "a”  in  one  place  in  memory  and 
those  points  designated  by  the  letter  "b"  in  another  place  in  memory.  However,  if  the  measure¬ 
ment  space  has  high  dimensionality,  as  in  multi-sensor  systems,  it  would  have  so  many  points  with¬ 
in  it  that  storing  all  the  associations  would  be  almost  impossible  if  not  a  waste  of  machinery  and 
money . 

Thus  we  divide  each  measurement  into  a  number  of  characteristics,  just  as  we  could  divide  ^ 
the  perception  of  a  chair  into  a  number  of  characteristics  which  describe  it.  For  the  chair,  these 
might  be: 

characteristic  1-2  —  number  of  legs,  3  or  4 

characteristic  3  —  back  or  no  back 

characteristic  4  --  if  back,  is  It  hard  or  soft 

characteristic  5  —  soft  or  hard  seat 

characteristic  6  —  arms  or  armless 

characteristic  7-14  --  color  of  upholstery:  red,  orange,  yellow 

green,  blue,  indigo  or  violet 

characteristic  15-18  —  material  of  upholstery  (cotton,  leather, 

or  plastic) 

*We  will  use  the  concepts  cluster  and  similarity  set  interchangeably.  When  the  emphasis  is  on  geo¬ 
metric  compactness  we  will  use  the  concept  cluster.  When  the  emphasis  Is  on  the  similarity  of  the 
points  within  a  cluster,  we  will  use  the  concept  similarity  set. 
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Kow  wc  try  to  find  the  association  between  the  letter  "a"  and  each  one  of  the  characteristics. 
If  "a"  is  always  associated  with  the  occurrence  of  the  l^h  characteristic,  then  a  measure  of  concoml- 
tarit  variation  between  the  letter  "a"  and  the  characteristic  would  attain  its  maximum.  If  It  Is 

••  never  as socla ted  with  the  occurrence  of  the  characteristic,  then  a  measure  of  concomitant  varia¬ 

tion  between  the  letter  "a"  and  the  characteristic  v/ould  attain  its  minimum.  To  find  the  total 
measure  of  concomitant  variation  between  a  measurement  and  the  letter  "a,"  we  could  appropriately 
sum  the  individual  measures  between  each  of  the  characteristics  and  the  letter  "a."  If  we  wanted 
to  compare  the  general  characteristics  of  the  subset  (associated  with  the  letter  a  )  with  any  o*'® 
of  the  measurements  belonging  to  it,  we  could  use  the  measure  of  concomitant  variation  Mckwards . 
We  could  2sk;  given  that  we  have  a  measurement  which  has  been  assigned  the  letter  "a,  is  each 
characteristic  more  likely  to  be  associated  with  the  letter  "a"  or  is  It  more  likely  not  be  associated 
with  the  letter  "a"?  To  answer  the  question,  we  could  look  at  the  measure  of  concomitant  variation 
between  the  1*^  characteristic  and  the  letter  "a."  If  It  Is  larger  than  its  neutral  concomitant-variation 
value,  then  we  could  answer  by  saying  that  the  l^h  characteristic  is  Indeed  associated  to  some  extent 
with  the  letter  "a";  if  It  is  less  than  its  neutral  concomltant-varlatlnn  value,  then  we  could  answer 
by  saying  that  the  l^h  characteristic  is  not  associated  with  the  letter  "a."  Using  such  a  measure  of 
concomitant  variation,  a  system  could  find  a  cluster  by  finding  a  subset  of  measurement  space  for 
which  the  general  characteristics  of  the  subset  are  sufficiently  related  to  Its  measurements.  This 
relatedness  is  a  relation  of  prediction  or  estimation. 

A  subset  of  measurement  space  can  be  called  a  cluster  if  we  are  able  to  predict  or  estimate 
the  kind  of  measurements  which  are  representative  of  It  from  the  characteristics  we  have  chosen. 

The  kind  of  system  we  seek,  then,  is  one  which  produces  an  output  from  which  It  can  predict  ih® 
put  measurement  or  some  similar  measurement.  Such  a  system  is  an  adaptive  system  with  a  predictive 

criterion. 


An  adaptive  system  Is  one  which  examines  Its  own  output,  evaluating  It  according  to  some 
built-in  criterion,  and  then  modifies  its  structure  In  a  way  which  tends  to  Improve  the  evaluation  of 
the  output.  A  predictive  criterion  Is  one  which  attaches  high  value  to  the  ability  of  the  system  to 
predict  or  estimate  the  input  pattern  or  some  similar  pattern,  given  only  the  output. 

We  not  proceed  to  design  a  system  of  the  type  we  Just  described.  It  Is  a  two -la  ye  red* 
adaptive  system  with  a  predictive  criterion.  First  the  following  definitions  are  in  order;  let 

.  /j  ,  I  e  1 , . .  .K  be  the  set  of  possible  quantized  measurements  from  the  l^^'  sensor  in  a  total 

of  K  sensors .  Measurement  space  G  is  the  set  of  all  possible  measurements;  G  «  Lj  x  L2  x, . .  .x  Lj^. 

ror  any  subset  A  of  G.  ACG,  define  A^  '  {  •  Construct  characteristic  or  feature  sets 

Ej,  I  »  1 ,  a .  .K  such  that  the  probability  of  is  strictly  greater  than  zero  P  (E^)  >  0,  and  for  every 


gcG  where  P  (g)  >0,  there  exists  a  vector  6  (g)  =^62'^such  that  {g}  =  ^  Let  the  pattern  set 

Abe  the  set  of  all  such  possible  vectors.  A  pattern  here,  is  considered  to  be,  not  the  measurement 
vector  itself,  but  a  binary-coded  N-tuple  of  the  measurement  vector. 

Of  course,  measurements  are  not  naturally  binary  but  usually  analog.  At  some  point  in  the 
recognition  process,  if  only  during  transmission  the  analog  measurement  is  quantized  and  converted 
tob^ary  form.  The  binary  code  Is  determined  by  using  information-theory  techniques  on  channel 
capacUy  and  noise  distribution.  However,  the  binary  N-tuples  which  we  have  called  patterns 
have  little  to  do  with  such  transmission  codes.  Rather,  the  binary  codes  which  form  the  patterns 
have  to  do  with  feature  or  characteristic  sets.  We  Illustrate  this  In  Figure  5,  where  the  large 
rectangular  area  represents  measurement  space.  The  fine  partition  on 
the  quantizing.  Each  measurement  which  occurs  within  a  particular  cell  of  the 

the  same  Quantized  value.  The  outlined  area  represents  a  connected  subset  o.  measurement  space . 
Each  point  within  the  subset  is  close  to  other  points  within  the  subset.  Such  a  subset  cati  be  called 
a  feature  or  characteristic  set.  since  It  characterizes  In  this  particular  case  what 
call  the  upper-left-hand  corner  of  measurement  space.  If  a  measurement  Is  made  which  belongs  to 
the  chara  ceteris  tic  subset,  we  can  code  the  l‘^^  component  of  the  pattern  corresponding  to  the  measure¬ 
ment  as  a  binary  1;  otherwise,  as  a  binary  0. 


N 


•  Ihb  term  layer  Is  used  In  accordance  »lth  the  terminology  Rosenblatt  developed  for  his  perc^tlor. 
devices.  In  the  perception,  a  layer  refers  to  a  grouping  of  associative  units  which  function  together. 
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This  kind  of  coding  is  quite  different  from  a  regular  binary  code.  Suppose  that  the  measure¬ 
ment  vector  has  two  components  and  that  each  component  is  bounded  by  0  and  4  inclusive.  Measure¬ 
ment  space  quantized  by  unit  intervals  is  illustrated  in  Figure  5.  A  regular  binary  code  for  the  quan¬ 
tized  measurement  can  be  formed  as  follows: 


(0,0) 

00  00 

(2.0) 

10  00 

(0,1) 

00  01 

(2.1) 

10  01 

(0,2) 

00  10 

(2.2) 

10  10 

(0,3) 

00  11 

(2.3) 

10  11 

(1.0) 

01  00 

(3,0) 

11  00 

(1.1) 

01  01  - 

(3,1) 

11  01 

(1.2) 

10  01 

(3.2) 

11  10 

(1.3) 

01  11 

(3.3) 

11 11 

The  first  two  binary  bits  represent  the  binary  code  for  the  first  component  while  the  last  two  binary 
bits  represent  the  binary  c^e  for  the  last  two  components. 

■  Figures  6,  7,  8,  and  9  illustrate  the  characteristic  sets  which  would  correspond  to  the  first, 
second,  third  and  fourth  «nary  bits  in  the  code. 

In  general  such  codes  do  not  have  corresponding  characteristics  which  actually  characterize 
a  connected  neighborhood  in  measurement  space.  In  fact,  the  2nd  and  4th  bit  Just  indicates  whether 
the  first  and  second  components  arc  even  or  odd.  Whether  a  quantized  measurement  is  even  or  odd 
does  not  have  the  slightest  bearing  on  recognizing  categories  such  as  corn  or  wheat. 


For  any  S-dimensional  real  vector  ''2  1  function  Sgn:  -  {-1 ,  l}®  be  defined 

(V<  whore  «  vj  -  ;  <  0  '  ‘  s  X  'i-  '''  ’> 


by  Sgn  (V' if  Vj  -  V  <0  '  *  * *  S  A  j 

binary  {-1  ,l}  for  the  output  categories .  The  form  of  the  sUoicture  for  the  first  layer  of  the 
adaptive  system  is  defined  by  the  equation  tj  ■  Sgn  (Q^S)  where  Q*  (q^^)  is  the  matrix  of  concomi¬ 
tant  variation  coefficients  and  6  is  the  input  pattern  vector.  The  first  layer  is  illustrated  in  Figure 

10. 

Each  q..  is  a  measure  of  the  concomitant  variation  between  the  l^h  pattern  characteristic  with 

the  J^h  similarity  set  characteristic.  If  it  generally  happens  that  when  a  pattern  has  the  charac- 
teriitlc.  the  similarity  set  into  which  it  is  classified  has  the  jth  similarity  set  charagterlstic,  then 
a  will  be  positive.  If  it  generally  happens  that  when  a  pattern  does  not  have  the  i  pattern 
SWaracterisUc,  the  similar?ty  set  into  which  it  is  classified  does  not  have  the  similarity  set 
characteristic,  then  q^  will  also  be  positive.  In  the  other  two  cases  q^^  will  be  negative.  Thus 

N  . 

V  6  q  Is  an  overall  measure  of  concomitant  variation  between  the  input  pattern  6  and  the  J 
ij  1  iJ 

similarity-set  characteristic.  If  this  sum  is  great  enough,  the  system  will  decide  that  the  similarity 
set  to  which  it  will  classify  6  will  have  the  Jih  similarity  set  characteristic:  otherwise  it  will 

decide  that  it  does  not. 

The  binary  {-1  ,l}  vector  n  is  the  output  of  the  system,  n  fs  «  function  of  the  input  6  and 
the  present  Q  matrix,  which  is  conditioned  by  past  Inputs.  The  system  tests  its  output  rj  by  forming 
a  prediction  or  estimation  of  6.  This  is  done  in  the  second  layer  of  the  system  whose  st^ctural 
fo?m  is  described  by  the  equation  ^  *=  Sgn  (Qq),  and  is  Illustrated  in  Figure  11 .  The  sector  t  is  the 
system's  prediction  or  estimation  of  6,  the  input  pattern.  The  output  is 
and  %  component  by  component.  The  system  then  modifies 

In  a  way  v!hlch  tends  to  improve  the  prediction  or  estimation,  t.  An  adaptation  rule  of  the  sort 
q  (t  +  ^1)  *=  q  ,  W  +  c  6j  Dj  l^i  "  1  produce  Just  such  a  change.  Supi'ose  6^  =  then  there 

would  be  no  cLge  since  16^  - =  0.  Suppose  6^  =  1  and^j  =  -l.  an  error  situation.  In  order 
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for  this  to  happen,  ^  is 

K-1 

appropriate  change  in  each 
^ik  **  ’’k  “  ^  decrease  q^j^  if 

M 


order  for  this  to  happen, 


small;  therefore,  the  s/stem  increases  each  term  in  the  sum  by  an 

k  *  1 , . .  .M.  To  increase  each  term  q^j^  rjj^  the  system  must  increase 
*  -1,  Suppose  *  -1  and  6^  *  1,  also  an  error  situation.  In 

is  large;  therefore,  the  system  decreases  each  term  in  the  sum  by 


k«l 


an  appropriate  change  in  each  q^j^,  k  «  To  decrease  each  term  q^j^  t)j^  the  system  must  de- 

drease  q^j^  if  rjj^  *  1  and  increase  q^j^  if  *  ■!.  We  may  summarize  these  changes  in  the  table  of 
figure  12  and  compare  them  with  the  change  c  6^  rjj  |®i  “  ^  I  Produces. 


An  adaptive  pattern-recognition  system  of  the  kind  Just  described  was  simulated  with  the 
GE-62S  computer,  using  K-band  radar  imagery  with  polarizations  HH,  HV,  W  taken  in  July  1966 
over  agricultural  areas  in  western  Kansas.  Since  automatic  image  digitization  equipment  was  not 
available,  a  few  line  scans  across  each  field  were  manually  taken  at  random  with  a  mlcro-densito- 
neter.  The  average  of  these  scans  was  used  as  an  estimate  for  the  average  radar  return  for  each 
field.  A  total  of  253  such  fields  consisting  of  land-usage  categories  of  bare  ground,  wheat  stubble, 
wheat  stubble-weeds,  wheat  stubble-mulch,  weeds,  com,  alfalfa,  grain  sorghum,  pasture,  and 
sugar  beets  were  examined. 


Measurement  space  in  this  example  consisted  of  all  possible  three-tuples,  the  first  compo¬ 
nent  being  the  average  return  from  HH,  the  second  component  being  the  average  return  from  HV,  and 
the  third  component  being  the  average  return  from  W.  From  our  data  the  minimum,  M^,  and  the  range, 

R,,  (maximum  minus  the  minimum)  for  the  1^^  component,  1*1,  2,  3  was  determined.  The  charac- 

I  t 


teristic  sets  were  defined  as: 


1  •  1 ,  2;  J  ■  0, . . .,5;  k  «=  I  +  1,  3;  m  ■  0,  ,,.,5.  The  geometric  configuration  of  the  characteristic 
sets  are  displayed  in  figure  13. 


The  adaptive  pattern-recognition  system  worked  as  follows.  For  each  iteration  a  measurement 
from  the  set  of  253  measurements  was  chosen  at  random.  Given  the  definition  of  the  characteristic 
sets,  the  system  determined  the  6  vector  which  corresponded  to  the  randomly  chosen  measurement. 
Via  the  first  layer  of  the  adaptive  system,  tj.  a  3  x  1  binary  {-1,1/  vector,  was  computed  from 
Sgn  Q'6.  Via  the  second  layer  of  the  adaptive  system  the  prediction  of  6  given  rj*  was  computed 
from  Sgn  Qt).  The  system  then  compared  6  and  Vand  reinforced  the  elements  q^  of  Q  appropriately 

and  started  another  iteration. 


The  reinforcement  parameter  c  was  chosen  at  .002.  Values  of  .009  and  .02  were  also  tried 
but  they  made  little  difference  in  the  final  classification.  (There  was  a  random  shift  of  at  most  seven 
points  with  the  different  values  of  c.)  The  cycling  of  iteration  after  iteration  was  kept  up  until  the 
total  classification  for  all  the  data  points,  based  on  the  updated  Q  matrix,  did  not  change  for  more 
than  30  cycles.  It  took  189  cycles  to  reach  this  situation. 

Figure  14  illustrates  a  scattergram  of  the  data,  and  figure  15  Illustrates  a  scattergram  of  the 
data  coded  by  their  respective  land-usage  categories  as  determined  from  ground  truth.  The  axes 
of  the  scattergrams  are  the  first  two  normalized  eigenvectors  of  the  covariance  matrix  for  the  data. 
These  axes  are  usually  called  the  first  two  principal  axes.  A  frequency  chart,  as  determined  after 
189  cycles,  of  similarity  sets  versus  the  land-usage  categories  is  shown  in  figure  16.  Similarity 
sets  II  and  III  are  closely  related,  their  difference  mainly  being  that  III  has  all  the  corn  while  11 
has  most  of  the  wheat  stubble-weeds.  For  illustration  and  basic  interpretation  purposes  we  group 
together  similarity  set  II  and  III  as  in  the  chart  of  figure  17.  If  we  now  group  together  those  cate¬ 
gories  which  appear  to  occur  in  mostly  the  same  similarity  sets,  we  will  obtain  the  chart  in  figure 
18.  Scattergram  of  the  similarity  sets  (clusters)  which  the  adaptive  system  determined  is  shown  in 
figure  19. 


Bare  ground  had  formed  a  cluster  of  low  returns;  wheat  stubble,  wheat  stubble-weeds, 
wheat  stubble-mulch,  and  pasture  formed  a  cluster  of  medium-low  returns  which  meshed  with  each 
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Figure  11 .  Second  Layer  of  Adaptive  System 


T 


•  ••  * 

••  •  • 

-  •  • 

•  •  •  •  •  • 

•  ••  •  •  ,  • 

.•  •  *•••/  •  ••  • 

••  iV..  \  •••  •. .  • 

*  •  *  *  *•*  fl*  *  ^  •** 

•  ••  *  f  •*  ,  ,  - 


•  •  • 


CR£S-t33 


fIGUaEM.  $CAIII«G»AM  OE  DATA  MOJICTIO  OH  H«ST  TWO  raiNCIPAl  AXES 


•  "A 


A*  A 


>■  ®  J  o"_ 
o  *1  /5i 


.  1  AA  0  ."'VSl  A 

•  ‘A  «•  ^  *  A*  <1^  o 

.  r*  r«w  t*  f.  «__ii — » — 

-  --t:  *  *t  * 


A  A 


•  ■<«((  CKOU« 

•  AlFAlF* 

■  (««iN  somH^ 

OCOKN 

*•!(» 

AwHCAT  s'ueBic  « 
w»«Ar  sn,BAf  t 

«h(AT  STUCAC  (WUMI 


CHES-IM 


#-nncft  DMA  PROJECTED  OH  EI«$T  TV.O  PBINCIPAl  AXIS 
IICUXE  IS  SCATTERGHAM  Of  COOED  DATA  rnojititi/ 


4 


- - P/ 

^  *1 
y*  » 

I 
I 
I 


✓••  I  •  ••  /• 

•*  •  .  •  / 

>  f  •  •  •'  ,  • .  I  .  • .  / .  • 

/•  V'"  *•.">••  ♦.  /••  • 

/  ••  ••  •  •  •  ••  •  «  «  ,  /  •  •  , 

is  r ..  / 

"-re.  .  .*.  •  r  y 

\  •  ••  / 

s*  •  •  / 

N*  •  •  ^ 

\  • 


\ 


FIGURC  19.  CLUSTtnS  THE  ADAPTIVE  PATTERN  RECOGNITION  SYSTEM  FOUND 
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Abstract 


A  number  of  the  constraints  with  which  remote  sensing  must  contend 
in  crop  studies  are  outlined.  These  include  sensor,  weather,  and  congru- 
encing  constraints;  the  nature  of  the  answers  demanded  of  the  sensor  system; 
and  the  complex  spatial  and  temporal  variances  of  both  crops  and  cropping 
system?,  in  areas  of  sub-continental  dimensions.  Attention  is  then  focused 
on  several  methods  which  may  be  used  in  the  statistical  analysis  of  multi¬ 
dimensional  remote  sensing  data. 

Crop  discrimination  for  radar  K-band  imagery  is  investigated  by  two 
methods.  The  first  method  uses  spatial  conditional  probabilities  in  conjunc¬ 
tion  with  a  nearest-neighbor-probability  approach,  and  the  second  method 
uses  the  standard  statistical  techniques  of  cluster  analysis  and  principal 
axes  representation. 

Results  indicate  that  crop  type,  height,  and  percent  of  cover  signi¬ 
ficantly  affect  the  strength  of  the  radar  return  signal.  Discrimination  is 
better  when  the  crops  are  mature.  Sugar  beets,  corn  and  very  bare  ground 
are  easily  distinguishable,  while  sorghum,  alfalfa  and  young  wheat  are 
harder  to  distinguish.  Distinguishability  will  be  improved  if  the  imagery 
is  examined  in  time  sequence  so  that  changes  between  times  of  planting, 
maturation  and  harvest  provide  additional  discriminant  tools. 
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Introduction 


Agriculturalists,  geographers  and  others  are  constrained  in  the  study 
of  crops  with  remote  sensors  by  the  kinds  of  agricultural  characteristics 
that  the  sensors  can  in  fact  detect.  The  sensors  used  in  the  various  por¬ 
tions  of  the  electromagnetic  spectrum  are  sensitive  to  dissimilar  energy- 
matter  interactions .  They  detect  different  characteristics  and  therefore 
convey  different  kinds  of  information.  In  the  visible  region,  molecular 
absorptions  produce  the  color  effects  which  convey  information  about  crop 
condition;  in  the  infrared  region,  the  cycle  of  thermal  response  under  an 
isolation  load  may  give  information  on  moisture  stresses  within  crops;  in 
the  radar  region  the  backscattered  return  is  primarily  related  to  surface 
roughness  and  dielectric,  which  may  be  related  to  crop  type,  percent  of 
cover,  moisture  content,  and  similar  variables. 

While  at  first  sight  radar  cannot  sense  "color"  for  its  information 
on  crop  condition,  it  may  detect  related  changes  in  crop  geometry  and 
moisture  and  thus  convey  in.formation  about  crop  type  and  state  which  are 
not  immediately  obvious .  Agriculturalists  may  wish  radar  and  other  remote 
sensors  to: 

(1)  detect  the  presence  of  different  crops--is  the  region  entirely 
homogeneous  or  are  there  di^erent  crops , 

(2)  determine  the  crops  actually  present, 

(3)  determine  the  boundaries  and  acreages  of  different  crops, 

(4)  determine  the  vigor  (state)  of  crops , 

(5)  determine  the  agent  responsible  for  any  loss  of  crop  vigor, 

(6)  predict  yields . 

This  paper  is  concerned  with  the  second  of  these  items:  determining 
the  crops  actually  present.  Studies  on  crop  discrimination  are  now  underway 
at  Purdue  University,  The  University  of  Kansas,  the  U.S.  Department  of 
Agriculture  Experiment  Station  at  Weslaco,  Texas,  and  elsewhere,  with  a 
view  to  the  ultimate  use  of  remote  sensor  systems  in  automated  or  partially- 
automated  data  collection  and  analysis  of  crop  type,  condition  and  yields . 
Such  data  would  be  valuable  in  the  United  States  and  especially  so  in  under¬ 
developed  regions  where  adequate  crop-yield  and  similar  information  is  not 
available. 


Because  radar  exhibits  the  combined  advantages  of  fine  resolution^ 
virtual  all-weather  capability,  and  Independence  of  the  sun  for  Illumination, 
and  moreover  is  without  serious  degradation  by  the  atmosphere,  its  abilities 
and  shortcomings  as  a  sensor  of  crops,  pastures,  and  natural  vegetation 
are  matters  of  no  small  moment  to  the  agricultural  community.  A  great  deal 
of  research  will  be  needed  to  clearly  define  these  abilities  and  shortcomings, 
since  little  work  has  been  carried  out  on  radar  as  a  sensor  for  agricultural 
purposes.  The  study  reported  here  is  one  of  a  handful  which  explores  these 
topics.  All  are  concerned  with  a  single  frequency  and  one  (less  frequently 
two  or  four)  polarization (s)  only,  and  thus  fall  well  short  of  an  ultimate 
capability  of  multi-frequency,  poly-polarizatlon  radars  that  are  feasible 
for  agricultural  surveillance.  The  underlying  concept  behind  synchronous 
data  acquisition  with  multi-frequency  polarization  radars  is  that  it  repre¬ 
sents  the  microwave  equivalent  of  multiband  spectral  reconnaissance  in 
the  photographic  and  thermal  IR  region,  which,  as  Hoffer  (1967)*has  shown, 
has  produced  some  encouraging  results  in  tests  at  Purdue  University. 


Temporal  and  Spatial  Components  of  Crop  Discrimination 

Ideally,  in  the  separation  of  crop  types  with  a  single  sensor,  we 
hope  that  there  will  be  a  narrow  spread  in  the  probability  distribution  of 
the  radar  reflectances  for  each  crop:  we  also  hope  that  the  probability 
distributions  for  the  radar  reflectances  of  each  crop  will  not  overlap.  This 
is  of  course  not  the  case  with  radar  or  any  sensor,  the  reason  being  that  in 
that  discrimination  between  crops ,  the  attribute  we  are  seeking — genetic 
difference  between  crops  must  be  approached  through  the  surrogate  of 
crop  reflectance.  The  surrogate  may  be  quite  Imperfectly  correlated,  in  a 
discriminatory  sense,  with  the  defining  genetic  attributes.  Furthermore, 
there  is  considerable  natural  variability  in  the  time  of  planting,  varieties 
planted,  and  conditions  of  crops  as  a  function  of  time  of  planting,  soil 
type  and  so  on.  Consequently  no  single  sensor  is  able  to  effect  clear-cut 
discrimination  at  a  single  time  of  ell  crops  in  all  fields.  It  becomes 

* 

Hoffer,  R.  M.,  Interpretation  of  Remote  Multispectral  Imagery  of  Agricul¬ 
tural  Crops.  Research  Bulletin  831,  Laboratory  for  Agricultural  Remote 
Sensing,  vol.  1,  Purdue  University  Agricultural  Experiment  Station. 
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necessary,  therefore,  to  increase  the  probability  of  correct  discrimination 
of  crops  through  the  use  of  several  frequencies  and  polarizations  in  con¬ 
cert,  both  synchronically  and  diachronically. 

It  is  only  through  empirical  testing  of  numerous  sites  in  the  United 
States  (and  elsewhere)  at  repeated  intervals  throughout  a  number  of  growing 
seasons  (probably  at  least  five)  that  we  will  be  able  to  sufficiently  define 
the  time-sequential  probability  distributions  of  the  radar  returns  for  different 
crops.  It  is  immediately  apparent  that  if  crops  are  spatially  and  temporally 
distributed,  then  the  distribution  of  the  reflectances  within  a  single  growing 
season  and  between  growing  seasons  will  demand  a  number  of  information 
channels  to  insure  identification,  which  may  become  unmanageably  large 
even  if  only  a  few  looks  through  the  growing  seasons  are  available.  This 
Immediately  creates  a  compromise  in  the  amount  of  Information  obtained 
at  any  time  and  the  number  of  times  it  is  obtained.  The  various  trade-offs 
between  number  and  kind  of  sensors  and  number  and  timing  of  data  collec¬ 
tion  by  the  sensors  has  yet  to  be  worked  on,  but  will  be  an  essential  ingre¬ 
dient  of  continuing  studies  in  the  U.S. 

From  this  brief  discussion  it  might  intuitively  be  felt  that  multiple 
looks  through  time  may  turn  out  to  be  very  important  in  remote  sensing,  and 
that  through  multiple  looks  it  may  be  feasible  to  reduce  the  number  of  infor¬ 
mation  channels  needed  at  each  time.  Consider,  for  example,  the  condition 
of  crops  and  their  times  of  planting  and  harvesting  in  the  United  States. 
There  are  many  differences  in  kinds  of  crops,  growing  season  and  growth 
patterns  between  west,  center  and  east;  between  humid,  subhumld,  and 
Irrigated;  and  between  south,  center  and  north,  to  name  only  the  more 
obvious  factors  involved  in  crop  geography.  This  means  that  from  March 
through  October,  every  week  is  a  remote  sensing  week,  since  critical 
changes  in  the  status  of  crops  are  occurring  somewhere. 

There  are  various  degrees  of  crop  Identification  accuracy  which  we 
may  want  from  remote  sensing  systems.  For  example,  we  may  wish  (1)  to 
positively  identify  the  crops  in  every  field  in  an  area  and  to  determine  the 
acreages  of  each  field  with  an  error  of  1%  or  less;  or,  less  difficult,  (2)  to 
identify,  with  some  acceptable  error,  the  kind  of  crops  in  each  field  of  a 
region,  and  to  estimate  total  acreages  of  each  crop  in  the  region  accurate 
to  +5%. 
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If  we  want  to  determine  the  exact  acreages  of  crop  being  grown  in 
an  area,  either  for  statistical  reporting  purposes  or  for  checking  on  wheat 
and  feed  grain  program  compliance,  we  have  a  situation  described  by  the 
first  case  above.  In  such  a  circumstance,  if  ironclad  identification  cannot 
be  made  in  a  single  pass,  then  multiple  passes  must  be  made,  with  the 
attendant  problems  of  congruencing:  storing  spatial-reference  coordin<ites 
for  all  points  on  the  ground  and  of  rendering  time-separated,  different- 
flight-path  data  in  congruent  planimetric  geometry.  This  is  no  mean  task! 

At  present  it  is  really  an  intolerable  burden  to  place  on  a  remote  sensor 
system. 

Alternatively,  as  in  the  second  case  above,  we  may  be  able  to  avoid 
the  necessity  for  many  passes  by  allowing  modest  errors  in’ identification 
and  by  accepting  regional  average  crop  acreages  with  up  to  a  5%  error. 

If  we  select  this  alternative,  however,  a  more  stringent  requirement  for 
spatiai  accuracy  must  be  met. 

As  we  have  already  mentioned,  one  of  the  major  problems  in  obtain¬ 
ing  high  identification  accuracy  with  any  automated  multi- sensor  system 
is  that  of  bringing  the  images  into  mutually  congruent  geometry.  Each 
resolution  cell  of  each  image  corresponds  to  a  small  area  on  the  ground. 

The  system  must  be  able  to  determine  the  correspondence,  in  planimetric 
coordinates,  between  each  resolution  cell  of  each  image  and  each  resolu¬ 
tion  cell  of  every  other  image,  and  between  each  resolution  cell  and  each 
small  area  on  the  ground.  This  correspondence  must  be  established  auto¬ 
matically  if  the  system. is  to  work  quickly  and  effectively. 

The  problem  of  image  congruence  has  been  attacked  by  researchers 
investigating  the  visible  and  infrared  region  of  the  spectrum.  Multiple 
wavelengths  and  mechanical  rotating  mirror-scanners  have  been  developed 
in  which  the  signal  received  is  split  into  different  wavelength  bands ,  each 
with  its  own  detector.  The  signals  are  then  recorded  simultaneously  on 
different  channels  of  a  magnetic  tape.  A  similar  technique  is  feasible  with 
radar  systems  operating  with  multiple  polarizations  and  wavelengths , 
although  direct  recording  of  the  video  calls  for  a  fairly  wide-bandwidth 
recorder.  The  use  of  some  form  of  intermediate  storage  so  that  only  the 
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actual  significant  bandwidth  of  the  signal  need  be  recorded  (sometimes 
hundreds  of  pulses  are  averaged  to  make  one  resolution  cell)  can  reduce 
this  problem. 

A  more  severe  problem  of  image  congruence  comes  from  multiple 
flights .  If  an  aircraft  or  spacecraft  carrying  the  radar  could  fly  exactly 
the  same  path  within  a  wavelength  or  so  each  time,  the  problem  would  not 
be  difficult.  Unfortunately,  aircraft  cannot  repeat  either  horizontal  or 
vertical  location  of  the  path  to  anything  like  wavelength  precision;  further¬ 
more,  operational  conditions  may  require  significant  differences  in  altitude 
which  result  in  distortions  between  images .  Even  the  best-programmed 
flight  lines  sometimes  deviate  so  that  there  are  horizontal  differences  as 
well  as  vertical  differences  . 

The  need  for  automatic  congruencing  equipment  is  clearly  urgent. 
Multiple  recording  during  a  single  pass  solves  only  the  problem  for  that 
pass.  No  system  known  produces  images  on  successive  flights  that  are 
initially  congruent.  Digital  correlation  and  subsequent  linear  or  non-linear 
correction  of  scale  could  in  theory  solve  the  problem.  The  mass  of  infor¬ 
mation  collected  in  any  extensive  aircraft  flight  program,  however,  prevents 
the  use  of  these  digital  techniques  for  congruencing  in  the  agricultural 
problem.  For  a  geological  problem,  such  techniques  might  be  realistic 
because  the  flights  need  not  be  repeated  so  often,  and  a  delay  of  a  year 
or  two  in  processing  the  data  does  not  necessarily  result  in  loss  of  its 
value.  Some  faster  means  must  be  found  fcr  the  more  timely  problems  in 
agriculture  before  the  automated  systems  can  be  truly  useful. 

There  is  another  approach  to  the  image  congruencing  problem  which 
relies  on  the  fact  that  there  are  geographic  boundaries  which  are  relatively 
permanent.  Homogeneous  regions  (determined  by  "natural"  statistical 
classes)  may  be  found  on  each  image,  independent  of  the  other  images. 
Boundaries  must  exist  between  homogeneous  regions.  Congruencing  may 
be  achieved  by  distorting  the  geometry  of  the  images  so  that  there  is  max¬ 
imum  correlation  between  the  boundaries  on  each  image  with  those  on 
every  other  image.  The  distorted  ("congruenced")  images  may  now  be  used 
for  classification  purposes.  The  end  product  of  the  classification  may  be 
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FIGURE  1.  RADAR  IMAGERY 
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brought  to  planimetric  correspondence  with  the  ground  by  examining  the 
boundaries  around  the  categories  and  distorting  the  geometry  of  the  cate¬ 
gorized  image  so  that  it  matches  the  natural  boundaries  on  a  map. 

If  we  did  not  need  to  know  exactly  where  each  crop  occurred  in  a 
region,  a  relatively  crude  job  of  congruencing  could  be  performed.  Such 
unsophisticated  congruencing  would  retain  just  enough  statistical  structure 
to  tell  us  the  acreage  of  each  crop,  independent  of  its  spati.,J  position,  or 
the  percentage  of  occurrence  for  each  crop  in  the  region.  Little,  if  any, 
work  has  been  done  in  this  area,  and  we  believe  it  warrants  further  thinking. 

In  the  following  sections  of  the  paper  we  outline  a  number  of  ways 
of  data  manipulation  of  radar  reflectances  from  crops  with  v/hich  we  are 
experimenting  at  the  University  of  Kansas.  Before  doing  so,  however,  a 
few  words  are  necessary  about  the  area  where  tl  e  crop  data  has  been 
obtained . 

The  region  where  the  crop  information  is  collected  lies  a  few  miles 
west  of  Garden  City,  in' western  Kansas  (Figure  1).  In  this  area  there  are 
over  400  fields  for  which  data  were  collected  at  the  time  of  radar  over¬ 
flights.  The  area  was  chosen  because  fields  are  generally  large,  slopes 
are  extremely  gentle,  soils  are  uniform,  both  dryland  and  irrigated  cropping 
is  practiced,  and  there  is  some  diversity  of  crop  types,  agricultural  tech¬ 
niques  ,  and  moisture  conditions .  In  this  test  area  we  can  hold  constant 
in  an  analysis  such  factors  as  topography,  soil,  and  field  size,  and  can 
obtain  the  maximum  variations  in  soil  and  crop  states,  moisture  levels, 
and  so  on,  in  order  to  quantitatively  explore  both  crop  discrimination  and 
the  contribution  of  different  parameters  to  the  radar  backscatter. 


Data  Analysis 

Sensor  systems  are  designed  to  obtain  infomation  about  an  inter¬ 
esting  statistical  structure  of  an  environment,  where  "interesting"  is 
defined  by  the  investigator.  The  sensor  system  must  reflect  or  transmit 
this  interesting  environmental  structure  onto  the  data  structure  if  anything 
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Is  to  be  gained  using  the  sensor  system.  Thus,  once  a  sensor  system  has 
been  conceived  and  designed,  the  question  that  must  be  answered  when 
It  Is  tested  is  "what  statistical  structure  of  the  environment  is  actually 
preserved  through  the  sensor's  'eyes' — in  particular,  is  the  interesting 
structure  preserved?" 

In  order  to  answer  this  question  some  way  of  examining  the  data, 
some  way  of  determining  the  structure  within  the  data,  must  be  devised. 
Then  the  structure  of  the  data  may  be  compared  with  the  known  structure 
of  the  environment,  and  the  question  can  be  answered. 

The  most  usual  type  of  environmental  structure  which  is  defined 
as  "Interesting"  is  that  of  a  class ificatory  nature.  The  Investigator, 
arranges  the  environment  into  a  set  of  mutually  excL.sive  categories .  It 
is  his  hope  that  the  sensor  system  will  preserve  enough  of  the  environment's 
structure  so  that  data  obtained  from  the  sensor  system  may  also  be  cate¬ 
gorized  into  a  set  of  mutually  exclusive  classes,  where  each  class  corres¬ 
ponds  to  a  unique  category  in  the  environment. 

In  order  to  determine  the  strur*-are  of  the  data,  the  following  oper¬ 
ations  are  valuable: 

(1)  the  data  set  must  be  represented  in  some  geometric  way  such 

as  a  scattergram,  where  the  data  structure  is  immediately  appar¬ 
ent.  This  representation  must  be  perceivable  directly  as  a 
"picture." 

(2)  a  statistically  optimum  classification  of  the  data  must  be  made, 
optimality  being  defined  by  the  investigator. 

(3)  the  natural  classes  of  the  data  structure  must  be  found  and 
compared  to  the  categories  in  the  environment. 

The  information  received  in  a  "picture"  of  the  data  structure  gives 
the  investigator  an  immediate  "feel"  for  what  is  actually  happening.  It 
gives  him  an  intuitive  base  which  he  can  expect  to  guide  him  to  more 
sophisticated  techniques.  Further,  it  may  quickly  tell  the  investigator 
whether  there  have  been  faults  in  the  sensor  system  while  the  system  was 
gathering  the  present  set  of  data.  If  so,  then  the  data  would  not  have 
to  be  further  processed  and  the  experiment  could  be  rerun.  The  information 
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gained  by  the  statisticaliy  optimum  classification  of  a  data  set  is,  of 
course,  most  fundamental  in  the  sensor  system  test.  From  it  the  inves¬ 
tigator  finds  out  how  good  the  system  really  is  In  preserving  the  classifi- 
catory  structure  of  the  environment.  Finally,  by  finding  out  the  natural 
classes  in  the  data  structure,  the  investigator  can  learn  "what  language 
the  sensor  system  speaks."  If  a  natural  class  in  the  data  structure 
actually  corresponds  to  two  environmental  categories ,  then  we  would 
expect  that  regardless  of  the  optimality  of  any  statistical  classifying 
rule,  these  two  environmental  categories  are  going  to  be  confused,  and 
the  Investigator  can  do  nothing  about  it  except  search  for  a  new  sensor 
system.  This  confounding  of  categories  must  then  be  taken  into  account 
in  any  Interpretation  of  data  from  the  sensor  system. 

In  the  case  where  the  sensors  take  data  in  a  sequential  way  and 
it  is  known  that  observations  which  are  close  to  each  other  in  the  sequence 
are  likely  to  be  in  the  same  category,  some  effort  should  be  made  to  pre¬ 
serve  the  similarity  of  neighboring  data  points  in  their  classification.  For 
imaging  sensors  in  particular,  this  means  that  two  neighboring  resolution 
cells  should  not  be  classified  independently.  An  exception  to  this  occurs 
when  the  image  data  are  like  that  of  Figure  2,  which  has  the  kind  of  struc¬ 
ture  illustrated  in  the  scattergram  in  Figure  3.  Here  there  is  complete 
isolation  betv/een  the  categories,  with  bare  ground  occupying  the  lower 
third  of  densities,  grain  sorghum  occupying  the  middle  third  of  densities, 
and  corn  occupying  the  upper  third  of  densities . 

Representing  the  data  structure  in  a  two-  or  three-dimensional 
scattergram  is  a  problem  when  the  measurements  the  sensor  system  takes 
are  n-dlmenslonal  (n  being  more  than  three) .  Some  criterion  for  choosing 
the  representation  must  be  determined.  Suppose  we  decide  that  we  would 
like  to  obtain  the  best  representation  possible  by  translating  the  data  and 
then  "squeezing  it  through  a  two-dimensional  hole"  V7ith  a  linear  trans¬ 
form.  By  "best"  we  will  mean  that  the  mean  squared  error  of  a  least 
squares  estimate  of  the  original  set  of  data,  obtained  from  the  data  which 
have  been  squeezed  through  a  two-dimensional  hole,  is  minimized. 
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FIGURE  2.  Simple  image,  where  the  density  for  each  resolution 
cell  is  indicated  by  the  number  appearing  in  the  top  of  each 

box  and  the  category  of  surface  cover  which  corresponds  to 
the  cell  is  indicated  by  the  two  letter  code  appearing  at  the 
bottom  of  each  box. 


bare  ground  grain  sorghum  corn 


density 


FIGURE  3  .  Scatterograms  of  simple  image  data:  The 
categories  are  all  isolated. 
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It  can  be  shown  that  this  representation  corresponds  to  finding  a 
plane  in  n-dimensional  space  where  the  data  projected  onto  this  plane 
have  maximuin  variance.  Mathematically  the  situation  can  be  represented 
as  finding  the  projection  operator  P  which  has  rank  2  and  a  translation 
vector  t  such  that  the  mean  squared  error  ^ ,  defined  by  =  E  ^^(x-t)  - 

P(x-t)^'  ^(x-t)  -  P(x-t))  j  ,  is  minimized.  The  minimization  is  achieved  when 

^  9nd  P  projects  onto  the  plane  spanned  by  the  two  eigenvectors, 

Vf,  Va^  with  largest  eigenvalues  of  e[(x-E[x])  (x-e[x])']  ,  the  covariance 
matrix,  ‘(see  Kendall  and  Stuart  1964).  We  will  call  this  representation 
the  principal  axes  representation.  The  coordinates  of  any  data  vector 
relative  to  the  principal  axes'  plane  can  be  found  easily  as  (x-E[x])'  v^ 
for  the  first  component  and  (x-E[x])'  for  the  second  component. 

Although  there  are  many  ways  of  classifying  data,  such  as  linear 
decision  theory,  which  constructs  hyperplanes  boundaries,  or  nearest- 
neighbor  search  procedures,  these  are  statistically  admissable  only  under 
severly  restrictive  conditions.  A  Bayes*  decision  rule  is  almost  always 
statistically  admissable  and  optimal,  in  addition  to  allowing  the  use  of  a 
loss  function  which  can  weigh  a  correct  classification  more  heavily  on  the 
more  Important  categories.  Therefore,  to  get  information  regarding  the 
bes^  possible  classification  which  can  be  made,  a  Bayes'  decision  rule 
may  be  employed. 

We  used  a  Bayes'  decision  rule  (see  Middleton  1960)  classifying 
each  resolution  cell  of  the  radar  Imagery  data  independent  of  its  neighboring 
cell.  This  was  done  because  we  did  not  have  a  Bayes'  decision  program 
taking  into  account  the  similarity  of  neighboring  resolution  cells.  In  order 
for  the  Bayes'  decision  rule  to  indicate  the  classificatory  structure  of  the 
data  rather  than  to  merely  "memorize"  the  data,  the  data  set  must  be  ran¬ 
domly  divided  in  half;  the  first  half  must  be  used  to  generate  the  statistics 
for  the  Bayes'  rule  and  the  second  half  must  be  classified  on  the  basis  of 
this  rule.  A  "goodness  criterion"  may  be  adopted,  such  as  the  total  per¬ 
centage  of  correct  classification  or  the  average  percentage  of  overall  cate¬ 
gories  of  correct  classification,  where  the  percentages  are  properly  weighed 
by  the  loss  function  adopted. 


Finding  the  natural  classes  in  the  data  structure,  as  seen  through 
the  sensor's  eyes,  is  perhaps  the  hardest  of  the  three  operations  suggested, 
since  the  problem  is  not  well-defined.  Certainly  some  information  about 
the  natural  classes  will  be  obtained  from  the  principal  axes  representation 
and  the  Bayes'  classification.  Categories  which  are  confused  will  definitely 
show  up.  However,  it  may  be  that  parts  of  the  environment  which  are  cate¬ 
gorized  by  the  investigator  as  the  same  will  show  up  in  different  natural 
classes,  indicating  that  further  distinctions  may  be  made. 

We  will  approach  the  natural-class  problem  of  image  data  by  two 
clustering  methods .  One  sequentially  forms  clusters  from  data  measure¬ 
ments,  called  centers,  which  have  the  highest  probability  in  any  small 
spatial  region  in  the  image  (a  spatial  conditional-probability  approach) ,  and 
includes  in  each  cluster  those  points  which  are  sufficiently  close  to  the 
center  in  measurement  space  and  whose  probabilities  do  not  differ  too  much 
from  the  probability  of  the  center  point.  The  other  forms  clusters  taxono- 
mically  on  the  basis  of  Euclidean  distance  in  measurement  space,  using 
the  multiple  linkage  model. 

Clustering  by  Spatial  Conditional  Probability 

The  conditional  probability  approach  first  quantizes  the  data  in  the 
following  way.  Let  the  density  in  the  ij^^  resolution  cell  be  = 

(Xj  ,X2^  . .  .xJI )  where  xj^^  is  the  output  from  the  k^^  sensor  when  the  k^^ 

sensor  is  observing  the  area  represented  by  the  ij^^  resolution  cell.  Let 

th 

R.  be  the  range  for  the  k  coordinate.  The  quantizing  function  Q  is 

...  1,  _  1  1,  R, 

defined  as  Q  (X)  =  (i, ,  L  , . . . ,  1  )  where  -"Vn  '  R.  <  x.  <  ^  ^  . 

1  ^  n  iO  j  ""  j  10 

In  forming  clusters,  priority  is  given  to  the  quantized  measurement 
vectors  for  which  the  conditional  probability  of  occurrence  in  an  image 
spatial  region  is  higher  than  the  joint  probability  of  its  occurrence  over  the 
whole  image.  For  example,  a  multi-spectral  image  may  be  partitioned  into 
a  number  of  connected  spatial  regions,  as  inFigure  4.  The  conditional 
probability  distribution  for  each  region  is  computed  and  compered  with  the 
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Figure  .4.  Multispectral  image 
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Joint  probability  distribution  for  the  entire  image.  All  the  measurement 
vectors  which  have  the  highest  conditional  probability  in  the  multi-spectral 
images  are  given  priority  in  forming  clusters.  Each  cluster  is  formed  around 
a  center--that  is,  a  quantized  measurement  vector  which  is  not  yet  in  a 
cluster  and  which  has  the  highest  conditional  probability  for  some  spatial 
region.  At  any  iteration,  for  a  point  X  to  be  included  in  the  cluster  just 
being  formed,  there  must  exist  a  point  Y  already  in  the  cluster  such  that: 


(1)  |x^  -  y^(  <  1 

(2)  P0()<€^P(C) 

(3)  P0C)<e2P(Z) 

(4)  POC)  <  ~  P(W) 

^2 


for  each  i 

where  C  is  the  center  and  ^  ^  is  a  prechosen 
parameter 

where  Z  has  the  highest  probability  of  any 
point  already  in  the  cluster  and  ^2  ®  pre- 

chosen  parameter. 

where  W  has  the  lowest  probability  of  any 
point  already  in  the  cluster. 


If  there  are  any  points  which  do  not  meet  those  criteria  they  are  put  into 
the  cluster  which  is  nearest  to  them. 

Clustering  by  Multiple-Linkage  Technique^ 


Multiple-linkage  cluster  analysis  has  been  widely  used  as  a  numer¬ 
ical  taxonomic  technique  in  the  biological  sciences  (see  Sokal  and  Sneath, 
1963).  However,  many  of  the  assumptions  made  in  classifying  biota  cannot 
be  made  in  analyzing  radar  images ,  because  there  is  no  evolutionary  rela¬ 
tionship  between  the  objects  being  categorized.  Nevertheless,  since  the 
measurement  vectors  which  can  be  clustered  together  are  in  some  sense 
similar,  the  technique  is  worthwhile  in  finding  the  natural  classes  of  a 

data  set. 

The  clustering  procedure  takes  the  following  form: 

(1)  a  matrix  of  the  similarity  coefficients  between  each  pair  of 
measurement  vectors  Is  constructed.  (For  the  multi-spectral  Images,  these 

are  the  measurement  vectors  In  each  resolution  cell.) 

(2)  the  matrix  is  scanned  and  all  pairs  of  measurement  vectors 

which  are  mutually  close  are  linked  together. 

(3)  a  new  matrix  is  computed  on  the  basis  of  the  paired  vectors  and 

single  vectors.  This  new  matrix  is  scanned  and  all  single  vectors  or  paired 
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vectors  which  are  mutually  close  are  linked  together.  This  is  continued 
iteratively  until  all  clusters  are  linked  together. 

Figure  5  illustrates  the  phenogram  obtained  from  the  multiple-linkage 
clustering. 

The  investigator  may  determine  the  natural  class  structure  of  the 
data  by  truncating  the  dendrogram  at  some  distance,  as  illustrated  by  the 
vertical  dashed  line  in  Figure  5.  Here  a  distance  coefficient  of  .4  gives 
a  reasonable  picture  and  the  natural  classes  are  (b) ,  (c) ,  and  (d) ,  as  indi¬ 
cated  by  the  left  margin.  The  dendrogram  may  also  be  used  as  a  guide  for 
further  analysis.  For  example,  the  sugar  beet  field  (202  BE)  which  linked 
with  the  corn  field  should  be  investigated  in  greater  detail  to  determine 
why  it  is  so  different  from  the  rest  of  the  sugar  beet  fields. 

There  are  some  important  differences  between  the  conditional- 
probability  clustering  technique  and  the  multiple-linkage  cluster  analysis . 

The  multiple-linkage  cluster  technique  analyzes  the  data  microscopically, 
linking  each  observed  measurement  vector  with  other  similar  observed 
measurement  vectors  and  simultaneously  forming  a  set  of  clusters.  In  order 
to  do  this,  an  N  X  N  matrix  must  be  stored  if  there  are  N  observed  measure¬ 
ment  vectors .  The  conditional  probability  clustering ,  a  macroscopic  tech¬ 
nique,  starts  with  the  most  important  cluster,  forms  the  clusters  sequentially, 
and  needs  a  storage  of  2N  places.  These  differences  reflect  themselves 
in  the  amount  of  data  the  techniques  can  handle.  The  conditional- probability 
clustering  technique  is  quick  and  easily  used  with  data  sets  containing 
thousands  of  measurement  vectors,  while  the  multiple- linkage  clustering 
is  slower,  and  data  sets  containing  even  hundreds  of  measurement  vectors 
are  at  the  limit  of  the  procedure.  However,  the  multiple-linkage  clustering 
compensates  for  these  disadvantages  by  giving  a  more  detailed  analysis  of 
the  data  set. 


E17 


FIGURE  5.  Section  of  a  dendrogram  showing  graphically  the  results  of 
cluster  analysis.  This  is  a  section  of  tne  actual  dendrogram  used  in  categorizing 
the  August  and  September  data.  The  section  of  21  data  points  was  selected  from 
the  original  231  because  it  shows  portions  of  three  of  the  four  major  categories 
clustered.  The  dashed  line  shows  the  break  point. 
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Statistical  Crop  Discrimination 

We  will  illustrate  this  threefold  approach  to  data  analysis  with  two 
sets  of  data  obtained  from  K-band  radar  Imagery  taken  over  agricultural  areas 
In  western  Kansas.  At  the  time  the  analysis  was  done,  image  digitization 
equipment  was  not  available  at  the  University  of  Kansas,  and  we  had  to 
satisfy  ourselves  with  only  the  average  radar  return  in  each  farming  field. 

This  field  average  was  obtained  by  taking  a  few  random  scans  over  each 
field  in  the  images  with  a  microdensitometer  and  using  the  average  of  the 
scan  lines  as  an  estimate  of  the  average  radar  return  for  the  field. 

The  first  set  of  images  was  obtained  in  a  time-sequential  fashion 
during  the  months  of  August  and  September  in  1965.  The  August  image  had 
only  one  polarization,  being  horizontally  transmitted  and  received,  while 
the  September  images  had  two  polarizations,  being  horizontally  transmitted 
and  both  horizontally  and  vertically  received.  All  three  images  were  taken 
over  the  same  test  site  near  Garden  City,  Kansas.  The  second  set  of  images, 
also  taken  over  this  area,  was  obtained  during  the  month  of  July  in  1966  and 
had  all  four  polarization  combinations .  We  wanted  to  find  out: 

(1)  the  agricultural  type  of  environmental  structures  the  radar 
imagery  system  did  preserve,  and 

(2)  whether  more  environmental  structure  could  be  preserved 

(a)  if  the  data  were  taken  in  time  sequence  with  few  polar¬ 
ization  combinations  or 

(b)  if  the  data  were  taken  at  the  same  time  with  all  the 
polarization  combinations. 

Figures  6  and  7  illustrate  scattergrams  with  the  principal  axes 
representation  for  these  data.  It  may  be  noticed  from  the  general  pattern 
that  as  the  agriculture  becomes  more  intense  (surface- cover  or  land-usage 
categories  going  from  bare  ground  to  sugar  beets) ,  the  radar  return  becomes 
greater.  This  correlates  most  directly  with  the  crop  percent  of  cover.  For 
both  sets  of  data,  we  immediately  see  the  overlapping  of  categories.  In 
the  July  data  the  bare  ground  and  wheat  stubble  overlap.  Then  the  other 
part  of  wheat  stubble  overlaps  with  alfalfa,  while  the  alfalfa  overlaps  with 


FIGURE  6.  Scattergram  of  August  -  September  1965  Radar  Imagery  field  averages  projected  onto  first 
two  principal  axes. 
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FIGURE  ?•  ScatterJram  of  July  1966  Radar  Imagery  field  .averages  projected  onto  first  two  principal  axes . 


the  corn,  and  grain  sorghum  overlaps  and  covers  both  alfalfa  and  com. 

Sugar  beets  stand  almost  isolated  by  themselves.  In  the  time-sequential 
data  bare  ground  stands  out  by  itself,  while  wheat  stubble  overlaps  with 
alfalfa  and  both  overlap  with  the  grain  sorghum.  The  corn  and  sugar  beets 
stand  by  themselves.  This  indicates  that,  for  the  July  data,  we  would 
expect  the  Bayes'  decision  rule  to  confuse  wheat  stubble,  alfalfa,  corn,  an 
and  grain  sorghum  with  each  other,  and  to  confuse  bare  ground  and  wheat 
stubble  with  each  other  while  classifying  sugar  beets  fairly  well.  For  the 
time-sequential  data  we  would  expect  the  Bayes'  decision  rule  to  confuse 
wheat  stubble,  alfalfa,  and  grain  sorghum  with  each  other,  while  classifying 
bare  ground,  corn,  and  sugar  beets  fairly  well.  This  is  indeed  true  if  we 
examine  the  tables  in  Figures  8  and  9,  indicating  the  results  of  the  Bayes' 
classification.  For  this  problem  the  loss  function  was  chosen  to  be  an 
equal  loss  function;  the  same  finite  negative  amount  was  lost  for  each 
correct  classification  and  nothing  was  lost  for  an  incorrect  classification. 
Notice  that  90%  of  the  data  grouped  as  shown  were  correctly  classified 
for  August-September,  while  only  '  8%  of  the  July  data  were  correctly  clas¬ 
sified.  We  would  expect  a  natural  class  approach  to  show  that  the  cate¬ 
gories  confused  in  the  Bayes'  classifying  technique  all  belonged  to  the 
same  natural  class  or  cluster. 

The  spatial  conditional- probability  approach  could  not  be  fully 
utilized  because  of  the  lack  of  image-digitization  equipment.  However, 
a  suitable  approximation  was  tried.  All  field  averages  which  were  measure¬ 
ments  of  the  same  land  usage  category  were  lined  up  together.  In  this  way 
the  spatial  closeness  of  similar  types  of  fields  in  the  data  sequence  could 
be  utilized  exactly  as  it  would  have  beer  if  the  data  sequence  represented 
resolution  cells  of  an  image.  Each  spatial  region  consisted  of  a  set  of 
10  sequential  vectors  in  the  sequence.  The  error  from  this  approximation 
results  mainly  from  the  fact  ^hat  clustering  procedures  of  this  type  require 
large  amounts  of  data  to  work  well,  and  using  field  averages  for  resolution 
cells  reduces  the  data  size  by  two  or  three  orders  of  magnitude.  Despite 
this  har.dicap,  the  conditional-probability  clustering  procedure  indicates 
that  the  categories  of  vegetation  which  were  confused  in  the  Bayes'  classi¬ 
fier  do  form  fairly  isolated  clusters  or  similarity  sets,  as  the  tables  in 
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were  completely  confused  with  one  another. 
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Figures  10  and  11  show.  Although  the  Bayes'  procedure  categorized  the 
August-September  data  better  than  the  July  data,  the  clusters  found  in  the 
July  data  corresponded  better  to  the  surface  cover  categories  than  those 
clusters  found  in  the  August-September  data.  This  was  due  to  the  fact  that 
there  was  no  physical  separation  between  the  confused  surface  cover  cate¬ 
gories  as  may  be  seen  from  the  scattergrams .  We  should  note  that  these 
clusters  were  formed  without  any  information  regarding  "ground  truth",  i.e. , 
knowledge  concerning  the  vegetation  category  which  each  measurement  vec¬ 
tor  actually  represents . 

The  Euclidean  distance- similarity  coefficient  from  the  multiple- linkage 
clustering  technique  were  found  to  discriminate  crop  types  nearly  as  well  as 
the  other  methods  of  categorization  (see  Figures  12  and  13) .  Little  discrimin¬ 
atory  potential  was  indicated  for  the  correlation  coefficient,  which  suggests 
that  reciprocity  holds  between  the  two  radar  cross-polarization  returns  and 
that  depolarization  of  the  radar  signal  from  crop  to  crop  is  so  slight  that  it 
cannot  be  used  as  a  discriminating  factor.  The  correlation  similarity  coef¬ 
ficient  was  not  used  in  later  analyses. 

In  an  attempt  to  glean  more  information  from  the  analysis  of  the  August- 
September  data,  category  (b)  (see  Figure  12)  was  subdivided  into  three  separ¬ 
ate  categories  'Figure  14) .  Note  that  category  (b^)  contains  mostly  grain  sor¬ 
ghum,  category  ^2)  contains  mostly  grain  sorghum  and  alfalfa,  and  category 
(bj)  contains  mostly  grain  sorghum  and  fall  wheat,  indicating  that  there  is 
some  discrimination  even  in  the  very  complex  portion  of  the  data  set.  Com¬ 
paring  this  grouping  with  Figure  6  shows  that  one  would  be  hard  pressed  to 
discriminate  regions  of  such  quality  from  the  scatterplot.  Partial  discrimin¬ 
ation  of  this  nature  leads  one  to  believe  that  if  data  were  available  in  other 
wavelengths  or  through  time,  better  discrimination  could  be  achieved. 

Obtaining  data  through  time  leads  to  the  question  of  the  best  time  to 
collect  the  data  if  better  discrimination  is  desired.  Looking  again  at  Figure 
14,  we  see  that  the  problem  lies  in  the  b^,  b^ ,  and  b^  categories.  To  arrive 
at  a  solution  to  the  problem  of  finding  the  optimum  tiine  to  gather  data  for 
crop  type  discrimination,  the  curves  in  Figure  15  v/ere  constructed.  The  solid 
lines  on  the  curves  are  based  on  actual  average  of  the  total  data  available 
at  the  University  of  Kansas  (July,  August,  September,  October,  and  November). 
We  have  based  the  dotted  lines  on  crop  growth  patterns  and  maturity  dates. 
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FIGURE  10.  Conditional  probability  clustering  for  August-September 
1965  radar  Imagery.  S  be  clusters  were  asked  for,  of  which  4  were  col¬ 
lapsed  together  since  they  had  the  most  similar  category  distributions 


together  since  they  had  the  most  similar  category  distributions . 
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SIX  CATEGORY  BREAKDOWN  USING  MULTIPLE  LINKAGE  CLUSTER  ANALYSIS 


FIGURE  14 
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Using  the  curves  to  determine  when  to  collect  data  which  would 
allow  better  discrimination  of  grain  sorghum  and  winter  wheat  (category 
bp,  we  see  that  data  taken  in  April  or  May  would  allow  us  to  achieve 
much  better  discrimination.  During  these  spring  months  the  grain  sorghum 
fields  are  being  prepared  for  the  planting  of  the  crop.  On  the  other  hand, 
winter  wheat  is  In  the  process  of  rapid  growth,  and  therefore  would  give 
a  much  higher  radar  return.  Likewise,  this  same  time  period  would  allow 
much  better  discrimination  between  grain  sorghum  and  alfalfa  (category  b2). 
Alfalfa,  like  winter  wheat,  is  in  the  early  rapid  growth  stage  during  the 
early  spring  months .  The  discrimination  of  alfalfa  from  winter  wheat  appears 
difficult  at  first  glance.  However,  the  curve  for  winter  wheat  in  Figure  15 
is  for  the  entire  winter-wheat  growing  season,  which  starts  in  August. 

Many  times  winter  wheat  will  not  be  immediately  replanted  in  the  crop  rota¬ 
tion  system;  that  is,  winter  wheat  will  be  planted  on  a  fallow  field  or  will 
be  followed  by  a  fallow  field.  Even  when  a  field  of  winter  wheat  is  replanted 
after  harvest,  discrimination  may  be  achieved  by  concentrated  sensing  in 
July  and  early  August,  during  the  bare-ground  phase  of  the  cycle. 

While  crop  discrimination  in  this  study  was  confined  to  radar 
imagery,  the  techniques  could  just  as  easily  be  applied  to  imagery  from 
other  remote  sensing  systems .  It  is  probable  that  the  best  discrimination 
would  be  obtained  by  combining  imagery  of  several  congruent  data  sets 
from  a  number  of  sensing  systems  taking  time  sequential  imagery  rather 
than  using  all  four  polarizations  in  a  single  time  period. 

Possible  Radar  Systems  and  Their  Problems 


The  results  of  the  study  presented  above  show  an  encouraging 
ability  to  discriminate  crops  with  radar,  especially  in  view  of  the  paucity 
of  data.  The  use  of  only  a  single-frequency  multiple-polarization  system 
at  various  times  and  for  various  polarizations  and  frequencies  is  urgently 
needed  before  an  operational  system  for  crop  identification  can  be  designed. 

Radar  systems  flown  to  date  for  agricultural  purposes  have  operated 
at  K-band  and  at  X-band.  Only  the  K-band  system  has  provided  the  multiple- 
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TIGURE  15,  Relative  radar  return  curves: 
The  relative  densities  are  based  on  HH 
polarization  K-Band  radar  image  densities. 
The  solid  lines  are  averages  of  the  return 
for  each  crop  type  based  upon  the  temporal 
data  available.  The  images  were  obtained 
over  a  three  year  time  period  and  therefore 
do  not  represent  the  average  return  for  a 
;  given  set  of  fields.  The  dashed  lines  were 
interpolated  using  known  crop  growth 
cycles,  maturity  dates  and  the  average 
return  from  categories  such  as  bare  ground 
and  wheat  stubble. 
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polarization  coverage  that  seems  so  useful.  As  indicated  previously/  the 
use  of  multiple  frequencies  should  have  many  of  the  advantages  for  identi¬ 
fication  of  crops  found  by  the  use  of  multiple  wavelengths  in  the  visible 
and  infrared.  The  exact  spread  of  frequencies  required,  however,  cannot 
be  determined  on  a  purely  theoretical  basis,  and  many  more  measurements 
will  be  necessary  before  a  full  specification  can  be  produced. 

Clearly  both  K-band  and  X-band  data  are  useful  in  crop  discrimin¬ 
ation,  For  crops  with  fairly  large  physical  structure,  such  as  corn,  wave¬ 
lengths  longer  than  the  3  cm  of  X-band  may  also  be  valuable.  Wavelengths 
as  long  as  10  or  20  cms  may  help  in  crop  discrimination,  and  even  longer 
wavelengths  may  be  required  for  discrimination  between  soils  or  between 
trees  in  forest  land.  Research  is  required  to  determine  both  the  spread 
and  number  of  frequencies  that  can  be  effectively  used. 

The  images  studied  here  were  produced  with  a  real-aperture  system. 
At  the  longer  wavelengths,  a  synthetic  aperture  system  will  be  necessary 
to  achieve  the  type  of  resolution  needed  in  crop  studies  from  aircraft. 
Spacecraft  platforms  operating  at  any  reasonable  frequency  will  require 
synthetic  aperture  systems  to  produce  imagery  of  usable  resolution. 

Synthetic-aperture  Images,  like  photographs  made'with  pure  laser 
light,  have  a  speckled  appearance,  because  each  point  is  likely  to  be 
governed  by  a  Rayleigh  distribution  and  to  be  only  one  element  of  that 
distribution.  Averaging  over  many  points  produces  the  mean  value  which 
is  most  significant  in  the  recognition  systems.  These  averages  may  be 
achieved  spatially,  but  they  may  also  be  achieved  in  frequency  by  trans¬ 
mitting  a  continuous  spectrum  in  the  vicinity  of  each  of  the  frequencies 
of  a  mu  It  I- frequency  system.  In  fact,  it  may  be  desirable  to  sweep  across 
an  entire  band  of  frequencies — for  example,  from  8  to  12  GHz,  if  it  turns 
out  in  the  experiments  that  there  is  sufficient  differentiation  of  crops 
within  this  band.  This  technique  is  being  tested  in  the  1  to  8  GHz  band 
at  the  University  of  Kansas  at  present,  although  no  measurements  have 
yet  been  made. 

Although  we  clearly  lack  a  scientific  basis  for  describing  a  system 
that  might  not  be  operational  for  ten  years,  the  meager  evidence  to  date 
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Indicates  that  the  system  might  look  something  like  this:  an  aircraft 
system  carrying  a  radar  operating  at  either  fixed  or  continuously-swept 
frequencies  over  a  range  of  perhaps  3  to  30  GHz.  The  radar  would  look 
to  the  side  and  produce  synthetic-aperture  images.  Information  would  be 
recorded  on  a  multi-channel  tape  recorder.  The  recorded  signals  would  be 
processed  by  a  suitable  congruencing  technique  so  that  they  could  be 
compared  not  only  with  each  other,  but  with  signals  from  the  same  area 
taker,  in  previous  flights.  The  measurement  vectors  so  produced  would 
be  processed  through  one  of  the  techniques  described  here  (or  others  to 
be  developed)  to  produce  outputs  of  value  to  the  agriculturalist.  The  out¬ 
puts  would  include  maps  in  which  individual  fields  would  be  located,  along 
with  their  crops  and  perhaps  some  indication  of  the  average  measurement 
vector  for  each  of  the  crops.  In  addition,  the  system  would  be  capable 
of  producing  integrated  outputs  giving  total  areas  of  various  crops  and 
various  stages  of  growth.  Whether  the  system  will  ever  be  able  to  point 
to  diseased  crops  and  actually  to  give  predictions  of  yields  over  major 
areas  remains  to  be  determined.  That  such  a  system  could  indeed  discri¬ 
minate  crops  seems  plausible  on  the  basis  of  evidence  at  hand,  but  this 
ability  certainly  has  not  been  proven. 

An  aircraft  system  such  as  this  could  tiy  at  perhaps  500  miles  per 
hour  to  obtain  continuous  coverage  in  the  United  States  in  a  reasonable 
length  of  time.  The  swath  width  might  be  20  to  40  miles,  and  the  area 
covered  would  be  determined  accordingly. 

Small  satellites  can  carry  single-frequency,  single -transmitted 
polarization  systems  at  the  present  state  of  the  art  (Moore  and  Simonett, 
1967)t  Such  systems  could  give  coverage  over  the  entire  world  and 
would  be  capable  of  the  sorts  of  discriminations  described  here,  although 
repetitive  coverage  might  improve  the  results  over  the  rather  meager  cover¬ 
age  used  in  the  present  analysis.  With  larger  spacecraft  capable  of 
supplying  perhaps  a  kilowatt  of  power,  the  full  broad-band  multipolarization 
system  could  be  crrried  in  space.  The  swath  width  to  be  used  in  space 
will  depend  upon  the  antenna  size  and  the  power  available ,  and  will  vary 
from  a  minimum  of  around  20  kilometers  with  the  systems  that  could  be 

GRES  Technical  Report  61-32  (see  Bibliography,  p.  G4) 
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carried  on  a  Nimbus-size  satellite  at  present  to  much  larger  numbers  pos¬ 
sible  if  more  power  and  space  are  available. 

The  potential  of  this  technique  is  clearly  high  for  crop  studies  the 
world  over,  but  many  of  the  predictions  are  based  upon  extreme  extrapola¬ 
tions  from  rather  limited  data.  Obviously,  a  great  deal  more  experimenta- 
Uon  is  called  for,  particularly  with  radars  operating  in  different  wavelengths 
and  with  multiple  coverage  at  frequent  Intervals  throughout  the  growing 
season.  Furthermore,  extensive  studies  at  present  have  been  limited  to 
one  site  in  western  Kansas  and  to  another  in  the  Phoenix  area,  where  some 
work  has  been  done  by  Goodyear  Aerospace  Corporation.  These  sites  are 
representative  of  certain  cropping  practices  and  certain  types  of  crops, 
but  flights  must  obviously  be  arranged  in  other  areas  of  the  world  where 
the  crops,  agricultural  practices,  or  both  are  different  from  those  observed 

in  the  United  States . 


Concluding  Remarks 

From  the  preceding  discussion  the  following  conclusions  emerge. 

(1)  In  any  remote  sensing  system  concerned  with  the  detection  of 
and  discrimination  between  crops,  there  is  a  severe  fall-off  in  the  ability 
of  the  system  to  provide  the  answers  requested  of  it,  when  constraints  of 
Increasing  severity  are  placed  on  the  sensor,  the  timeliness  of  data  collec¬ 
tion,  the  congruencing  of  data  sets  through  time,  or  the  spatial  or  data- 
discrimination  requirements  of  the  system.  Thus,  rather  than  asking  too 
much  of  a  sensor  system  and  suffering  disappointment  because  it  cannot 
provide  the  required  information,  it  is  aoproprlate  to  define  for  each  remote 
sensing  problem  various  alternative-constraints  strategies  and  loss  funcllans 

so  that  rational  choices  may  be  made. 

(2)  Though  only  a  small  amount  of  data  were  available,  we  have 

clearly  shown  that  identification  accuracies  In  the  ninetieth  percentile  can 
be  attained  by  using  the  methods  Illustrated  in  this  paper.  Further  work 
needs  to  be  done  on  clustering  techniques  v/hlch  can  analyze  data  in  as 
detailed  a  fashion  as  the  multiple-linkage  clustering  procedure,  yet  with 
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the  ability  of  the  conditional-probability  clustering  method  to  handle  large 
amounts  of  data  quickly.  No  attempt  has  been  made  here  to  associate 
crop  type  with  the  spatial  texture  distribution  it  produces  on  the  image, 
yet  this  might  well  alia  /  for  identification  accuracies  approaching  99%. 
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OBSERVATIONS  ON  THE  GEOMORPHOLOGY  AND  LAND  1!fSE 
OF  PART  OF  THE  WASATCH  RANGE,  UTAH 

R.  M.  Peterson 

AN/APQ-97  multiple  polarization  radar  imagery  obtained  in  October, 
1965  on  NASA  flight  number  90  was  examined  to  document  its  potential  for 
obtaining  the  following  types  of  information:  physiography,  surface  mat¬ 
erial,  geomorphic  history,  mineralized  areas,  drainage,  land  use,  and 
transportation . 

Both  HH  and  HV  radar  polarizations  were  examined  with  the  naked 
eye  and  with  a  5  power  microscope.  2X  enlargements  of  the  radar  imagery 

were  also  examined;  sometimes  with  a  reading  glass. 

Thorough  study  of  the  radar  imagery  resulted  in  a  wealth  of  informa¬ 
tion,  much  of  which  is  new  and  which  does  not  appear  in  reports  or  maps 
of  the  area.  Examples  are  some  ancestral  stream  patterns  with  resultant 
information  on  the  origin  of  the  Wasatch  Range,  fracture  zones  and  possibly 
undiscovered  mineralized  areas . 

After  a  study  of  the  radar,  information  obtained  from  the  imagery  was 
compared  with  data  on  topographic  maps,  geologic  maps,  and  1:63,360 
scale  aerial  photo  mosaics .  It  is  significant  that  even  though  the  scale  of 
the  radar  imagery  is  only  approximately  1:160,000,  it  was  necessary  to 
consult  maps  of  1:24,000  scale  to  observe  some  of  the  same  features.  In 
some  cases  comparison  of  first-round  radar  data  with  published  data  gave 
clues  to  more  features  to  search  for  on  radar  imagery  and  resulted  in  new 
discoveries.  Examples  are  two  previously  unreported  major  fracture  zones 
and  a  joint  set  that  appear  to  control  location  of  mineral  deposits  in  the 

Park  City  area. 

The  following  maps  were  prepared  by  tracing  information  from  the 
radar  imagery  onto  transparent  overlays:  physiographic  maps  (2  types) , 
landform  classification  maps,  surface  material,  ancestral  drainage, 
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present  drainage ,  land  use,  and  transportation. 

The  study  area  is  a  rectangle  approximately  12  by  30  miles  extending 

east  from  Salt  Lake  City  to  the  east  side  of  Rhodes  Valley,  Utah  (see  Fig.  1). 
The  major  portion  of  the  study  area  is  occupied  by  the  Wasatch  Mountains 

whose  steep  escarpment  overlooks  Salt  Lake  City, 

The  Wasatch  Range,  the  westernmost  range  of  the  Middle  Rocky 
Mountains,  trends  north  and  south  and  lies  at  the  east  side  of  the  Great 
Basin.  The  local  relief  between  the  Great  Salt  Lake  Valley  and  the  highest 
part  of  the  range,  which  is  about  two  miles  from  the  Wasatch  Front,  is 
over  6,800  feet  (the  highest  peaks  in  this  area  are  Mt.  Olympus  at  9,026 
feet  and  Twin  Peaks  at  11,330  feet),  (Fig.  1).  The  drainage  divide  lies 
approximately  10  miles  east  of  the  front  and  at  a  lower  elevation  than  peaks 
near  the  front.  There  is  a  gentle  eastward  slope  of  the  Wasatch  Range  from 
near  the  front  to  the  broad  valleys  of  Parleys  Park  and  Rhodes  Valley,  which 

are  in  the  Wasatch-Uinta  transition  area. 

Rhodes  Valley  (also  known  as  Kamas  Valley  or  the  Kamas  Sag)  lies 

at  the  extreme  western  end  of  the  Uinta  Range  and  is  a  wide,  flat-floored 
valley  at  an  elevation  of  approximately  6,400  feet.  To  the  west  of  Rhodes 
Valley  the  low,  rolling  West  Hills  are  covered  predominately  by  sagebrush  and 
scrub  oak  and  are  usable  pri.marily  for  grazing .  The  only  other  flat  areas  in  the 
study  area  are  Parley's  Park,  at  an  elevation  of  6,600  feet,  and  the  narrow 
flood  plains  and  terraces  along  the  Provo  and  the  Weber  Rivers  and  their 

tributaries . 

Natural  Vegetation 

The  natural  vegetation  on  the  floors  of  Parleys  Park  and  Rhodes 
Valley  was  grass;  cottonwood,  river  birch,  willow  and  box  elder  trees 
grew  along  streams.  The  lower  hills  are  covered  by  sagebrush  and  some 
juniper,  June  grass  and  remnants  of  the  native  Indian  rice  grass.  Scrub 
oak  and  sagebrush  with  snowberry  and  some  maple  grow  between  7,000  and 
8,000  feet.  At  8,000  to  10,000  feet  are  coniferous  forests  of  Engleman 
spruce  with  some  Alpine  fir,  Douglas  fir,  and  lodgepole  pine.  The  higher 
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parts  of  the  study  area  are  above  timber  line  which  ranges  from  9,500  to 
10,500  feet  (Hawkes ,  1959,  pp.  9-11). 

Drainage 

Several  streams  have  cut  scenic  canyons  into  the  Wasatch.  Within 
the  study  area  most  of  these  stream  valleys  and  the  higher  parts  of  the  range 
are  covered  by  dense  forests  and  have  snow  for  several  months  each  year. 
Because  of  the  spectacular  alpine  scenery,  snow  cover  and  proximity  to 
Salt  Lake  City,  three  nationally  famous  ski  resorts  are  located  within  the 
study  area.  These  are  the  Alta  at  the  head  of  Little  Cottonwood  Canyon, 
Brighton  at  the  head  of  Big  Cottonwood  Canyon,  and  the  Park  City  3ki  area 
near  the  mining  town  of  Park  City,  (see  Fig.  1). 

Drainage  of  the  entire  study  area  is  to  the  Great  Salt  Lake.  Several 
canyons,  cut  through  the  Wasatch  Front  by  westflowing  streams,  are  visible 
on  the  imagery.  From  north  to  south  these  are  Parleys,  Mill  Creek,  Neffs, 
Big  Cottonwood,  Little  Cottonwood,  and  Bells  Canyons.  The  heads  of  most 
of  these  canyons  contain  cirques  from  which  Pleistocene  valley  glaciers 
advanced  down  the  canyons,  in  cases  as  far  as  the  Wasatch  Front. 

The  east  side  of  the  region  is  drained  by  two  rivers  that  head  in  the 
Uinta  Mountains.  The  Weber  River  has  its  origin  on  the  north  side  of  the 
Uintas  and  flows  west  across  the  north  side  of  Rhodes  Valley,  then  in  its 
northward  course  it  receives  tributaries  which  drain  most  of  the  central 
part  of  the  study  area.  The  Provo  River  heads  on  the  south  flank  of  the 
Uintas,  flows  through  the  south  end  of  Rhodes  Valley,  through  the  narrow 
Upper  Provo  Canyon  across  Heber  Valley  and  then  to  Great  Salt  Lake  via 
Utah  Lake  and  the  Jordon  River. 

Land  Use 

The  land  use  in  this  area  depends,  to  a  large  extent,  on  elevation, 
slope  of  the  land,  availability  of  water,  and  location  of  mineral  deposits. 
Land  use  of  the  study  area  is  shown  in  Figure  2  which  was  prepared  from 
radar  imagery.  Mountainous  areas  here  serve  five  main  uses:  One  is  the 
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function  of  watershed  and  storage  area  for  water.  This  is  especially  impor¬ 
tant  for  the  arid  Wasatch  Oasis  which  extends  along  the  west  side  of  the 
Wasatch  Range.  Bailey  (1941,  p.  194)  states  that  at  least  80%  of  the  usable 
stream  flow  in  Utah  comes  from  mountain  lands  over  7,000  feet  in  altitude 
and  that  each  acre  of  irrigated  land  is  dependent  upon  water  from  approxi¬ 
mately  seven  acres  of  watershed  land.  The  Wasatch  Oasis,  which  contains 
most  of  Utah's  farms,  cities,  and  population,  would  not  exist  if  there  were 
no  watersheds  such  as  those  in  this  study  area.  Grazing  is  a  second  use 
of  the  mountain  land.  Large  herds  of  sheep  and  some  beef  cattle  are  grazed 
in  many  parts  of  the  area  although  the  hills  between  Rhodes  Valley  and 
Parleys  Park  are  most  heavily  grazed.  Forestry  is  a  third  use  of  these 
mountains  although  this  is  not  as  important  as  it  was  in  earlier  times. 

Mining,  a  fourth  land  use,  has  been  conducted  near  Park  City  and  in 
Little  Cottonwood  Canyon  since  1868.  Over  450  million  dollars  of  gold, 
silver,  copper,  lead,  and  zinc  have  been  mined  in  the  Park  City  Mining 
District,  although  many  mines  have  closed  in  recent  years.  A  fifth  land 
use  in  the  mountainous  areas  is  recreation.  Ski  resorts  here  draw  large 
numbers  of  people  from  all  over  the  United  States  each  year.  In  the  canyons 
east  of  Salt  Lake  City  there  are  many  weekend  cottages  and  summer  homes 
in  addition  to  numerous  camps  and  picnic  facilities . 

The  predominant  use  of  land  in  the  valleys  is  for  agriculture .  The 
majority  of  the  farmers  live  in  nucleated  settlements  which  are  readily 
apparent  on  the  radar  imagery.  In  Rhodes  Valley  farmsteads  not  in  the 
villages  tend  to  be  on  the  main  north-south  road.  The  imagery  shows  the 
relationship  between  well  drained  land  and  the  location  of  the  villages  and 
farms;  most  of  the  farms  are  on  the  better  drained  sides  of  Rhodes  Valley  and 
Parleys  Park.  The  lower  sides  of  the  valleys  have  stronger  return  from  the 
wet  land  covered  with  marsh  grass. 

Inspection  of  the  radar  imagery  shows  many  clusters  of  cottages , 
mine  buildings  and  mine  dumps.  Examination  of  both  HH  and  HV  polarizations 
reveals  more  cultural  features  than  either  one  polarization  alone.  For 
example,  several  mine  dumps  appear  on  the  imagery.  Near  Keatley  south 
of  Park  City,  mine  dumps  and  railroads  leading  to  the  mines  are  visible  on 
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imagery.  South  of  Parleys  Park  the  waste  dumps  on  Quarry  Mountain  can 
be  seen. 

By  comparing  both  polarizations  and  combining  information  from  each, 
it  is  possible  to  prepare  a  transportation  map  from  radar  imagery  (Fig.  3). 

Shown  are  only  the  hard-surface  main  roads,  not  jeep  roads  which  are  visible 
in  many  places . 

Geomorphic  history.  The  geomorphology  of  the  present  Wasatch 
Pange  is  the  result  of  a  fairly  complex  history.  As  summarized  by  Granger 
,  (1952)  and  Threet  (1959)  the  following  sequence  occurred.  At  the 
beginning  of  the  Laramide  orogeny,  in  middle  and  late  Cretaceous  time, 
north-south  compression  formed  the  east-west  folds  of  the  Central  Wasatch. 
Thrust  sheets  from  the  west  overrode  the  site  of  the  present  Wasatch  and 
elastics  eroded  from  highlands  were  deposited  in  adjacent  basins.  After 
a  quiet  period  the  Cretaceous  folds  were  further  compressed  in  Paleocene 
or  Eocene  time;  at  about  the  same  time  the  Little  Cottonwood,  Alta,  Clayton 
Peak  and  Park  City  stocks  were  intruded. 

The  Wasatch  Range  was  first  elevated  by  folding,  along  north-soviUi 
axes,  which  began  in  the  late  Eocene  or  Oligocene .  At  this  time  the  major 
relief  features  were  formed,  the  courses  of  the  Weber  and  Provo  Rivers  were 
established  and  lavas  and  andesitic  pyroclastics  were  extruded. 

Regional  warping  along  old  trends  in  late  Oligocene  or  early  Miocene 
was  followed  by  erosion  of  broad  surfaces.  One  of  these  broad  erosional 
surfaces,  the  Weber  Valley  surface,  may  be  preserved  in  Big  Cottonwood 
Canyon  (Granger,  et^.  1952,  p.  36)  where  there  are  elevated,  broad 
valley  profiles.  The  former  erosion  surface,  represented  by  an  extension 
of  the  spurs  on  the  sides  of  Big  Cottonwood  Canyon  has  been  dissected. 

This  erosion  probably  began  before  Pleistocene  glaciation,  possibly  in  late 
Pliocene  time  (Granger  et  al,. ,  1952,  p.  36). 

Block  faulting  and  broad  warping,  which  began  in  late  Miocene  or 
Pliocene,  was  continued  in  late  Pliocene  and  Pliestocene.  Associated  with 
the  uplift  was  dissection  of  old  erosional  surfaces,  div  ''rsion  and  capture  of 
streams,  and  erosion  of  Parleys  Park  and  Rhodes  Vaik>  During  the  Pleistocene 


the  high  parts  o£  the  Wasatch  were  glaciated  and  Lake  Bonneville,  the  fresh 
water  ancestor  of  Great  Salt  Lake,  filled  part  of  the  Great  Basin. 

n.<.nmnrnhlc  Regions .  The  area  studied  can  be  divided  into  the 
following  regions:  (1)  The  small  portion  of  the  Great  Basin  where  Salt  Lake 
City  is  situated,  (2)  The  mountainous  terrain  of  the  complex  Wasatch  Range, 

(3)  The  valleys  of  Parleys  Park  and  Rhodes  Valley,  and  (4)  The  Andesite 
Hills.  The  Uinta  Range  is  not  discussed  in  this  paper  because  only  the 
western  tip  of  this  range  is  shown  cn  the  imagery. 

1.  ThP  q.U  Lake  CitvAre_a.  Shown  on  the  radar  imagery  is  part  of 

salt  Lake  City  and  small  neighboring  cities.  Several  stream  terraces  are 
visible  where  the  distributaries  from  Big  and  Little  Cottonwood  Canyons 
have  deepened  their  channels.  Much  of  Salt  Lake  City  is  built  on  the  floor 
of  Pleistocene  Lake  Bonneville.  The  waves  of  Lake  Bonneville,  which  was 

approximately  the  size  of  present  day  Lake  Superior,  cut  ‘  " 

Wasatch  Front.  A  delta  built  into  Lake  Bonneville  at  the  mou  h  of  Big  Cotton 
wood  canyon  is  seen  on  the  imagery.  Also  shown  is  a  lateral  moraine  at  the 
mouth  of  Little  Cottonwood  Canyon  and  an  end  moraine  at  the  °  ^ 

canyon.  Fault  scarps  are  shown  which  cut  across  the  moraines  and  delta 

between  Bell  and  Big  Cottonwood  Canyons. 

2  Th,  lAissstch  Range.  This  subregion  includes  the  Wasatch  Range 

from  the  Wasatch  Front  to  Parleys  Park  and  the  spur  of  the  Wasatch  Range 

south  of  Park  City  and  east  to  the  Provo  River.  ,  , 

The  complex  structure  consists  mostly  of  sediments  folded  along 
east-west  axes:  these  folds  are  truncated  on  the  west  by  the  Wasatch  fault, 
on  the  south  side  of  the  study  area  a  series  of  plutons  has  intruded  a  large 
Ist-west  anticlinal  fold  which  seems  to  be  the  continuation  of  the  Uinta  fold 
Thrust  faults  and  high  angle  faults  occur  throughout  the  area . 

The  topography  of  the  Wasatch  Range  is  youthful.  Most  slopes  ar 
near  the  angle  of  repose.  Stream  valleys  are  characterized  as  having  V- 
shaped  profiles  in  the  non-glaciated  portions,  steep  gradients,  and  valley 
floors  barely  wider  than  the  streams  in  them.  Road  building  in  these  canyons 
has  always  been  a  problem  because  of  the  narrowness  of  the  canyons  an 
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the  almost  vertical  cirque  walls  at  the  heads  of  canyons.  Hard  surfaced 
roads  extend  up  the  canyons  on  the  west  side  of  the  range,  but,  other  than 
jeep  roads,  no  roads  (or  railroads)  cross  the  area  shown  on  this  imagery. 

The  radar  Imagery  shows  that  the  drainage  of  the  Wasatch  Range  is 
well  developed.  No  uplands  remain  other  than  the  narrow  ridge  divides. 
Several  stream  patterns  with  a  variety  of  controls  can  be  recognized  on  the 
imagery.  The  drainage  pattern  of  the  major  v;est-flowing  streams  is  basic¬ 
ally  dendritic  but  some  angularity,  caused  by  faults  and  joints,  can  be 
distinguished  in  stream  orientations.  The  attitude  of  strata  influences 
some  stream  orientations;  this  is  especially  evident  in  Thaynes  Canyon, 

Iron  Canyon,  and  Dutch  Draw  south  of  Parleys  Park. 

The  influence  of  structure  on  drainage  is  evident  by  the  stream  pattern 
around  the  Alta  Stock  and  Clayton  Peak  stock. 

The  topography  of  the  higher  parts  of  the  Wasatch  can  best  be  des¬ 
cribed  as  alpine.  The  high  Wasatch  Range  is  glaciated  to  such  an  extent 
that  one  of  the  most  striking  features  observable  on  radar  imagery,  or  on 
geologic  maps ,  is  the  large  proportion  of  the  area  covered  by  glacial  moraine . 
Glacial  moraine  is  conspicuous  on  radar  imagery  because  of  its  character¬ 
istic  graytone  a^  speckled  texture.  Moraine  covers  the  cirque  floors,  and 
in  some  cases  extends  for  miles  down  the  canyons. 

Shown  on  the  imagery  are  many  slpine  features  such  as  U-shaped 
valleys,  end  moraines,  lateral  moraines,  ground  moraines,  cirques,  tarns, 
horns,  cols,  and  truncated  spurs.  Several  hanging  valleys  can  be  seen  in 
Big  Cottonwood  Canyon.  A  landslide  is  shown  on  the  imagery  in  Big  Cotton¬ 
wood  Canyon. 

3 .  Parleys  Park  and  Rhodes  Valley  .  Parleys  Park  is  an  alluvial- 
filled  valley  at  approximately  6,600  feet  elevation.  The  valley  is  almost 
level  with  a  gentle  slope  to  the  northeast.  Much  of  the  valley,  especially 
the  east  side,  is  swampy  and  covered  with  marsh  grass  as  indicated  by  dark 
gray  tone  on  the  radar  imagery.  On  Quarry  Mountain  at  the  south  side  of 
the  valley,  waste  dumps  from  large  quarries  can  be  seen  on  the  imagery. 

Rhodes  Valley  is  a  wide  alluvial  valley  with  several  nucleaded  farming 
villages.  River  terraces  are  visible  along  both  the  Provo  River  and  the  Weber 
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River.  The  andesite  Hills,  locally  called  the  West  Hills,  lie  to  the  west 
of  Rhodes  Valley. 

4.  The  Andesite  Hills  .  The  Andesite  Hills  lie  in  the  eastern  part 
of  the  study  area  between  Parleys  Park  and  Rhodes  Valley  and  south  of 
Rhodes  Valley.  The  higher  hills  lie  at  about  7,800  feet  elevation  compared 
with  6,400  in  Rhodes  Valley.  Most  slopes  are  moderate,  although  in  the 
western  portion  there  are  many  gentle  slopes.  Local  relief  is  generally 
moderate;  less  than  500  feet.  The  region  is  characterized  as  being  in  a 
mature  stage  of  an  erosion  cycle. 

The  andesitic  pyroclastic  rock  of  early  Tertiary  age  (Geologic  Map 
of  Northeastern  Utah,  1961 ,  compiled  by  Stokes  and  Madsen)  has  a  charac¬ 
teristic  graytone  and  texture  on  radar  imagery  plus  a  characteristic  erosional 
pattern.  The  vegetation  cover  of  grass,  sagebrush  and  scrub  oak,  charac¬ 
teristic  of  these  hills,  seems  to  be  responsible  for  some  of  the  signal  return. 

The  West  Hills  (generally  the  region  west  of  Rhodes  Valley  and  east 
of  Parleys  Park,  Richardson  Flat  and  Keetly)  have  a  drainage  divide  near  the 
west  side  of  the  region  so  the  streams  flowing  east  are  long  and  have  gentle 
gradients  compared  with  the  streams  on  the  west  side  of  the  region.  The 
drainage  pattern,  generally  dendritic,  is  an  indication  of  the  homogeneous 
andesite . 
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Rock  Units 

Rocks  within  the  study  area  range  in  age  from  Precambrian  to 
Quaternary.  Following  is  a  generalized  list  of  stratigraphic  units  in  the 
Wasatch  Mountains  and  vicinity. 

Igneous  Rocks 

Early  Tertiary?  Andesitic  pyroclastic  rocks 
Diorite  porphyry;  Park  City  stock 
Diorite;  Clayton  Peak  stock 
Granodiorite;  Alta  stock 
Quartz  monzo'iite;  Little  Cottonwood  stock 

Quaternary 

Alluvial  deposits 
Glacial  deposits 

Tertiary 

Fowkes  Formation 
Knight  Conglomerate 
Crazy  Hollow  Formation 

Cretaceous 

Echo  Canyon  Conglomerate 
Wanship  Formation 
Frontier  Sandstone 
Aspen  Shale 
Bear  River  Formation 
Kelvin  Conglomerate 

Jurassic 

Preuss  Sandstone 
Twin  Creek  Limestone 
Nugget  Sandstone 

Triassic 

Ankareh  Formation 
Thaynes  Formation 
Woodside  Shale 

Permian 

Park  City  Formation 

Pennsylvanian 

Weber  Quartzite 
Round  Valley  Limestone 


no 


Mississippian 

Doughnut  Formation 
Humbug  Formation 
Deseret  Limestone 
Gardison  Limestone 
Fitchville  Formation 

Cambrian 

Maxfield  Limestone 
Ophir  Formation 
Tintic  Quartzite 

Precambrian 

Mutual  Formation 
Mineral  Fork  Tillite 
Big  Cottonwood  Formation 

Many  of  the  above  rock  units  can  be  recognized  and  traced  on  radar 
imagery.  Examples  are  the  Weber  Quartzite,  the  Thaynes  Formation, 

Ankareh  Formation,  the  Nugget  Sandstone,  and  various  igneous  rocks. 

An  extensive  area  of  volcanics  south  and  west  of  Rhodes  Valley 
which  are  mapped  as  early  Tertiary?  andesitic  pyroclastics  (Geologic  Map 
of  Utah,  compiled  by  Stokes  and  Madsen,  1961)  are  recognizable  on  radar 
imagery  by  characteristic  gray  tones  and  textures,  plus  the  group  of  char¬ 
acteristics  one  would  use  to  recognize  the  same  type  of  rocks  on  aerial 
photographs . 

An  especially  interesting  feature  of  imaging  radar  is  its  ability  to 
delineate  glacial  moraines.  Not  only  is  ground  moraine  readily  apparent, 
but  so  also  are  lateral  moraines  (in  Bonanza  Flats,  for  example)  and  terminal 
moraines  (examples  in  Thaynes  Canyon,  Dutch  Draw,  and  at  the  mouth  of 
Bells  Canyon) .  The  interface  between  moraines  and  alluvium  is  very  appar¬ 
ent  (for  example  in  Thaynes  Canyon  and  Dutch  Draw).  General  boundaries 
of  the  alluvium  which  fills  river  valleys  and  small  basins  (as  at  Parley 
Park,  Rhodes  Valley,  and  the  valley  of  the  Provo  River)  can  be  readily 
mapped  from  radar  imagery. 

Figure  4,  which  is  a  generalized  map  showing  surface  materials, 
was  prepared  by  tracing  boundaries  of  alluvium,  moraine,  and  rock  classes 
from  the  radar  imagery  onto  a  transparent  overlay . 
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Structure 

In  the  west  part  of  the  study  area  the  major  folds  of  the  Wasatch 
have  east-west  axes.  The  north  limb  of  the  large  anticlinal  fold  that 
appears  to  be  the  westward  extension  of  the  Uinta  Arch  can  be  seen  in 
Mill  Creek  Canyon  and  Big  Cottonwood  Canyon.  The  south  limb  of  this 
anticline  is  now  mostly  part  of  the  Little  Cottonwood  stock.  Steep  north¬ 
ward  dips  and  northeast- striking  beds  are  visible  on  the  imagery  north 
of  Big  Cottonwood  Canyon.  Inasmuch  as  Little  Cottonwood  Canyon  lies 
at  the  outer  range  of  the  radar  scan  where  the  radar  beams  had  a  low  inci¬ 
dent  angle,  the  high  north  wall  of  Little  Cottonwood  Canyon  caused  a 
radar  shadow  in  the  lower  end  of  the  canyon  making  this  portion  of  the 
Imagery  unusable.  A  few  miles  up  the  canyon,  however,  the  imagery  shows 
rock  that  has  erosional  characteristics  of  massive,  Igneous  rock.  This 

agrees  with  geologic  maps  which  show  the  Little  Cottonwood  granitic 
stock  at  this  location. 

In  addition  to  the  Little  Cottonwood  granite  stock,  already  mentioned, 
the  Alta  stock,  the  Clayton  Peak  stock,  and  the  Park  City  stock  are  dis- 
cernable  on  the  Imagery.  Since  these  are  mineralized  regions  which  have 
had  considerable  economic  importance,  the  ability  to  recognize  such  regions 
on  radar  is  very  significant. 

Many  of  the  folds  in  the  study  area  can  be  observed  on  radar  imagery. 
Because  the  attitudes  of  strata  can  be  ascertained  from  radar  imagery,  the 
orientation  of  major  folds  can  be  determined  from  imagery.  An  example  is 
the  northeast  trending  anticline  of  Park  City. 

Several  faults  in  the  Wasatch  Range  can  be  detected  by  methods 
similar  to  those  used  for  fault  detection  on  aerial  photographs.  In  some 
cases  abrupt  changes  in  gray  tones  and  textures  characteristic  of  rock  units 
suggest  faulting.  When  lithologic  units  are  recognized  on  radar,  offsetting, 
abrupt  termination,  or  abrupt  changes  in  altitudes  of  rock  units  can  be 

used  as  evidence  for  faulting.  Faults  with  physiographic  expression  are 
most  easily  detected. 

The  most  conspicuous  feature  on  the  Wasatch  Range  is  the  Wasatch 
Front  fault  escarpment.  There  is  relative  relief  at  the  front  of  6810  feet 
(4520  feet  elevation  near  Sandy,  opposite  the  mouth  of  Little  Cottonwood 
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Canyon;  11,330  feet  at  Twin  Peaks  north  of  Little  Cottonwood  Canyon) . 

In  addition  to  the  numerous  faceted  spurs  along  the  front,  the  imagery 
also  shows  recent  fault  scarps  that  cut  across  terraces  and  deltas  asso¬ 
ciated  with  Lake  Bonneville . 

Along  the  Wasatch  Front  in  Parleys  Canyon  and  Mill  Creek  Canyon 
the  imagery  shows  sediments  striking  generally  northeast  and  dipping 
steeply  north.  In  Neffs  Canyon  changes  in  strike  that  are  discernable 
on  radar  correspond  to  faults  shown  on  the  geologic  map  by  Granger  e^ 
(1952).  Attitudes  of  beds  are  apparent  at  many  other  places  on  the 
imagery , 

In  a  geomorphic  study  of  a  mountainous  region  such  as  the  Wasatch 
Range ;  structure  is  a  chief  element  in  the  development  of  the  geomorpho¬ 
logy.  It  is  therefore  necessary  to  know  the  location  and  nature  of  at 
least  the  major  faults  and  the  dominant  joint  patterns  in  order  to  under¬ 
stand  the  geomorphology.  A  few  other  examples  of  faults  detectable  by 
radar  are  given  below. 

In  the  upper  part  of  Mill  Creek's  Porter  Fork  the  Lower  Strand 
Mountain  branch  of  the  Mount  Raymond  thrust  fault  can  be  discerned  by 
a  chaijcje  in  altitude  of  beds  and  by  lineation  along  the  fault 
To  the  i^ast  the  Mount  Raymond  thrust  can  be  traced  by  the  position  of 
the  Weber  Quartzite  thrust  slices,  which  lie  above  the  thrust  planes,  and 
by  lineations  along  the  faults.  The  Silver  Fork  normal  fault  is  very  appar¬ 
ent  where  it  terminates  the  Thaynes  and  Arkereh  Formations  and  the  Nugget 
Sandstone  in  Reynolds  Gulch  north  of  Big  Cottonwood  Canyon.  This  fault 
can  be  traced  on  the  imagery  north  to  the  Mount  Raymond  thrust  which 
was  offset  by  the  Silver  Fork  fault.  Ground  truth  confirming  this  observa¬ 
tion  is  reported  by  Granger  ej^  ^. ,  (1952,  p.  24). 

In  Mill  D  Canyon,  a  Big  Cottonwood  tributary  of  Reynolds  Gulch, 
the  imagery  shows  a  zone  of  north-south  faulting  that  terminates  the 
Paleozoic  sediments  to  the  east.  The  sediments  east  of  Mill  D  Canyon, 
which  are  for  the  most  part  of  Mississippian  Deseret  and  Madison  Lime¬ 
stones,  exhibit  an  east-west  joint  pattern  that  is  very  pronounced  on  radar. 
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Possible  faults  or  fracture  zones  which  do  not  appear  on  published 
geologic  maps  are  suggested  by  lineations  on  the  radar  imagery.  South  of 
Parleys  Park  a  lineation,  trending  generally  northeast- southwest,  extends 
from  the  west  side  of  Quarry  Mountain  southward  up  the  west  side  of  Iron 
Canyon  and  along  the  west  side  of  the  intrusion  immediately  east  of 
Brighton.  Dikes,  mines  and  prospect  pits  north  of  this  intrusion  lie 
along  this  zone.  Another  possible  fault  or  fracture  zone,  suggested  by 
a  lineament  on  the  imagery,  lies  along  the  east  side  ol  the  aforementioned 
stock.  This  lineament  extends  from  the  west  sld^  ol  Clayton  i'^eak  notch- 
eastward  near  the  east  side  of  the  Alta  Stock. 

The  Park  City  Mining  District 

Park  City  is  one  of  the  great  base  metal  mining  districts  of  Utah 
(Wilson,  1954,  p.  182)  and  has  been  one  of  the  most  important  mining 
districts  in  the  United  States.  Silver,  gold,  lead,  zinc  and  copper 
are  the  chief  minerals  mined  there .  The  large  mine  dumps  are  visible 
on  radar  imagery,  especially  when  HH  and  HV  polarizations  are  compared. 

The  mineralized  zones  at  Park  City  are  related  to  the  intrusive  stocks 
of  quartz  monzonite  which  caused  contact  metamorphism  in  the  adjacent 
sediments,  especially  in  the  limestones.  (Lindgren,  1933,  p.  584), 

Lode  deposits  are  often  between  limestone  and  quartzite,  but  may 
be  in  quartzite.  Bedded  ores  are  sulphides  that  have  replaced  limestones 
(the  Park  City  Formation  or  Thaynes  Formation) . 

An  examination  of  the  Park  City  West  Geologic  Map  (Crittenden, 
et^. ,  1966)  shows  that  mines  tend  to  be  at  contacts  between  the  Woodside 
shale  and  Thaynes  Formation  or  between  the  Weber  Quartzite  and  the  cal¬ 
careous  Park  City  Formation.  The  mineralizing  solutions  apparently  moved 
along  channels  formed  by  faults  and  joints.  The  Park  City  West  geologic 
map  shows  a  set  of  parallel  faults  (shown  on  geologic  maps  by  both  contin¬ 
uous  and  dashed  lines)  trending  east-northeast  and  another  set  trending 
north-northeast . 

A  set  of  lineaments  that  may  be  faults  or  joints,  visible  on  radar 
but  not  shown  on  the  geologic  map,  trends  northwest-southeast. 
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It  is  very  significant  that  the  large  mines  (shown  on  radar  by  the 
large  waste  dumps)  tend  to  be  located  (1)  at  the  contact  of  a  calcareous 
formation  and  a  relatively  impermeable  formation  and  (2)  where  the  linea¬ 
ments  visible  on  radar  intersect  a  fault.  It  might  be  possible  to  locate 
new  mineral  deposits  with  a  knowledge  of  the  new  lineaments  shown  on 
this  radar  imagery. 

Ancestral  Drainage 

A  significant  feature  of  radar  imagery  is  its  ability  to  depict  former 
stream  courses.  The  imagery  shows  distinctly  the  former  courses  of  the 
tributaries  in  Upper  Provo  Canyon  (Figures  5  and  6). The  angles  at  which 
these  streams  join  the  Provo  River  form  acute  angles  that  point  east¬ 
ward,  or  up  the  present  course  of  the  Provo.  Through  almost  the  entire 
length  of  the  Upper  Provo  Canyon  the  streams  apparently  flowed  east  to 
Rhodes  Valley  and  northward  to  the  Weber  River.  Rhodes  Valley  is  a  grossly 
oversized  valley  for  Beaver  Creek  which  flows  through  its  central  portion. 
This  example  of  stream  piracy  has  been  known  since  Anderson  described 
it  in  1915  and  it  was  chosen  by  Longwell  and  Flint  in  their  physical  geology 
text  to  illustrate  stream  piracy,  but  it  is  so  striking  on  radar  imagery  that 

it  suggests  radar's  value  as  a  useful  tool  for  mapping  ancestral  stream 
patterns . 

The  radar  imagery  shows  a  well-developed  dendritic  stream  pattern 
in  the  headwaters  of  Mill  Creek.  The  present  stream  divide  transects 
the  upper  part  of  this  stream  pattern  which  results  in  an  abnormal 
appearing  drainage  pattern  east  of  the  headwaters  of  Mill  Creek.  It  is 
postulated  that  water  from  Parleys  Park  may  have  formerly  drained  west¬ 
ward  through  Mill  Creek  and  possibly  other  westward  flowing  streams. 

It  should  be  noted  that  the  streams  at  the  head  of  Mill  Creek  Canyon  flow 
at  high  angles  to  the  strike  of  the  Paleozoic  and  Mesozoic  sediments  which 
have  been  folded  and  thrust-faulted.  There  is  also  a  possibility,  suggested 
by  drainage  shown  on  radar  imagery,  that  other  westward-flowing  streams, 
such  as  Big  Cottonwood,  formerly  had  tributaries  that  extended  farther 
eastward  (Figures  5  and  6),  To  illustrate  the  comparative  value  of  radar 
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in  detecting  old  stream  courses  an  air  photo  mosaic  of  upper  Mill  Creek,  with 
a  scale  larger  than  the  scale  of  the  radar  imagery  was  examined;  the  ground 
clutter  on  the  photo  obscured  the  ancestral  drainage  pattern. 

There  has  been  much  speculation  as  to  why  the  Weber  and  Provo 
Rivers,  which  drain  the  eastern  Wasatch  and  the  western  Uintas,  flow 
across  the  Wasatch  Range.  One  explanation  is  Powell's  antecedence 
hypothesis  of  a  river  maintaining  its  course  while  cutting  through  a  slowly 
rising  mountain  range.  A  later  hypothesis  was  that  of  superposition, 
whereby  the  river  courses  were  superimposed  from  Tertiary  strata  which 
unconformably  covered  the  Laramide  structures .  However,  in  the  Wasatch 

Range  there  are  no  remnants  of  the  Tertiary  cover  near  the  cross  drainage 
of  the  Weber  and  Provo  Rivers . 

Hunt,  in  his  1956  report  on  the  Colorado  Plateau,  used  a  combina¬ 
tion  of  Powell's  antecedence  and  Gilbert's  superposition  in  what  he  called 
antepositlon. 

Hinze  (1913)  and  others  postulated  that  consequent  streams  on  the 
west  side  of  the  Wasatch  might  have  cut  across  the  range  to  capture  streams 
which  drained  to  such  lowlands  as  the  Uinta  or  Bridger  Basins.  According 
to  this  concept  rivers,  such  as  the  Weber  and  Provo,  that  captured  suffi¬ 
cient  drainage  were  able  to  cut  gorges  through  the  Wasatch  Range. 

Eardley  believes  that  the  Weber  and  Provo  Rivers  were  established 
by  Oligocene  time  and  that  they  maintained  their  courses  as  the  Wasatch 
Range  was  elevated  first  by  folding,  then  by  faulting.  He  believes  that 
the  Wasatch  Range  has  been  tilted  by  uplift  of  the  range  along  the  Wasatch 
Front.  An  absolute  upward  movement  on  the  east  side  of  the  Wasatch  Fault 
would  have  tilted  the  Wasatch  some  3*  or  4“  along  a  hingeline.  Eardley 
believes  an  axis  of  tilt  lies  at  the  head  of  the  Provo  Canyon;  that  below 

this  axis  the  Provo  has  cut  a  gorge  and  above  the  axis  the  Heber  Valley  is 
filled  with  alluvium. 

Threat  (1959)  questions  the  location  of  the  hinge  and  its  existence, 
because  (1)  Eardley's  axis  of  tilt  coincides  with  the  intersection  of  the 
present  profile  of  the  Provo  River  and  the  hard  Oquirrh  Quartzite  and 
(2)  it  is  not  known  that  Heber  Valley  is  deeply  alluviated. 

The  postulated  reversal  of  drainage  in  upper  Mill  Creek  Canyon 
supports  Eardley's  hypothesis  regarding  uplift  and  tilting  of  the 
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Wasatch  Range  and  the  location  of  the  hingeline.  It  is  postulated  that 
Mill  Creek  was  a  consequent  westward  flowing  stream,  that  at  that  time 
Parleys  Park  was  still  filled  by  the  Tertiary  conglomerates  and  andesites 
which  today  cover  hills  adjacent  to  Parleys  Park.  Uplift  of  the  Wasatch 
Range  caused  tilting  along  an  axis  which  is  located  at  or  near  the  present 
drainage  divide  between  Mill  Creek  and  east-flowing  streams .  This  loca¬ 
lity  Is  20  miles  north-northwest  of  Eardley's  Provo  Canyon  locality  and 

the  same  distance  east  from  the  Wasatch  Front. 

When  the  area  east  of  the  hinge  was  lowered  by  tilting ,  drainage 
in  the  headwaters  of  ancestral  Mill  Creek  was  reversed.  That  the  ances¬ 
tral  stream  pattern  has  been  altered  by  glaciation  indicates  that  the  tilting 
was  prePlelstocene.  Since  the  Ultlng,  according  to  this  hypothesis,  the 
level  of  Parleys  Park  has  been  lowered  by  erosion  of  Tertiary  conglomer¬ 
ates  and  andesites  which  formerly  must  have  filled  the  lowland. 

Landform  Classification 

Radar  imagery  appears  to  have  value  as  a  tool  in  studying  genetic 
landscape  evoluation  and  for  quantitative  geomorphic  classification. 
Although  much  research  reamains  to  be  done  in  this  area  findings  to  date 

are  encouraging. 

As  part  of  a  study  on  Integrated  landscapes  in  Puerto  Rico,  David 
Schwarz  and  Major  Roland  Mower  compared  landforms  discernable  on 
radar  Imagery  with  the  quantitative  landform  classification  of  the  Island 
by  Young  (1953).  Young's  classification  was  based  on  three  parameters 
represented  by  three  digits.  The  digits  show  local  relief,  percentage 
of  steep  land,  and  percentage  of  flat  land.  With  this  criteria  for  classi¬ 
fication,  shown  in  Table  1 ,  Young  devised  the  classification  of  major 
landforms  shown  in  Table  2.  In  many  respects  Young's  classification 

resembles  that  of  Hammond  (1964) . 

Young's  landform  classification  was  developed  specifically  for 
Puerto  Rico  and  was  based  on  topographic  maps.  In  their  study  of  Puerto 
Rico  with  classified  -adar  imagery,  Schwarz  and  Mower  superposed  Young's 
landform  regions  on  the  radar  imagery.  They  found  that  they  could  Identify 
Young's  landform  regions  on  the  radar  imagery  and  in  some  cases  were 
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able  to  add  detail  to  them.  For  example,  in  an  area  described  by  Young 
as  rugged  hill  land  (classified  233)  with  one  corner  of  the  region  described 
as  rolling  hills  (213),  Schwaiz  and  Mower  distinguished  two  different  types 
of  topography,  rather  than  one,  in  the  larger  part  of  the  region. 

They  also  found  that  Young's  classification,  which  was  based  on 
topographic  maps,  was  generalized  and  in  many  places  missed  much  detail 
shown  on  radar.  For  example,  in  Schwarz  and  Mower's  regions  C-2,  C-3, 
and  C-4,  the  imagery  shows  many  lowland  areas  between  the  rugged  hills 
and  the  rolling  hills  which  were  not  included  in  Young's  classification. 

Other  examples  can  be  seen  in  this  part  of  Puerto  Rico  where  the  geomor- 
phic  classification  based  on  topographic  maps  generalized  much  more  than 
a  classification  based  on  radar  imagery. 

Schwarz  and  Mower  concluded  that  radar  imagery  could  be  used  to 
refine  and  add  detail  to  Young's  classification  and  that  radar  imagery  is 
adequate  for  some  aspects  of  regional  generalization.  Discrimination 
between  landform  regions  such  as  flat  lowlands  and  low  hills  could  be 
accomplished  with  relative  ease.  They  did  experience  difficulty  in  differ¬ 
entiating  between  various  degrees  of  slopes  in  mountainous  areas. 

Included  in  this  paper  (Fig .  7)  is  a  landform  map  of  part  of  the< 
Wasatch  Range  using  Hammond's  1964  classification,  for  which  criteria 
are  given  in  Table  3,  Hammonds  classification  is  included  here  in  Table 

Although  some  difficulty  was  experienced  in  adapting  Hammonds' 

« 

broad  classification  to  the  small  study  region,  with  modifications  the 
classification  was  found  to  be  workable. 

Hammond  based  his  classification  on  the  three  parameters  of  slope, 
loccT  relief,  and  profile  type.  The  symbols  B3b,  according  to  Hammond's 
classification,  describe  an  area  which  is  50-80%  gently  sloping  (B);  has 
300-500  feet  of  local  relief  (3);  and  50-75%  of  the  gentle  slope  is  in  low¬ 
land  (b) .  Hammond,  who  considered  areas  of  about  six  square  miles  as 
the  smallest  units  to  be  classified,  stated  that  since  his  map  is  small 
scale  and  is  based  upon  a  simple  classification,  it  is  a  highly  abstract 
version  of  reality. 

Because  the  Wasatch  Range  map  is  drawn  to  a  much  larger  scale 
than  Hammond's  maps,  some  modifications  in  the  classifications  were  made. 
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The  local  relief  and  profile  type  classes  remain  the  same,  but  the  slope 
classes  were  changed.  Hammond  selected  the  value  of  eight  per  cent 
inclination  as  the  upper  limit  of  "gentle"  slope,  but  for  the  Wasatch  Range, 
where  little  ot  the  area  is  gently  sloping,  slope  was  modified  to  refer  to 
relative  slope  rather  than  to  any  absolute  slope  gradient.  Therefore,  in 
this  context,  C  slopes  are  steeper  than  B  slopes  and  not  as  steep  as  D 
slopes;  areas  characterized  as  B  generally  have  gentle  slopes  and  D 
values  have  extreruely  steep  slopes. 

Physiographic  Maps  from  Radar 

Radar  imagery  portrays  terrain  features  in  such  a  way  that  the 

imagery  is  a  generalized,  shaded,  relief  map.  Not  only  is  there  no  sym¬ 
bolization  problem,  as  on  topographic  and  physiographic  maps,  but  there 
are  continuous  graytone  values  on  the  radar  imagery  to  convey  a  fairly 
detailed  representation  of  terrain  over  large  areas .  With  the  addition  of 
some  place  names  a  print  of  radar  imagery  portrays  general  physiography. 

The  two  physiographic  maps  in  this  paper  (Figures  8  and  9) 
portray  the  same  area  in  the  wasatch  Range.  Figure  8  was  drawn  by 
William  Brooner  who  used  a  quill  pen  and  India  ink  to  sketch  topography 

on  tracing  paper  laid  over  the  radar  imagery. 

Figure  9 ,  was  drawn  by  Mrs .  William  Brooner  from  a  2X  enlarge¬ 
ment  of  radar  Imagery.  Various  graytones  were  achieved  by  using  a  brush 
and  different  dilutions  of  India  ink.  Some  features,  such  as  ridge  crests, 

were  drawn  with  pen  and  ink.  Although  this  illustration  incorporates 

details  from  the  radar  Imagery,  it  is  more  of  an  artist's  concept  than 
Figure  8. 

It  is  necessary  to  emphasize  that  both  of  these  physiographic 
diagrams  are  first  attempts  by  persons  with  only  limited  cartographic 
experience  to  show  that  the  technique  of  tracing  generalized  terrain  maps 
from  radar  is  feasible  and  may  have  potential.  It  may  be  that  generalized 
terrain  maps  could  be  produced  quicker  and  for  less  cost  from  radar,  or 
a  combination  of  radar  and  space-borne  photography  ,  than  from  large  scale 
aerial  photographs  and  topographic  maps . 
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TABLE  1 

Criteria  for  Young’s  Classification  of  Landforms 

First  Digit:  local  relief  in  meters 
1:  0  -  112 
2:  112  -  242 
3:  242  -  770 

Second  Digit:  percent  of  steep  slope* 

1:  1  -  5 
2:  5-25 
3:  25  -  100 

Third  Digit:  percent  of  flat  land"*" 

1:  100  -  60 
2:  60  -  3 
3:  3  -  0 


*  Steepland  -  slopes  greater  than  59% 
+  Flatland  -  slopes  less  than  3% 
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TABLE  2 

Major  Landform  Classification  by  Young 

LOWLANDS 

(111)  Flat  lowlands 

(112)  Rolling  lowlands  with  some  flat  land 

(113)  Rolling  lowlands 

(122)  Rough  lowlands  with  some  flat  land 

HILL  LANDS 

(212)  Rolling  hill  land  with  some  flat  land 

(213)  Rolling  hill  lands 
(223)  Rough  hill  lands 

(232)  Rugged  hill  lands  with  some  flat  land 

(233)  Rugged  hill  lands. 

MOUNTAIN  LANDS 

(313)  Rolling  mountain  lands 
(323)  Rough  mountain  lands 
(333)  Rugged  mountain  lands 


TABLE  3 

Criteria  for  Hammond*s  Classification  of  Landforms 


First  Letter:  percent  of  area  gently  sloping* 
A  >80 
B  50-80 
C  20-50 
D  <20 

Second  letter:  local  relief  in  feet 
1:  0  -  100 
2;  100  -  300 
3:  300  -  500 
4:  500  -  1000 
5:  1000  -  3000 
6:  3000  -  5000 

Third  letter:  profile  ty.ie  -  percent  of  gentle 
slope  in  lowland  or  on  upland 

a  >  75  (in  lowland) 

b  50  -  75  (in  lowland) 

c  50  -  75  (on  upland) 

d  <  75  (on  upland) 


*  Gentle  slope  -  Q%  or  less  slope 


TABLE  4 

Major  Landform  Classification  by  Hammond 


PLAINS 

A1  Flat  plains 
A2  Smooth  plains 
B1  Irregular  plains,  slight  relief 
B2  Irregular  plains 

TABLE  LANDS 

B3c,d  Tablelands,  moderate  relief 
B4c,d  Tablelands,  considerable  relief 
B5c,d  Tablelands ,  high  relief 
B6c,d  Tablelands,  very  high  relief 

OPEN  HILL  AND  MOUNTAINS 
02  Open  lo  v  hills 
03  Open  hills 
04  Open  high  hills 
05  Open  low  mountains 
06  Open  high  mountains 

HILLS  AND  MOUNTAINS 
D3  Hills 
D4  High  hills 
D5  Low  mountains 
D6  High  mountains 
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TABLE  4 

Major  Landform  Classification  by  Hammond 


PLAINS 

A1  Flat  plains 
A2  Smooth  plains 
B1  Irregular  plains,  slight  relief 
B2  Irregular  plains 

TABLE  LANDS 

B3c,d  Tablelands ,  moderate  relief 
B4c,d  Tablelands,  considerable  relief 
B5c,d  Tablelands,  high  relief 
B6c,d  Tablelands,  very  high  relief 

OPEN  HILL  AND  MOUNTAINS 
02  Open  low  hills 
03  Open  hills 
04  Open  high  hills 
05  Open  low  mountains 
06  Open  high  mountains 

HILLS  AND  MOUNTAINS 
D3  Hills 
D4  High  hills 
D5  Low  mountains 
D6  High  mountains 


Summary 


The  observations  in  this  report  are  preliminary.  Much  work  re¬ 
mains  to  be  done  in  the  way  of  regional  description  and  analysis  which 
will  be  followed  by  field  checks. 

However,  this  report  indicates  that  radar  imagery  has  utility’  for 
regional  geomorphic  analysis  and  description  and  for  land  use  studies. 

The  lack  of  ground  clutter  on  radar  allows  one  to  see  geomorphic  features 
that  are  often  obscured  on  aerial  photographs . 

The  topography  of  the  Wasatch  Range  is  well  displayed  on  radar 
imagery.  Especially  striking  are  glacial  features,  some  fault  scraps,  and 
stream  patterns.  Radar  shows  unmapped  lineaments,  possibly  representing 
faults,  that  control  development  of  certain  landforms  and  apparently 
deposition  of  mineral  deposits.  A  postulated  ancestral  stream  pattern  with 
reversal  of  drainage  in  the  headwaters  seems  to  yield  fundamental  data 
on  the  origin  of  cross  cutting  rivers  and  tilting  in  the  Wasatch  Range. 

The  8  thematic  maps  prepared  from  radar  imagery  used  in  this  study 
are:  land  use,  transportation,  surface  materials,  landform  classification, 
present  drainage,  ancestral  drainage,  and  2  physiographic  maps.  In  the 
next  phase  of  research  other  maps,  such  as  vegetation  maps,  will  be 
prepared,  regional  descriptions  will  be  expanded,  and  field  checks  will- 
be  performed. 
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Figure  7.  Landform  classification  map  of  the  Wasatch  Range,  Utah.  Derived  from  K-band  radar 
imagery  except  for  elevations  . 


Figure  8.  Physiographic  Map  of  the  Wasatch  Range,  Utah. 


Picture  9.  Physiographic  Map  of  the  Wasatch  Range,  Utah. 
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